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MEASURING  THE  COMPONENTS  OF  THE  EARTH'S  GEOMAGNETIC  FIELD  BY 
THE  PROTON  RESONANCE  METHOD 
V.  D.  Loraanyy,  A.  G.  Prokofyev,  and  B.  M.  Yanovskiy 
PP.  3-15 


Introduction 

All  magnetic  observatories  in  the  USSR  at  present  use 
classical  methods  to  measure  H  (the  horizontal  component)  and 
Z  (the  vertical  component).  The  error  in  H  as  so  measured  is 
2-3  Y«  while  that  in  Z  for  latitudes  between  50  and  60*  is  5- 
10  y.  One  of  the  proton-resonance  methods  (free  nuclear  in¬ 
duction)  can  measure  the  field  to  1  part  in  10“*,  which  corre¬ 
sponds  to  not  more  than  0.5  y  at  the  above  latitudes;  this  re¬ 
presents  a  great  improvement.  The  method  gives  the  strength 
of  the  total  field,  so  either  component  can  be  measured  alone 
only  if  the  other  is  first  balanced  out  by  a  field  of  compar¬ 
able  homogeneity.  Helmholtz  coils  can  provide  this  if  the 
homogeneous-field  region  is  larger  than  the  detector  (a  vessel 
containing  water).  Here  we  describe  an  apparatus  for  measur¬ 
ing  Z  and  give  some  values. 

1,  Theoretical  Principles  of  the  Free  Nuclear 
Induction  Method 

Any  atomic  nucleus  has  a  mechanical  moment  p  and  a  mag¬ 
netic  moment  p;  p/p  =  ty  is  called  the  gyromagnetic  ratio,  the 
sign  indicating  the  direction  of  p  relative  to  p.  The  protons 
in  -(say)  water  have  no  preferred  orientation  in  the  absence  of 
an  external  field,  but  a  field  H  (steady)  gives  then  an  addi¬ 
tional  energy  AE  a  pH,  which  causes  p  to  take  the  direction  of 
H.  Quantum  theory  indicates  that  the  projection  of  p  on  the 
direction  of  H  can  take  2m  +  1  values,  in  which  m  is  the  magne¬ 
tic  quantum  number;  m  =  %  for  the  proton,  so  it  has  two  orien¬ 
tations  (along  and  opposed  to  the  field).  Boltzmann's  theory 
indicates  that  the  proportion  parallel  to  the  field  is  greater 
than  that  antiparallel,  so  there  is  an  excess  moment  p^  in  the 
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direction  of  H  (the  substance  is  paramagnetic1) .  The  Magnetic 
moment  per  unit  volume  is  J0  =  xH,  in  which  x  i ,s  the  nuclear 
paramagnetic  susceptibility: 

Afr*  I) 

x_  . 


in  which  N  is  the  number  of  nuclei  in  unit  volume,  k  is  Boltz¬ 
mann's  constant,  and  T  is  absolute  temperature. 

The  interaction  between  H  and  p  is  accompanied  by  inter¬ 
actions  between  the  nuclei  individually  and  collectively.  The 
latter  (interaction  with  the  lattice)  causes  the  temperature  t' 
of  the  nuclear-spin  system  to  be  the  temperature  t"  of  the  latt¬ 
ice  in  the  usual  state;  but  t. 1  becomes  larger  than  t"  when  a 
magnetic  field  is  applied,  and  a  certain  time  T,  (the  longitud¬ 
inal  relaxation  time)  is  needed  for  the  system  to  return  to 
thermodynamic  equilibrium.  This  means  that  the  vector  does 
not  rise  to  JQ  =  xH  instantly  but  instead  rises  in  accordance 
with  J  =  J0[l  -  exp(-t/T)3,  in  which  t  is  time. 

The  interaction  between  the  spins  is  caused  by  the  local 

magnetic  field  H'  produced  by  adja¬ 
cent  nuclei,  so  spins  initially  mo¬ 
ving  in  phase  in  response  to  H  after 
a  time  T^  differ  by  Tp  oC 

l/Aui  and  is  called  the  transverse 
relaxation  time.  The  commonest  me¬ 
thod  of  exciting  free  nuclear  in¬ 
duction  is  to  switch  of  an  accessory 
strong  magnetic  field  H  very  rapid¬ 
ly,  this  polarizing  field  is  first 
applied  perpendicular  to  Ht ,  the 
field  to  be  measured.  Then,  after 
a  time  of  about  Tj,  ,  when  "5  is  app¬ 
roaching  its  equilibrium  value, 
this  H  is  removed  rapidly  (so  rapid¬ 
ly  that  J  does  not  alter  in  orien¬ 
tation).  The  oriented  protons  are 
then  subject^  to  Ht,_>the  couple  on 
the  moment  u.  from  H,.  being  C  = 

L U ^ * 

The  equation  of  motion,  in 
the  absence  of  relaxation,  is 


Fig.  1. 


for  Mp  =  Ypi' 


~HT  =  C~  [pi/ir  ]  =  -  (PjHr  |,  CD 

V/e  multiply  both  sides  scalarly  by  and  get 
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Pi(dPi/dt)  =  0,  so  P  lies  at  90*  to  dP,  the  modulus  of  P  being 
constant,  so  the  mocion  of  f  is  that  of  precession  about  the  di¬ 
rection  of  Hf.. 

Let  ♦  b£  the  angle  between  P  and  (Fig.  1);  then  dP^ 
can  be  put  as  P^sin  so 

~jij-  «  Pi  sin  ®  ™  sin  I, 

in  which  ut  is  the  angular  velocity  of  P.  We  replace  dPi/dt 
here  to  get 


i^//riin  I 

%it  la  I 


-£-//r=T//r, 

Pi 


in  which  U)  is  the  Larmor  precession  frequency. 

Equation  (1)  should  be  put  as  follows  if  relaxation 
effects  are  present,  for  these  resemble  damping;  replacing 
by  J,  we  have 


The  conditions  for  t  =  0  are 

/,=o,s^=#=o, 

and  the  solution  to  the  equations  is 

_  t 

Jx  =  J¥e  T‘  cos 

_  t 

Jy  =  J$e  7,sln(»f-f-<p0), 

2-n  which  to  =  yH|.  and  (fQ  is  the  initial  phase.  Then  J  precesses 
around  Ht  with  a  frequency  cu  and  has  a  decreasing  amplitude;  a 
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cylindrical  coil  having  its  axis  in  the  plane  of  rotation  then 
has  a  current  induced  in  it.  The  amplitude  of  mis  signal  doe 
not  depend  on  the  orientation  of  the  coil  in  the  xCy  plane  but 
it  doe3  depend  on  the  inclination  in  the  zOy  plane,  for  this 
affects  the  effective  area  intersecting  the  flux.  The  signal 
decays  gradually,  but  its  amplitude  remains  adequate  for  seve¬ 
ral  seconds  provided  that  the  relaxation  time  is  long. 

The  magnetic  induction  Bx  corresponding  to  Jx  is  4n«Jx, 
so  the  flux  linkage  ^  is  BxSN(3sin(a) ,  in  which  S  is  the  effec¬ 
tive  area  of  one  turn,  a  is  the  angle  between  Jx  and  H^,  B  is 
the  filling  factor,  and  N  is  the  number  of  turns.  We  assume 
that  sin  a  =  1  and  p  =  1,  so  the  induced  emf  is 

£ =  -  as  ~  NJ,m  lin  *>t. 


The  power  corresponding  to  this  signal  is 

p=t(£~>  yr)’ "  <» 

in  which  r  is  the  circuit  resistance.  In  our  case  the  para¬ 
meters  were  r  =  I85/I,  N  =  2550,  S  =  6.2  cm^,  and  H  =  1100  oe^ 
These  give  with  (2)  that  P  =  2.1  x  10--^  w  and  that  E  =  (2rP)'‘‘ 

=  9  pv,  so  the  signal  is  rather  weak;  even  an  H  of  1100  oe 
gives  a  signal  comparable  with  the  noise,  but  narrow-band  fil¬ 
ters  eliminate  most  of  this.  The  signal-to-noise  ratio  is 
about  7  or  ?  in  the  present  system. 

2.  The  Apparatus 

The  apparatus  for  detecting  and  measuring  the  frequency 
given  by  Z  has  two  main  parts,  namely  the  Helmholtz  coils  with 
the  sample  and  the  electronic  equipment  for  measuring  the  frequ 
ency.  Figure  2  shows  the  block  diagram;  the  Helmholtz  coils 
1  balance  out  the  horizontal  component.  The  sample  (distilled 
water)  4  lies  at  the  center  of  the  coils  and  is  enclosed  in  the 
pickup  coil  5;  this  last  is  itself  enclosed  in  the  polarizing 
coil  2.  The  emf  in  coil  3  has  the  frequency  u);  it  is  ampli¬ 
fied  by  the  preamplifier  5  and  then  is  doubled  in  frequency  by 
the  multiplier  6,  after  which  it  is  again  amplified  and  is  fed 
to  the  mixer  7,  where  it  is  mixed  with  a  signal  of  frequency  ufc 
from  the  calibrated  generator  8.  The  two  frequencies  are  si¬ 
milar,  and  the  beat  frequency  =  4>0  -  u>  is  detected  by  the  de 
tector  9i  this  is  passed  through  the  matching  stage  10  to  the 
oscillographs  11  and  12.  The  MPO-2  oscillograph  records  the 
beat  frequency  on  film,  together  with  time  marks  provided  by 
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the  generator. 


Fig.  2. 

The  Helmholtz  coils  (Fig.  3)  consist  of  two  flat  discs  1 
of  plexiglass  360  mm  in  diameter  each  bearing  112  turns.  These 
discs  are  fixed  to  the  plexiglass  cube  2  (edge  180  nun),  the 
cube  being  supported  on  the  cylinder  3  mounted  on  the  goniometer 
stage  4.  The  leveling  screws  5  are  used  to  set  the  axis  of  the 
coils  horizontal,  the  stage  providing  a  means  of  setting  in  the 
correct  direction.  The  cube  contains  the  polarizing  coil, 
within  which  there  is  the  pickup  coil,  which  is  wound  on  the 
tube  containing  the  distilled  water.  The  polarizing  coil  has 
its  axis  normal  to  the  axis  of  the  Helmholtz  coils  and  to  the 
axis  of  the  pickup  coil.  The  above  dimensions  give  the  Helm¬ 
holtz  coils  a  field  constant  of  2.8  oe/a,  so  the  horizontal 
component  at  Leningrad  (0.15  oe)  is  balanced  by  a  current  of 
about  55  ma. 

The  polarizing  coil  is  a  straight  solenoid  lying  in  a 
horizontal  plane;  it  has  2200  turns  of  copper  wire  (diameter  1 
mm)  and  its  constant  is  125  oe/a. 

The  pickup  coil  consists  of  2 550  turns  (diameter  0.26  mm) 
and  forms  part  of  a  tuned  circuit  in  the  amplifier;  it  is  tuned 
to  2075  c/s  and  its  3  is  15-17. 

The  specimen  (distilled  water)  is  contained  in  a  sealed 
plastic  cube  of  side  2.5  cm  (volume  16  cc).  This  size  corre¬ 
sponds  to  the  region  of  homogeneity.  Four  of  the  sides  act  as 
former  for  the  pickup  coil,  so  p  «  1. 

The  signal  produced  by  J  is  amplified  by  system  giving  a 
good  signal-to-noise  ratio  that  is  stable  and  not  sensitive  to 
external  interference.  The  four-tube  system  has  an  over-all 
gain  of  5  x  10';  the  first  tube  (6ZhlP)  is  triode-connected ,  to 
minimize  noise  in  this  stage.  The  main  source  of  noise  is  the 
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fluctuations  in  t]ie  coil,  though;  Nyquist's  formula  gives  this 
noise  voltage  as  U“  =  4kTRAf,  in  which  k  is  Boltzmann's  constant. 
T  is  temperature  (*K),  4f  io  bandwidth,  and  3  i. ",  resistance. 


Fig. 

Clearly,  one  way  of  minimizing  the  noise  is  to  make  the  ampli¬ 
fier  of  small  Af,  but  too  small  a  value  is  undesirable,  for  the 
permissible  range  of  variation  in  4  is  then  small.  The  A f 
chosen  was  ±  60  c/ s  for  the  above  central  frequency  of  2075 
c/s;  thenyu  is  about  1  uv.  The  contribution  from  the  oZhlP 
make?  the  total  only  1.5  mv  at  most,  /sich  means  a  signal-to- 
noise  ratio  of  7-8. 

The  high  gain*nakes  a  severe  demand  on  the  screening, 
for  any  interference  enterin'-  in  the  early  stages  may  distort 
or  even  mask  the  signal.  All  the  input  circuits  were  care¬ 
fully  screened,  each  stage  being  contained  in  a  separate  screen¬ 
ed  box,  the  tubes  and  variable  capacitors  being  enclosed  in 
special  screens.  Screened  cables  and  couplers  were  used 

•aoplific  atio*  coefficient 
of  5  •  105 
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throughout,  all  screens  bein g  grounded.  The  poi -.rizing  coil 
was  short-circuited  and  grounded  when  not  carrying  current;  it 
then  acted  as  a  screen  for  the  pickup  coil. 

The  frequeiicy-converting  unit  includes  the  doubler,  the 
amplifying  filter,  the  mixer,  and  a  cathode  follower.  The 
doubler  reduces  Afc/fc  for  a  given  Afc  (error  of  frequency  mea¬ 
surement);  it  consists  of  a  6ZhlP  whose  plate  circuit  contains 
the  primary  of  a  transformer,  whose  second  is  connected  to  a 
full-wave  rectifier  (DGTs-23  diodes).  The  signal  at  the  out¬ 
put  of  the  doubler  is  reduced  in  amplitude  and  Bomewhat  distor¬ 
ted,  so  2fc  is  isolated  and  amplified  by  a  stage  containing  a 
6N9S.  The  mixer  compares  2fc  with  the  standard  frequency  (from 
a  quartz  crystal)  of  4166  c/s;  at  the  plate  of  the  mixer  there 
appears  the  beat  (difference)  frequency  ffc.  A  detector  con¬ 
taining  a  DGTs-27  diode  is  used  to  isolate  f^.  The  cathode 
follower  feeds  the  loop  oscillograph  with  the  beat  frequency, 
the  bias  for  the  tube  being  produced  by  a  resistance  of  200/i, 
only  the  alternating  component  reaching  the  MFO-2  oscillograph 
(via  a  20  pF  capacitor).  The  quartz  oscillator  is  designed 
to  maintain  a  fixed  frequency  within  0.0001%,  the  fundamental 
standard  being  a  100  kc  crystal,  whose  frequency  remains  con¬ 
stant  to  this  order  after  it  has  been  running  for  two  hours. 

The  frequency  of  4166.6  c/s  is  generated  by  successive  division 
by  4  and  6  from  100  kc,  the  filtered  and  amplified  product 
being  fed  to  the  mixer. 

The  100  kc  quartz  crystal  is  enclosed  in  a  thermostatic 
vessel;  it  forms  the  first  stage  in  a  self-excited  oscillator, 
where  it  replaces  the  usual  tuned  circuit  in  the  grid  circuit. 
This  stage  uses  half  a  6N1F.  The  plate  circuit  is  tuned  to 
120  kc  to  improve  the  stability  and  is  also  decoupled.  The  po¬ 
sitive  half-cycles  are  passed  to  the  4:1  division  stage  via  a 
circuit  consisting  of  n  DGTs-25  diode,  a  capacitor,  e.nd  a  resis¬ 
tor. 


3.  Results 

The  horizontal  component  H  is  balanced  out  by  setting 
the  axis  of  the  Helmholtz  coils  in  the  plane  of  the  magnetic  me¬ 
ridian  and  in  the  horizontal  plane.  Any  such  setting  is  made 
with  a  certain  error,  so  we  must  consider  the  effect  on  the  mea¬ 
sured  Z  of  any  error  in  balancing  H. 

Let  the  coil  make  an  angle  a  with  the  horizontal  plane 
(Fig.  4)  and  an  angle  180  -  3  with  the  plane  of  the  magnetic  me¬ 
ridian;  let  the  actual  field  of  the  coil  be  ,  which  differs 
from  H  by  AH  «  H  -  Hj,  Then  the  residual  field  at  the  center 
is  given  by 

Z'*  =  (i?  +  //,  sin  a)*  +  (H  —  Hx  cos*  •cos?)1'  +  //Jsln1?, 
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which,  if  we  neglect  qu  anti ti of  the  fourth  or  .e  >•  in  dH,  a, 
and  P  (which  are  all  snail),  we  have 


Z' —  Z  =  HX  sin  *  + 


T+fzV  +  P'). 


The  largest  terra  is  the  first,  but  this  is  re  ■••■lily  eliminated 
by  turning  the  coil  through  lnO*  and  taking  a  second  reading 
Z";  a  and  p  then  change  sign,  so  H.sin  c.  does  so  as  well. 
Then  Z'  -  Z  and  Z"  -  Z  give  us  the  relative  error  £Z/Z  as 
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I  (SH?  I 

7  ~z*~  +  7?r  (**  +  ?*)• 


For  example,  if  the  error  in  balancing  H  is  it'd/ H  =  0.01 


Fig.  k. 

and  if  a  =  P  =  1*,  we  have  dZ/Z  -  0.  jCC>5 ,  of  which  1/7  ‘rises 
from  JH  and  6/7  fro:?:  a  find  p. 

The  H/Z  found  in  the  USS.-i  range  fro:.:  0.001  (in  the  polar 
regions)  to  0.75,  so  the  dZ/h  for  the  above  conditions  range 
from  0  to  O.0CO18.  The  error  from  all  causes  can  be  kept  be¬ 
low  0.01?*  if  H  is  balanced  out  to  within  a  few  ?.  and  the  angles 
are  correct  within  1*. 

The  current  required  by  the  coils  at  their  present  site 
is  about  55  ma;  the  exact  value  may  be  determined  by  experi¬ 
ment. 

Let  the  current  be  I  for  exact  compensation,  which  gives 
a  beat  frequency  f.  .  Then  currents  1^  =  I  -  41  and  IP  =  I  + 

41  give  Z'  =  Z  +  4Z  and  f  =  f ,  +  df,  and  so  on;  this  enables 
one  to  find  I,  and  such  that  f^  =  0,  which  is  easily  obser¬ 
ved.  Then  the  correct  current  is 

!  =  (!]_  +  I  ?)/!’. 


8 


Our  measurements  gave  I  as  53«5  ma,  -which  agrees  well 
with  the  value  calculated  above;  the  precession  signal  was  re¬ 
corded  at  intervals  of  1  ran  from  40  to  70  mu,  -uni  ojso  at  inter¬ 
vals  of  0.5  n.-j  near  the  point  of  compensation. 


Fig.  5. 

Measurements  of  Z  were  made  at  the  Kavgolovo  laboratory 
of  the  All-Union  Metrology  Research  Institute,  which  is  outside 
the  city,  on  15  May  I960;  recordings  were  made  of  the  beat  fre¬ 
quency  and  also  of  time  marks  separated  by  0.01  sec  (Fig.  5  is 
typical).  The  19  pairs  of  measurements  are  listed  in  tne  table. 
Readings 'were  taken  simultaneously  fro  an  M-14  magnetometer; 
recordings  were  made  of  trie  variations  in  K  with  a  variation  in¬ 
strument.  The  Z-vnriometer  provided  corrections  and  so  enabled 
us  to  establish  the  error  of  the  method;  the  H  variations  pro¬ 
vided  a  check  on  the  deviations  from  exact  compensation.  One 
scale  division  on  the  M-14  is  equivalent  to  10  y,  and  the  read¬ 
ings  can  be  made  to  within  0.1  division,  so  the  correction  for 
variation  can  be  made  to  1  y.  The  H-variometer  recording 
showed  no  major  changes  during  the  experiment,  and  the  compen¬ 
sation  remained  within  the  permitted  limits. 

The  table  shows  that  the  standard  deviation  is 

»  * 

in  which  n  is  the  number  of  observations ;  then  f v,  =  14.151  - 
0.016  c/s.  The  compensation-current  tests  showed  that  f0  > 

2fc,  for  zero  beat  was  obtained  for  Ip  and  1 2  that  "/ere  respec¬ 
tively  larger  and  smaller  than  Ip  so  f0  =  2f(Z  +  42),  in  which 
f(Z  +  42)  is  the  precesaion  frequency  for  currents  Ip  and  I^. 
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>♦  ■*  4  4  -* 

Pi(dPt/dt)  ■  0,  so  F  lias  at  90*  to  dP,  the  modulus  of  P  being 
oonstanti  so  the  motion  of  P  is  that  of  precession  about  the  di¬ 
rection  of  Ht. 

Let  +  ba  the  angle  between  &  and  ift  (Fig.  1);  then  dP^ 
can  be  put  as  P^sin  so 

■■  sin  I, 

in  which  w  is  the  angular  velocity  of  P.  We  replace  dPj/dt 
here  to  get 


in  which  tu  is  the  Larmor  precession  frequency. 

Equation  (1)  should  be  put  as  follows  if  relaxation 
effects  sure  present,  for  these  resemble  damping;  replacing 
by  J,  we  have 


The  conditions  for  t  s  0  are 

and  the  solution  to  the  equations  is 

_  t 

Jx-J^  *  CO.  (.*  +  *), 

_  t 

y,=v’T«in  (•<+!,), 

— /.(i  — .  *0 , 

in  which  U)  =  yHt  and  Vo  is  the  initial  phase.  Then  precesses 
around  Ht  with  a  frequency  w  and  has  a  decreasing  amplitude;  a 
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cylindrical  coil  having  its  axis  in  the  plane  of  rotation  then 
has  a  current  induced  in  it.  The  amplitude  of  this  signal  does 
not  depend  on  the  orientation  of  the  coil  in  the  xOy  plane  but 
it  does  depend  on  the  inclination  in  the  zOy  plane,  for  this 
affects  the  effective  area  intersecting  the  flux.  The  signal 
decays  gradually,  but  its  amplitude  remains  adequate  for  seve¬ 
ral  seconds  provided  that  the  relaxation  time  is  long. 

The  magnetic  induction  IL  corresponding  to  Jx  is  4*JX* 
so  the  flux  linkage  $  is  BxSN(3sin(a) ,  in  which  S  is  the  effec¬ 
tive  area  of  one  turn,  a  is  the  angle  between  Jx  and  H^,  P  is 
the  filling  factor,  and  N  is  the  number  of  turns.  We  assume 
that  sin  a  =  1  and  p  =  1,  so  the  induced  emf  is 

£=--l-£=±L  fin  mt. 


The  power  corresponding  to  this  signal  is 

in  which  r  is  the  circuit  resistance.  In  our  case  the  para¬ 
meters  were  r  =  185^*  N  =  2550,  S  =  6.2  cm2,  and  H  =  1100  oe, 
These  give  with  (2)  that  P  =  2.1  X  10-1^  w  and  that  E  =  (2rP)^‘ 

=  9  mv,  so  the  signal  is  rather  weak;  even  an  H  of  1100  oe 
gives  a  signal  comparable  with  the  noise,  but  narrow-band  fil¬ 
ters  eliminate  most  of  this.  The  signal-to-noise  ratio  is 
about  7  or  3  in  the  present  system. 

2.  The  Apparatus 

The  apparatus  for  detecting  and  measuring  the  frequency 
given  by  Z  has  two  main  parts,  namely  the  Helmholtz  coils  with 
the  sample  and  the  electronic  equipment  for  measuring  the  frequ¬ 
ency.  Figure  2  shows  the  block  diagram;  the  Helmholtz  coils 
1  balance  out  the  horizontal  component.  The  sample  (distilled 
water)  4  lies  at  the  center  of  the  coils  and  is  enclosed  in  the 
pickup  coil  5;  this  last  is  itself  enclosed  in  the  polarizing 
coil  2.  The  emf  in  coil  3  has  the  frequency  UJ;  it  is  ampli¬ 
fied  by  the  preamplifier  5  and  then  is  doubled  in  frequency  by 
the  multiplier  6,  after  which  it  is  again  amplified  and  is  fed 
to  the  mixer  7,  where  it  is  mixed  with  a  signal  of  frequency 
from  the  calibrated  generator  8.  The  two  frequencies  are  si¬ 
milar,  and  the  beat  frequency  =  <n0  -  tf  is  detected  by  the  de¬ 
tector  9?  this  is  passed  through  the  matching  stage  10  to  the 
oscillographs  11  and  12.  The  MPO-2  oscillograph  records  the 
beat  frequency  on  film,  together  with  time  marks  provided  by 
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the  generator. 


Fig.  2. 

The  Helmholtz  coils  (Fig.  3)  consist  of  two  flat  discs  1 
of  plexiglass  3&0  nun  in  diameter  each  bearing  112  turns.  These 
discs  are  fixed  to  the  plexiglass  cube  2  (edge  180  nun),  the 
cube  being  supported  on  the  cylinder  3  mounted  on  the  goniometer 
stage  4.  The  leveling  screws  5  are  used  to  set  the  axis  of  the 
coils  horizontal,  the  stage  providing  a  means  of  setting  in  the 
correct  direction.  The  cube  contains  the  polarizing  coil, 
within  which  there  is  the  pickup  coil,  which  is  wound  on  the 
tube  containing  the  distilled  water.  The  polarizing  coil  has 
its  axis  normal  to  the  axis  of  the  Helmholtz  coils  and  to  the 
axis  of  the  pickup  coil.  The  above  dimensions  give  the  Helm¬ 
holtz  coils  a  field  constant  of  2.8  oe/a,  so  the  horizontal 
component  at  Leningrad  (C. 15  oe)  is  balanced  by  a  current  of 
about  55  ma. 

The  polarizing  coil  is  a  straight  solenoid  lying  in  a 
horizontal  plane;  it  has  2200  turns  of  copper  wire  (diameter  1 
mm)  and  its  constant  is  125  oe/a. 

The  pickup  coil  consists  of  2550  turns  (diameter  0.26  mm) 
and  forms  part  of  a  tuned  circuit  in  the  amplifier;  it  is  tuned 
to  2075  c/s  and  its  3  is  15-17. 

The  specimen  (distilled  water)  is  contained  in  a  sealed 
plastic  cube  of  side  2.5  cm  (volume  16  cc).  This  size  corre¬ 
sponds  to  the  region  of  homogeneity.  Four  of  the  sides  act  as 

former  for  the  pickup  coil,  so  1. 

The  signal  produced  by  J  is  amplified  by  system  giving  a 
good  signal-to-noise  ratio  that  is  stable  and  not  sensitive  to 
external  interference.  The  four-tube  system  has  an  over-all 
gain  of  5  x  10*;  the  first  tube  (6ZhlP)  is  triode-connected ,  to 
minimize  noise  in  this  stage.  The  main  source  of  noise  is  the 
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fluctuations  in  coil,  though;  Nyquist's  formula  gives  this 
noise  voltage  as  u^  =  4kTRAf,  in  which  k  is  Boltzmann's  constant. 
T  is  temperature  (*K),  Af  is  bandwidth,  and  R  ir.  resistance. 


Fig.  3. 

Clearly,  one  way  of  minimizing  the  noise  is  to  make  the  ampli¬ 
fier  of  small  Af,  but  too  small  a  value  is  undesirable,  for  the 
permissible  range  of  variation  in  Z  is  then  small.  The  Af 
chosen  was c/ s  for  the  above  central  frequency  of  2075 
c/s;  thenyu  is  about  1  mv.  The  contribution  from  the  6ZhlP 
make?  the  total  only  1.5  pv  at  most,  which  means  a  signal-to- 
noise  ratio  of  7-8. 

The  high  gain*makes  a  severe  demand  on  the  screening, 
for  any  interference  entering  in  the  early  stages  may  distort 
or  even  mask  the  signal.  All  the  input  circuits  were  care¬ 
fully  screened,  each  stage  being  contained  in  a  separate  screen¬ 
ed  box,  the  tubes  and  variable  capacitors  being  enclosed  in 
special  screens.  Screened  cables  and  couplers  were  used 

*anplific  atiqfe  coefficient 

of  5  •  105  ' 
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throughout!  all  screens  being  grounded.  The  poi  -.rizing  coil 
was  short-circuited  and  grounded  when  not  carrying  current;  it 
then  acted  as  a  screen  for  the  pickup  coil. 

The  frequency-converting  unit  include:*  the  doubler,  the 
amplifying  filter,  the  mixer,  and  a  cathode  follower.  The 
doubler  reduces  ttjt,  for  a  given  Afc  (error  of  frequency  mea¬ 
surement);  it  consists  of  a  oZhlP  whose  plate  circuit  contains 
the  primary  of  a  transformer,  whose  second  is  connected  to  a 
full-wave  rectifier  (DGTs-23  diodes).  The  signal  at  the  out¬ 
put  of  the  doubler  is  reduced  in  amplitude  and  somewhat  distor¬ 
ted,  so  2fc  is  isolated  and  amplified  by  a  stage  containing  a 
6N9S.  The  mixer  compares  2fe  with  the  standard  frequency  (from 
a  quartz  crystal)  of  41 66  c/s;  at  the  plate  of  the  mixer  there 
appears  the  beat  (difference)  frequency  ffe.  A  detector  con¬ 
taining  a  DGTs-27  diode  is  used  to  isolate  f^.  The  cathode 
follower  feeds  the  loop  oscillograph  with  the  beat  frequency, 
the  bias  for  the  tube  being  produced  by  a  resistance  of  20041, 
only  the  alternating  component  reaching  the  MPO-2  oscillograph 
(via  a  20  pF  capacitor).  The  quartz  oscillator  is  designed 
to  maintain  a  fixed  frequency  within  0.00015*,  the  fundamental 
standard  being  a  100  kc  crystal,  whose  frequency  remains  con¬ 
stant  to  this  order  after  it  has  been  running  for  two  hours. 

The  frequency  of  4166.6  c/s  is  generated  by  successive  division 
by  4  and  6  from  100  kc,  the  filtered  and  amplified  product 
being  fed  to  the  mixer. 

The  100  kc  quartz  crystal  is  enclosed  in  a  thermostatic 
vessel;  it  forms  the  first  stage  in  a  self-excited  oscillator, 
where  it  replaces  the  usual  tuned  circuit  in  the  grid  circuit. 
This  stage  uses  half  a  6N1P.  The  plate  circuit  is  tuned  to 
120  kc  to  improve  the  stability  and  is  also  decoupled.  The  po¬ 
sitive  half-cycles  are  passed  to  the  4:1  division  stage  via  a 
circuit  consisting  of  a  DGTs-25  diode,  a  capacitor,  and  a  resis¬ 
tor. 

3.  Results 

The  horizontal  component  H  is  balanced  out  by  setting 
the  axis  of  the  Helmholtz  coils  in  the  plane  of  the  magnetic  me¬ 
ridian  and  in  the  horizontal  plane.  Any  such  setting  is  made 
with  a  certain  error,  so  we  must  consider  the  effect  on  the  mea¬ 
sured  Z  of  any  error  in  balancing  H. 

Let  the  coil  moke  an  angle  a  with  the  horizontal  plane 
(Fig.  4)  and  an  angle  180  -  3  with  the  plane  of  the  magnetic  me¬ 
ridian;  let  the  actual  field  of  the  coil  be  H^,  which  differs 
from  H  by  AH  =  H  -  .  Then  the  residual  field  at  the  center 

is  given  by 

Z'*  =s (Z  -f  Hx  sin *)*  +  (//  —  //,  cos  «• cos  p)*  +  //’sin*?, 
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v/hich,  if  we  neglect  quantities  of  tho  fourth  nr>."  in  4H,  u, 
and  (?  (which  are  all  small),  we  have 


The  largest  term  is  the  first,  but  this  is  readily  eliminated 
by  turning  the  coil  through  180*  and  taking  a  second  reading 
Z";  a  and  (3  then  change  sign,  so  H,sin  a  does  so  as  well. 
Then  Z'  -  Z  and  Z"  -  Z  give  us  the  relative  error  £Z/Z  as 


AZ 

Z 


1  (!«)•,  1  H*  t  .  ,  «v 
7  z*“  +  77T (•  +P  )• 


For  example,  if  the  error  in  balancing  H  is  JH/H  =  0.01 


Fig.  4. 

and  if  a  =  3  =  1*,  we  have  43/3  -  0.00055 i  of  v/hich  1/7  arises 
from  4 H  and  6/7  from  a  and  p. 

The  H/Z  found  in  the  USS.d  range  from.  0.001  (in  the  polar 
regions)  to  0,75>  so  the  4Z/Z  for  the  above  conditions  range 
from  0  to  0.0C018.  The  error  from  all  causes  can  be  kept  be¬ 
low  0.01%  if  H  is  balanced  out  to  within  a  few  %  and  the  angles 
are  correct  within  1*. 

The  current  required  by  the  coils  at  their  present  site 
is  a,bout  55  ma;  the  exact  value  may  be  determined  by  experi¬ 
ment. 

Let  the  current  be  I  for  exact  compensation,  which  gives 
a  beat  frequency  fv.  Then  currents  1^  =  I  -  41  and  IP  =  I  + 

41  give  Z'  =  Z  +  42  and  f  =  f0  +  4f,  and  so  on;  this  enables 
one  to  find  I,  and  1^  such  that  f^  =  0,  which  is  easily  obser¬ 
ved.  Then  the  correct  current  is 

!  =  (!]_  +  I  Z)/2. 
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Our  measurements  gave  I  as  53*5  m» ,  which  agrees  well 
with  the  value  calculated  above;  the  precession  signal  was  re¬ 
corded  at  intervals  of  1  ma  from  40  to  70  ma,  and  a ilo  at  inter¬ 
vals  of  0.5  na  near  the  point  of  compensation. 


.  .  Fi?*  5. 

Measurements  of  Z  were  made  at  the  Kavgolovo  laboratory 
of  the  All-Union  Metrology  Research  Institute,  which  is  outside 
the  city,  on  15  May  I960;  recordings  were  made  of  the  beat  fre¬ 
quency  and  also  of  time  marks  separated  by  0.01  sec  (Fig.  5  is 
typical).  The  19  pairs  of  measurements  are  listed  in  the  table. 
Readings  were  taken  simultaneously  fro  ;  an  M-14  magnetometer ; 
recordings  were  made  of  the  variations  in  K  with  a  variation  in¬ 
strument.  The  Z-v-ariometer  provided  corrections  and  so  enabled 
us  to  establish  the  error  of  the  method;  the  H  variations  pro¬ 
vided  a  check  on  the  deviations  from  exact  compensation.  One 
scale  division  on  the  M-14  is  equivalent  to  10  y,  and  the  read¬ 
ings  can  be  made  to  within  0.1  division,  so  the  correction  for 
variation  can  be  made  to  1  y.  The  H-variometer  recording 
showed  no  major  changes  during  t lie  experiment,  and  the  compen¬ 
sation  remained  within  the  permitted  limits. 

The  table  shows  that  the  standard  deviation  is 

•/  =  V^r-0-016  v*. 

in  which  n  is  the  number  of  observations;  then  f^  =  14.151  - 
0.016  c/s.  The  compensation-current  tests  showed  that  f0  > 

2fc ,  for  zero  beat  was  obtained  for  1^  and  1 2  that  were  respec¬ 
tively  larger  and  smaller  than  so  f0  =  2f(Z  +  iZ) ,  in  which 
f(Z  +  £Z)  is  the  precession  frequency  for  currents  and  I £• 
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M 

0k«Kv*.ton< 

m 

AS 

-Ls 

K 

(Vk)* 

4  0* 

4 1»* 

4.m 

1 

16X0 

12X0 

14X0 

7,70 

-0.20 

-0,16 

14,23 

2 

M 

12X3 

14X2 

7,70 

HE 

-0,16 

14,16 

tn 

3 

M 

12,71 

14,41 

7,70 

-0X0 

-0,16 

14.25 

OulO 

0,0100 

4 

12X6 

14X7 

7,70 

HiIm: 

-0,16 

14,21 

0.0036 

5 

'J.V 

KWH 

7X0 

HiYuj 

MlfrU 

14.16 

6 

if 

14.18 

7,80 

Mill  t 

0,00 

14,18 

7 

y./ 

/ST 

7,40 

Mil  ft 

0,06 

14,13 

8 

■ 

[jrt 

7X0 

mYV 

MifQ] 

14,19 

0,0006 

9 

1  TX 

TCT 

■  uPt 

7,40 

M'lM 

0.01 

0,0196 

10 

'  Tpf 

■ 

7,40 

■tu! 

0X6 

14.01 

0.14 

0,0196 

11 

"irr 

'  Tft 

■ 

7.40 

^rXT; 

0X6 

14.10 

EEJ 

0,0075 

12 

16,71 

■  tT 

7X0 

■ijri  i 

14,16 

rm 

0,0001 

13 

16.71 

'  ^Tv 

■  ijr 

7X0 

moW 

0.0J 

14,16 

0,01 

0X001 

14 

IDlj 

:*r 

■  Tmv 

7,80 

0.00 

14X6 

0,10 

0.0100 

15 

16X6 

12X6 

14,12 

7X0 

0,00 

MiMJ 

14.12 

0,0009 

16 

16X7 

12X1 

14,14 

7X0 

oxo 

0.00 

14,14 

fiVol 

0.0001 

17 

Era 

12X4 

14,10 

7.50 

iM  WJl 

mnm 

0,0005 

18 

IM 

12X6 

14X3 

7X0 

0.00 

14,23 

iyu 

0,00:>l 

19 

16.77 

4 

12X3 

14X0 

7X0 

0,00 

14X0 

0.05 

0.0025 

Mta* 

7.60 

14,151 

S  0,0874 

But  f(Z  +  A2)  V/-..S  1  =r<  <?r  '  -  ’ i*  f(Z)  (-  fc>»  ■—  1  ^  >p''' v;  ’  v  ;-L':h 
means  that  ffc  ff-ust  bo  .;-a, :.r->ct<-J  r rot.,  f to  .»■•'«.■  .1  f  c ,  -e 

fc  -  1/2(4166.667  -  14,151)-  (2676.258  ±  0,008)  c/s- 


but  V  =  x  ’-C;‘  (in  ••iu:is'iBec'1  )  no 


•y  6.283 1 8 •  2076.26  /xotxsc  o  1  a  o\  ia- 

z °  "IMR-TV"  =  (48765-8  ±  °-2> * 10  *• 
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USING  NATURAL  ELECTROMAGNETIC  VARIATIONS  FOR  INVESTIGATION 
OF  THE  EARTH'S  UPPER  LAYERS 

M.  N.  Berdichevskiy,  B.  Ye.  Bryunelli,  A.  Ye.  Lantsov, 
and  O.  M.  Reepopov 

pp.  49-55 


Introduction 


The  theory  developed  by  Tikhonov  /6/»  Cagniard  and 

others  has  been  used  at  Leningrad  University  and  at  the  All- 
Union  Geophysics  Research  Institute  in  new  magnetotelluric 
survey  methods,  namely  magnetotelluric  profiling  (MTP)  and  mag¬ 
netovariation  surveying  (MVS)  /%4/*  MTP  consists  in  record¬ 
ing  variations  of  period  10-60  sec  (short-period  variations) 
in  the  geoelectric  field  and  in  the  horizontal  components  of 
the  geomagnetic  field  At  one  point  only.  The  mean  ratio  of 
these  variations  indicates  the  total  conductance  of  the  sec¬ 
tion,  which  gives  an  indication  of  the  relief  in  the  reference 
(high-resistance)  horizon.  MVS  is  analogous  to  the  widely 
used  telluric-current  methods  /27,  except  that  here  the  geo¬ 
magnetic  variations  are  recorded  instead  of  the  geoelectric 
ones.  The  maximal  variations  in  the  horizontal  components 
are  associated  with  increased  values  of  the  total  conductance, 
which  are  commonly  related  to  troughs  in  the  crystalline  base¬ 
ment  . 


These  methods  have  become  of  practical  value  in  studies 
on  the  upper  layers  of  the  crust  as  a  result  of  the  design  of 
a  sensitive  short-period  magnetometer  /5/  by  Bryunelli  at  the 
Department  of  the  Physics  of  the  Earth's  Crust  at  Leningrad 
University.  The  first  field  trials  of  magnetotelluric  methods 
were  made  in  1957  in  the  Tyumen  region  Wi  and  the  work  was 
extended  in  1958  by  the  prr  >«nt  authors  to  the  Druskininkag- 
Prenay-Sovetsk  profile  (about  250  km  long)  in  the  Polieh- 
Lithuanian  depression,  which  represents  a  major  trough  in  the 
crystalline  basement  filled  with  sediments  of  Paleozoic  age. 

The  eentral  part  contains  Permian  gypsum  beds  within  the 
|  Paleozoic  sediments |  these  beds  are  of  very  high  resistance  _ j 

"Trumans kaya  Oblas$ 
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Fig.  1.  Magneto tell uric  results  for  the  Druskininkay- 
Sovetsk  profile:  1)  mean  geomagnetic  field  Kjj  from  MVS; 

2)  total  conductance  S  of  sediments  from  MTP;  3)  recip¬ 
rocal  of  the  mean  telluric  field f  from  telluric-current 
methods;  4)  surface  of  Permian  gypsum  from  resistance 
surveys;  5)  surface  of  crystalline  rocks  from  MTP;  6) 
total  conductance  of  sediments  from  vertical  and  frequency 
sounding  techniques;  7)  surface  of  crystalline  basement 
From  seismic  surveys;  8)  deep  of  the  crystalline  base¬ 
ment  from  drilling. 

and  represent  an  obstacle  to  the  use  of  vertical  electrical 
sounding  for  examining  the  crystalline  basement.  Figure  1 
shows  a  section  along  the  above  profile,  which  shows  that 
the  thickness  of  the  sediment  cover  rises  from  400  m  near 

l_  J 
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Druskin inkny  to  P100  ni  near  I 
Uovetak  (from  drillings). 

The  apparatus  for 
the  magne to telluric  surveys 
consisted  of  a  fo"r-channel 
oscillograph  coupled  to  a 
switching  box,  two  grounded 
electrical  measuring  lines, 
two  H-magnetometers ,  and  a 
radio  station  with  remote- 
controlled  switch.  The 
lines  were  at  right  -ingles 
and  had  an  arbitrary  azi¬ 
muth,  while  the  magneto¬ 
meters  were  set  parallel  to 
the  lines  in  shallow  pits. 

The  observations  were*  made 
simultaneouciy  *-t.  two  points, 
magnetometers;  4)  radio  station;  as  i  n  telJ  nric-currer.t 
5)  remote-controlled  switch.  methods,  the  operations 

being  synchronized  with 
remote-controlled  switches  (Fig.  2). 

1.  Fro-ie.  -in g  and  h'eruJ.ts  from  KTP  i-tecordlngt  . 

The  processing  consists  of  measuring  the  apparent  ampli¬ 
tudes  and  periods  of  1 5 -TO  qunsisinusoi  dal  pulses  ir.  Ex ,  Ey ,  i!*, 
sr.'t  riy ,  from  which  are  deduced  the  principal  axes  of  the  polari¬ 
zation  ellipses  for  the  telluric  and  magnetic  rields.  The 
mesn  periled  of  the  aseoci;- 1>- J  variations  is  given  by 


Fig.  2.  -.poor.itus  for  magneto- 
telluric  studies:  1)  oscillo¬ 
graph;  2)  .-•-•itching  b<  x;  3) 


TXEX  -f  T y£y 

Ex  •+•  Ey 


(D 


in  which  T  j. :•  the  mean  period  for  r,  (in  .. ec)  and  T  is  the 
.  x  j 

for  R  . 

J 

The  resulting  K,  it,  and  T  are  inserted  in  the  general 
MT?  formula  for  horizontally  homogeneous  stretificot ion  /h/’ 


i'=796[d+/_rj,  «> 

in  which  S  is  the  total  longitudinal  conductance  (mho)  and  pn 
is  the  specific  resistance  of  the  reference  horizon  (A-m)  . 

A  statistical  treatment  of  the  individual  results  eliminates 
the  dependence  or.  azimuth  and  so  improves  the  accuracy  of  S . 

The  pn  of  {cl)  was  taken  su  700  A -m  in  the  calculations;  this 
was  the  value  given  by  our  magne totelluric  soundings  on  a 
nearly  area.  Figure  3  shows  a  magnetotell uric  sounding  curve,  j 
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the  ordinate  being  the 
apparent  resistance  pt 
A-m  as  given  by  the 
formula 


ir. 


Pt~0,2T 


m- 


in  which  Ex  is  the  ampli¬ 
tude  of  the  telluric 
variations  (pv/km) ,  rly 
is  the  amplitude  of  the 
geomagnetic  variations 
(y),  and  T  is  the  period 
(sec) , 

The  values  of  E 


Fig.  3- 
curve . 


Kagnetotelluric  sounding 


x» 

Hy,  and  T  were  deduced  by 
harmonic  analysis  of  the 
pulses,  the  latitudinal 
component  of  II  and  the 
longitudinal  component  of 
E  being  used.  The  result¬ 
ing  curve  approaches 
closely  the  right-hand  branch  of  the  theoretical  curve  for  a 
two-iayer  section  lying  above  a  high-resistance  base.  A 
method  analogous  to  that  already  described  £Z7  was  used  to 
interpret  this  curve  and  to  determine  the  reference  horizon. 

Figure  1  shows  the  curve  for  the  S  of  the  sediments, 
which  rises  from  ?7  mho  near  Druskininkay  to  800  mho  near 
Sovetsk.  Repeat  measurements  indicate  that  the  error  in  3 
ranges  from  1  to  7%,  with  a  i^ean  of  3>«,  with  values  of  15-30% 
at  a  few  points  \ here  the  hign-resistance  gypsum  deposits  tail 
off  and  the  wave  pattern  becomes  very  complicated.  The  MTP 
results  agree  well  with  those  from  logo  in  deep  boreholes  and 
from  previous  electrical  surveys,  which  shows  that  the  crystal¬ 
line  basement  acts  as  the  reference  horizon  in  this  region  in 
spite  of  the  high-resistance  gypsum  beds. 

The  MTP  results  were  given  a  quantitative  interpretation 
on  the -basis  of  drillings  at  oovetsk,  Prenay,  and  Druskininkay, 
which  extended  down  to  the  crystalline  rocks.  The  mean  0e 
(longitudinal  resistance)  of  the  sediments  between  boreholes 
was  deduced  by  linear  interpolation;  the  thickness  of  the 
sediments  was  deduced  from 

h  =  pe-S. 

Figure  1  shows  the  results  in  the  form  of  the  relief  of 
the  basement.  The  results  agree  with  those  from  seismic 
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surveys . 


2.  Processing  and  Results  from  MVS  Recordings 

The  ellipses  aethod  used  for  telluric  currents  is 
employed  here;  the  points  derived  from  transforming  the 
magnetic  variation  vectors  lie  closely  around  arcs  of  ellipses, 
which  shows  that  the  fluctuating  horizontal  components  of  the 
geomagnetic  field  at  points  p  and  q  are  related  by 

es  aHXp  -f  bHup, 

H„  =  eHXp  +  df/„t 

in  which  a,  b,  c,  and  d  are  bonstants.  Theory  shows  that 
these  are  dependent  on  the  directions  of  the  x  and  y  axes  as 
well  as  on  the  electrical  parameters  of  the  rocks  around  p 
and  q.  The  principal  parameter  used  in  MVS  is  Kg,  which  we 
call  the  mean  magnetic  strength;  it  is  the  square  root  of  the 
ratio  of  the  synchronous  closed  hodographs  for  the  magnetic 
field  and  is  given  by 


in  which  A  and  B  are  the  semiaxes  of  the  ellipse  and  R  is  the 
radius  of  the  circle  whose  transform  is  that  ellipse.  The 
efror  in  the  determination  of  this  parameter  is  found  to  be 
2-4%  in  our  numerous  repeat  measurements. 

Figure  1  shows  Kg  for  the  Druskininkay -Sovetsk,  profile; 
it  also  shows  Kg,  the  mean  telluric  field.  Both  are  taken  as 
unity  for  the  region  around  Druskininkay.  Kg  increases  as 
the  basement  becomes  deeper  (as  the  theory  of  MVS  would  indi¬ 
cate),  becoming  1.42  at  Sovetsk,  the  total  anomaly  thus  being 
42%,  which  is  at  least  10  times  the  possible  experimental 
error.  Some  details  of  the  structure  can  be  seen  on  the 
curve  for  Kg,  in  particular  the  margins  and  the  central 
(deepest)  part  of  the  depression.  The  relative  anomaly  in 
the  alternating  geomagnetic  field  is  much  less  than  that  in 
the  telluric-current  field  (the  42%  change  in  Kg  corresponds 
to  a  sixfold  change  in  Kg),  so  MVS  is  likely  to  be  of  value 
in  regional  studies  for  areas  in  which  telluric-current 
methods  cannot  be  used.  These  results  confirm  that  the 
main  theoretical  deductions  concerning  these  methods  are 
correct. 

In  conclusion,  we  oonsider  the  future  prospects  for 
the  above  two  methods.  The  MTP  method  is  of  the  greater 
practical  interest;  it  seems  likely  to  replace  telluric- 
current  methods  in  many  regions  and  may  well  greatly  reduce 
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the  ooet  of  electricil  surveys,  for  it  eliminates  the  need  for 
reference  grids  and  base  measurements.  Moreover,  RTP  gives'  a 
more  precise  indication  of  the  structure,  for  the  total  conduct¬ 
ance  is  determined  directly.  In  particular,  it  substantially 
simplifies  the  stratigraphic  assignment  of  the  reference  horizon, 
and  it  may  improve  the  quantitative  calculations  so  necessary 
in  geology. 

The  MVS  method  may  prove  of  great  value  in  regional 
studies  on  the  permafrost  regions  of  Siberia,  for  insulating 
horizons  have  very  little  effect  on  the  alternating  geomagnetic 
field.  Moreover,  MVS  may  be  more  productive  than  telluric- 
current  methode,  for  it  can  be  operated  with  simple  automatic 
equipment. 
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PLOTTING  MAGNETOTELLUrilC  SOUNDING  CURVES  ACCORDING  TO 
RECORDINGS  OF  SHORT-PERIOD  VARIATIONS  IN  THE  .i*  JURAL 
GEOMAGNETIC  FIELD 

pp.  56-66  A.  V.  Kovtan 

1.  The  following  assumptions  form  the  basis  of  the 
theory  of  magnetotelluric  sounding  (NTS).  Fluctuations  in  the 
electric  currents  in  the  upper  atmosphere  produce  plane  electro¬ 
magnetic  waves  whose  time  dependence  is  exp(iwt).  The  Earth 
is  assumed  to  be  flat  and  to  consist  of  homogeneous  layers,  so 
Maxwell's  equations  can  be  solved  (apart  fron  a  constant  factor 
representing  the  amplitude  of  the  incident  wave),  the  impedance 
at  the  surface  of  the  Earth  being  then 


in  which  Ex  and  H  are  the  mutually  perpendicular  horizontal 
field  components  at  the  surface  A ,6,2^ ,  these  being  dependent 
solely  on  frequency  and  on  the  structure  of  the  Earth.  The 
individual  waves  E^w)  and  Hy(w)  are  isolated  fro.i  simultaneous 
recordings  of  Ex  and  Hy  to  find  the  relation  of  Z  to  •*>,  which 
is  compared  with  theoretical  curves  to  determine  the  nature  and 
parameters  of  the  structure. 

These  natural  variations  of  periods  from  5  sec  upwards 
have  ‘been  used  to  examine  the  electrical  structure  at  great 
depths  since  Bryunelli's  high-sensitivity  short-period  magneto¬ 
meter  was  developed  at  the  Department  of  the  Physics  of  the 
Earth’s  crust  at  Leningrad  University,  thi3  be  in.”  suit  able  for 
use  in  the  field  /l, 2, £7.  The  following  deductions  are  based 
on  results  from  two  expeditions  made  jointly  with  the  All-Union 
Geophysics  Research  Institute  in  conjunction  with  the  Borok 
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geophysics  station*. 

Leaving  aside  any  discussion  of  the  theoretical  assump¬ 
tions  and  of  the  technical  methods  of  recording  the  variations 
(for  which  see  /1,2,8/),  I  now  discuss  the  practical  aspects  of 
the  isolation  of  particular  frequencies  from,  the  recordings,  the 
object  being  to  construct  MTS  curves. 

2,  The  short-period  variations  may  be  divided  into  per¬ 
sistent  (Pc)  and  transient  (Pt)  types;  the  first  type  may  per¬ 
sist  for  hours  with  only  slight  changes  in  period  and  amplitude, 
the  usual  ranges  being  respectively  10-120  sec  and  1-2  y.  The 
Pc  are  most  common  in  the  morning  hours;  they  are  not  always 
strong,  and  as  a  rule  they  have  only  one  main  period  on  any 
given  day.  The  Pt  take  the  form  of  short  bursts  of  oscilla¬ 
tions  of  decreasing  amplitude;  they  are  best  seen  during  the 
night  hours.  Even  the  Pc  are  not  steady-state  oscillations, 
though,  which  makes  it  difficult  to  isolate  particular  frequen¬ 
cies. 


Harmonic  analysis  of  the  recordings  is  essential  if  the 
recordings  do  not  indicate  a  steady  state.  Let  the  incident 
fields  E°  and  H°  represent  a  plane  pulse  of  arbitrary  but  finite 
duration;  the  Yields  recorded  at  the  surface  are  Ex(t)  and 
Hy(t),  which  can  be  represented  as  Fourier  integrals: 

*v 

E*  W*8  i  j  *** E*  w  dm> 

V 

—  SC 

in  which  •*  denotes  quantities  representing  spectral  densities 
at  the  Earth's  surface: 


x- 

?*(•)  =  J  e~MEx{t)dt, 

m 

«,(•)=  | 

—  m 


‘About  10  years  have  passed  since  Tikhonov's  and 
Cagniard's  papers  appeared,  but  f-. To  curves  have  not  yet  come 
into  use,  though  the  simpler  and  more  rapid  magne totelluric 
profiling  (hTP)  /! ,2/  h  as  been  used. 
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The  ratio  of  these  spectral  densities  is  the  impedance  of  the 
plane  wave  at  the  surface: 

•  (5) 

/*,<•) 

That  is,  any  varying  oscillation  should,  strictly  speak¬ 
ing,  be  handlfd  by  reference  to  the  Fourier  integrals  for  Ex(t) 
and  H  (t)  if  Z  is  to  be  found.  Numerical  calculations 
show  that  the  build-up  time  for  the  electric  field  is  virtually 


Fig.  IV  BuiTd-up  of  an  elec¬ 
tromagnetic,  field  over  a  homo¬ 
geneous  half-space:  1)  steady- 
state  field;  2)  shape  of 
field  pulse  II  if  the  incident 
wave  terminates  at  point  1; 

0  =  c/t  (ratio  of  conductivity 
to  dielec  trip  constant  for  a 
homogeneous  Earth. 


independent  of  the  conductivity;  this  time  is  proportional  to 
T  *  2%/u>  if  the  Earth  may  be  treated  as  a  uniform  half-space, 
and  the  steady  value  is  reached  within  10$»  after  a  time  T/3 
from  the  start  of  the  plane  wave.  Transient  effects  for  the 
magnetic  field  over  a  conducting  medium  are  negligible,  for  the 
correction  for  the  transient  component  is  small  relative  to  the 
steady-state  component.  Figure  1  illustrates  the  build-up  of 
an  electromagnetic  field  above  a  homogeneous  body. 

3.  A  correct  conception  of  the  buildup  enables  one  to 
avoid  tedious  calculations  on  Fourier  integrals  in  many  cases; 
the  'visible  amplitudes'  may  be  used  to  isolate  frequencies. 

This  method  is  as  follows.  The  recordings  are  examined  for  a 
train  of  pulses  whose  shape  is  close  to  sinusoidal,  the  dis¬ 
tance  between  adjacent  maxima  being  taken  as  the  period;  the 
visible  amplitude  of  the  train  is  then  twice  the  amplitude  for 
that  ppriod  (Fig.  1).  It  is  necessary  to  be  very  careful  in 
the  use  of  this  method,  as  calculations  on  transient  effects 
show,  for  the  following  reasons. 

1)  The  method  is  not  applicable  to  single  isolated  field 
deviations. 

2)  The  method  is  applicable  only  to  the  middle  part  of 
a  train  of  two  or  three  waves  differing  only  slightly  in  ampli¬ 
tude  and  period;  the  first  and  last  waves  *:.ust  be  excluded. 
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3)  Transient  effects  usually  alter*  the  shape  of  the 
electric-field  wave;  any  lack  of  oorrespondance  between  the 
wave  shapes  for  the  two  fields  (provided  that  the  two  lines  have 
identical  characteristics) ,  and  especially  any  difference  in  the 
distances  between  maxima,  is  a  direct  indication  of  transient 
effects,  so  such  waves  must  be  excluded.  The  Pc  type  is  the 
best  for  this  purpose,  for  the  waves  can  usually  be  treated  as 
sinusoidal,  in  spite  of  some  variation  in  amplitude  and  period. 
It  is  seldom  possible  to  use  the  Pt  type  for  this  purpose. 

4,  The  above  method  was  used  with  simultaneous  recor¬ 
dings  made  around  the  Baltic  coast  in  the  summer  of  1958  at 
points  along  the  line  from  Prenay  to  Sovetsk.  The  recordings 
were  of  inadequate  duration  and  related  mainly  to  one  time  of 
day,  so  the  spectrum  of  the  oscillations  was  rather  restricted; 
only  short  sections  of  the  frequency  curves  could  be  construe-’ 
ted.  More  interesting  results  were  obtained  from  continuous 


Fig.  2.  MTS  curves! 
a)  NS,  b)  E'tf  as  con¬ 
structed  1)  from  vis¬ 
ible  amplitudes,  2) 
by  harmonic  analysis, 

3)  from  films  with  T  = 
11  sec  after  correction 
for  sensitivity. 


recordings  at  Borok  from 
January  to  October  1959 
Figure  2  shows  curves  con¬ 
structed  a)  for  the  north 
electric  component  and  east 
magnetic  component,  and  b) 
for  the  east  and  north  com- 
ponents  as  above.  H ere^/T and 

P  infl-m,  being  given  by 


(4) 


in  which  T  is  in  sec,  Ex  and  H  being  measured  respectively  in 
mv/km  and  y.  The  recordings  Were  continuous  and  were  made 
with  a  strip  speed  of  6  mm/min,  which  provided  a  means  of  exa¬ 
mining  T  between  10  and  700  sec.  The  lack  of  agreement  bet¬ 
ween  curves  a  and  b  is  a  result  of  inexact  setting  of  the  mag¬ 
netometers  relative  to  the  directions  used  for  the  electric 


components.  Recalculation  (correction)  was  possible  for  T  >  200 
sec,  and  this  gave  the  same  (  for  the  two  directions.  This  type 
of  error  merely  causes  a  parallel  displacement  of  the  curve, 
which  does  not  interfere  with  a  study  of  the  structure.  PhyGi- 
cal  meanings  can  be  given  to  the  average  parameters  determined 
independently  from  the  two  curves. 

5.  As  a  rule,  the  f  given  by  different  recordings  for 
a  given  period,  and  even  sometimes  by  one  recording,  are  not  the 
same;  the  spread  in  the  f  derived  by  means  of  visible  amplitudes 
may  be  as  much  as  40-60%,  and  it  is  not  a  consequence  of  inade¬ 
quacy  in  the  theory.  The  spread  is  largely  random  and  results 
particularly  from  errors  in  the  electric-field  measurements. 

Calculations  show  that  the  maximum  error  in  the  g  of  (4) 
may  be  as  much  as  80%  in  some  cases,  the  sources  of  this  being 
errors  in  the  measurements  of  E*  and  Hy  (these  may  be  10-20% 
for  a  mean  amplitude  of  10  mm),  errors^in  the  calibration  (not 
less  than  5%),  and  errors  resulting  from  the  frequency  response 
of  the  instruments  (at  least  5%).  This  means  that  no  physical 
significance  can  be  assigned  to  individual  values;  only  the 
general  trend  (from  statistical  treatment)  is  significant.  For 
example,  the  EW  curve  of  Fig.  2  represents  the  result  of  recor¬ 
dings  of  total  duration  over  100  hr.  The  mean-square  error  of 
£in  some  cases  exceeded  15%,  iri  part  because  few  readings  were 
used.  A  major  improvement  hero  can  come  only  from  the  use  of 
more  sensitive  equipment. 

There  are  also  systematic  errors  in  the  results;  for 
example,  the  frequency  characteristics  of  the  lines  were  known 
only  for  T  £  400  sec.  Adjustments  were  made  for  T  >  400  sec 
by  means  of  theoretical  formulas  for  the  magnetometers,  which 
may  have  given  some  systematic  errors  in  the  g.  Another  source 
of  systematic  error  occurs  in  the  calibration,  for  which  square 
waves  are  used.  The  differing  frequency  characteristics  of  the 
recording  instruments  for  E  and  H  cause  an  appreciable  distor¬ 
tion  of  the  calibration  pulses  (especially  for  the  magnetic  mea¬ 
surements).  The  amplitude  of  the  pulse  at  the  start  may  differ 
from  that  at  the  end  by  1-2  mm.  It  is  usual  to  take  mean 
values  and  to  assign  the  resulting  sensitivity  to  the  frequency 
for  which  the  response  of  the  instrument  is  maximal;  but  this 
is  correct  only  if  the  maximum  in  the  spectral  density  for  the 
pulse  porr^ponds  to  the  frequency  of  maximum  response,  which 
is  not  always  so.  Some  of  the  recordings  at  Borok  were  made 
with  the  frequency  response  of  Fig.  3a  (t  •  11  sec)  for  the  mag¬ 
netic  channel,  for  example.  Here  pulses  of  duration  1  rain  were 
used  in  the  calibration;  Fig.  3b  shows  the  spectral  density  of 
the  pulse,  whose  maximum  lies  at  180-200  sec  for  a  duration  of 
60  sec;  the  frequency  of  maximum  response  corresponds  to  30  sec, 
though.  This  feature  was  neglected  in  the  initial  processing, 
so  all  f  from  this  recording  deviated  greatly  from  the  rnean,£. 
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Pig.  3«  Change  in  pulse  spectrum  after  passage  through 
a  magnetic  measuring  line:  a)  frequency  response  of 
the  line  (T  =  11  sec),  b)  spectral  density  of  a  square 
pulse;  1)  before,  2)  after  passage  through  the  line. 

Correction  for  ?  =  11  sec  brought  all  the  (  very  close  to  j*  (Fig. 
2).  This  example  shows  that  care  is  needed  in  order  to  elimi¬ 
nate  errors  of  this  type. 

6.  The  Pc  type  rather  rarely  has  T  *  100  sec,  so  it  is 
necessary  to  consider  bursts  of  oscillations  in  order  to  exa¬ 
mine  the  larger  T.  The  main  difficulty  here  is  to  select  a 
section  for  harmonic  analysis,  for  isolated  oscillations  of  fair 
amplitude  are  not  very  common.  Incorrect  termination  of  a 
burst  can  give  rise  to  large  errors  in  the  spectral  density. 

The  cause  of  the  errors  is  as  follows.  Let  U3  suppose 
that  the  incident  wave  is  an  arbitrary  function  of  time;  the 
spectral  density  of  the  incident  field  has  a  particular  leaning 
for  each  instant,  and  so  the  signal  recorded  at  the  Earth's  sur¬ 
face  from  a  fraction  of  the  incident  wave  has  a  different  but 
related  spectral  density.  However,  the  main  difficulty  is  that 
the  corresponding  electric  and  magnetic  signals  are  not  readily 
isolated. 

’  The  following  example  illustrates  this.  Let  the  incident 
wave  be  sinusoidal  (Fig.  1)  and  stop  at  the  time  corresponding 
to  point  1.  The  Ex  recorded  at  the  surface  will  persist  past 
this  point  (on  account  of  the  transient  response),  but  the  corre¬ 
sponding  magnetic  field  will  vanish  almost  instantly,  for  the 
transient  response  for  this  is  negligible.  It  is  clear  that, 
if  (3)  applies  precisely  to  the  Hy  and  Ex  of  pulse  II,  the  spec¬ 
tral  densities  of  pulses  (0-1)  in  Ex  and  Hy  Jo  not  satisfy  this 
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relation;  the  ratio  of  the  spectral  densities  of  the  two  has  no 
physical  significance.  The  longer  the  pulse,  the  smaller  the 
error  introduced  by  transient  effects  and  the  closer  the  spect¬ 
ral  densities  to  the  true  values;  but  it  i?  difficult  to  esti¬ 
mate  the  error  without  calculations  on  the  build-up  times  for 
electromagnetic  processes  in  the  particular  medium.  This  means 
that  a  pulse  should  not  be  terminated  artificially  unless  there 
is  some  especially  compelling  reason. 

Gol'tsman's  method  C^y^J  is  convenient  for  the  analysis 
of  single  pulses,  for  it  enables  to  replace  the  integration  of 
(2)  by  a  summation  without  incurring  major  errors.  The  corre¬ 
sponding  theorem  is  as  follows  fl) : 

The  spectral  density  S(u>)  of  any  function  f(t)  having  a 
spectrum  restricted  to  the  region  |w»|  <.  may  be  put  as 

_  • 

£(•)»&*  (5) 

kT-m 

in  which  f(kdt)  are  values  of  f(t)  taken  at  intervals  /It  =  n/u^. 

The  limiting  frequency  w0  is  readily  determined  from  the 
shape  of  the  pulse  in  the  following  way:  At  is  made  such  that 
the  narrowest  peak  has  4-6  points  on  it.  The  summation  of  (5) 
presents  no  difficulty  if  the  number  of  points  is  small;  vector 
addition  of  the  terms  is  sufficient  if  there  are  only  20-30 
terms.  Moreover,  for  numbers  of  points  less  than  73  one  can 
find  the  spectral  density  at  equally  spaced  points  by  means  of 
nomograms  or  a  harmonic  analyzer  lk7 . 

Several  pulse  trains  differing  in  length  and  spectral 
composition  were  chosen;  the  spectral  densities  were  such  as 
to  provide  values  of  f  for  T  from  36  to  4000  sec.  Figure  4 
shows  one  of  these.  The  error  in  £  may  be  40-50%,  on  account 
of  large  errors  in  the  amplitudes,  so  the  maximum  error  was  con¬ 
sidered  for  each  case;  only  those  frequencies  for  which  £  was 
likely  to  be  in  error  by  less  than  1 5>*  were  used.  The  resul¬ 
ting  points  (Fig.  2)  lie  close  to  the  curve  given  by  the  visible 
amplitude  method. 

7.  The  E V/  curve  (Fig.  2)  was  constructed  for  the  range 
of  T  from  10  to  4000  sec  in  these  two  .vays,  the  range  for  the  IIS 
curve  being  10  to  1000  sec.  These  curves  do  not  provide  a  reli 
able  interpretation  of  the  structure  around  3orok  without  refer¬ 
ence  to  additional  information  on  the  upper  layers  (ones  extend¬ 
ing  down  to  depths  of  about  3  km).  The  conductance  of  the 
upper  layers  near  Borok  is  found  reliably  as  1600  mho;  the  ge¬ 
neral  structure  at  depth  can  be  deduced.  The  slope  of  the  low- 
frequency  asymptote  indicates  that  the  ratio  of  the  resistivity 
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Q f  fcbf  crystalline  basement  to  the  resistivity  of  the  conducting 

lay arc  ia  ovar  300, 

The  bsbavipr  of  f  for  T  email  indicates  that  the  conduc¬ 
ting  layers  have  a  complicated  structure.  The  fell  in  f  for  T 
l*rge  (which  needs  to  be  checked)  would  correspond  to  a  sharp 
increase  in  conductivity  at  a  depth  of  300-300  km. 


Fig.  4*  Spectra  of  field  pulses i  a)  train  of  dura¬ 
tion  37  min,  amplitudes  read  at  intervals  of  1  min; 
b)  spectral  densities  of  ?  and  tT  (the  broken  line  is 
the  curve  for  ff  corrected  for  the  frequency  response 
of  the  instrument). 

8.  The  above  results  indicate  that  the  following  condi¬ 
tions  must  be  satisfied  in  order  to  obtain  reliable  information: 

1)  The  amplitude-  and  phase-response  curves  of  the  in¬ 
struments  rscording  E  and  H  must  be  carefully  examined; 

2)  Calibrations  should  bo  made  with  pulse  trains;  if 
square  waves  are  used,  the  conversion  factors  for  instruments 
differing  in'  frequency  response  must  be  examined  carefully; 

3)  The  visible-amplitude  method  may  be  used,  provided 
allowance  is  made  for  all  distortions  that  may  be  introduced  by 
transient  effects;  the  results  must  be  given  a  statistical 
treatment  in  order  to  minimize  random  errors; 

4)  Waves  of  T  100  sec  are  best  examined  by  analysis 
of  short  pulse  trains  of  high  amplitude,  Fourier  integrals 
being  used  in  conjunction  with  the  above  theorem  in  order  to 
simplify  the  calculations, 

I  am  indebted  to  students  R.  Latypova,  0.  Pokrovskiy,  and 

Rudnava  for  assistance  with  the  calculations. 
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CALCULATION  07  THE  LOW- FREQUENCY  DIPOLE  FIELD  IN  THE  FAN  ZONE 
L.  B.  auinuko  and  0.  P.  Sholpo 

pp.  67-77 


Introduction 

W#  have  previously  shown  that  functions  ?  and  7*  en¬ 
able  one  to  express  the  field  components  in. the  far  tone  for  a 
dipole  at  the  surfaoe  of  a  horisontally  layered  Barth: 


(0,1) 

at 

7*  (?)«=  —  -Jr  J  A  (mr)  "X+Z  • 

(0,2) 

la  whloh  r  is  the  distanoe  from  the  source  to  the  point  of 
observation,  a  is  a  coaplex  parameter  dependent  on  frequency 
and  on  the  electrical  structure  of  the  section  q  *  ar,  and  J0 
and  Ji  are  Bessel  functions  of  the  first  kind*  Wf  have  also 
shown  that  the  earlier  expressions  [7J  for  7  and  f  for  real 
argunents  nay  be  extended  to  coaplex  arguments. 

Here  we  deduoe  the  aaln  properties  of  f  and  F',  as  well 
as  those  of  the  related  functions  Q0  and  0  we  give  sons 
calculations  on  frequeney  sounding  curves  For  the  far  sons 
above  a  three-layer  section.  It  is  convenient  in  what  follows 
to  use  these  expressions: 


^<♦>  =  -£  +  0. ,<*). 


(0,3) 

(0,4) 
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Also,  SQ  sad  0.1  arc  defined  by 


(0,5) 


(0,6) 


For  simplicity,  ws  uss  ths  notation  b  •  is,  so  b  *  |b|exp(i^), 
in  which  }v|  |  <  *T.  Standard  integrals  /3»V  give  us  that 

•  (°.7) 

(»*»(»■*> — 

•  (A 


Soas  properties  of  0Q  and  3_1  are  considered  below. 

1 .  Integral  flcpreeentatione 

Expressions  (0,7)  end  (0,8)  in  principle  enable  one  to 
calculate  the  d  functions  froa  Struve's  function  and  Bessel 
functions  of  the  second  kind|  such  calculations  are  convenient 
if  )rb)  <  7,  while  asysstotie  expressions  are  of  practically 
useful  accuracy  for  |rb|  >  14.  This  sakes  it  convenient  to 
use  Quadratures  of  the  integral  representations  for  the  range 
7  A  |rb|  <  14.  The  integral  representation  of  Struve's 
function  737 


t  *r'  ,_i 

Mw-r.w+r^4jr(|j  J  +  V 

—  •y-fP<arg«<  all’)) 


give  us  that 


0,w“  J  ‘"7777i! 


ii,D 


J 
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th*  condition*  for  0  and  a rg  a  being  a*  boforo.  It  la  clearly 
■or*  convenient  to  p*rfor*  integration  with  respect  to  th*  real 
axi*|  if  art  q  <  0,  th*  l*t*cral  r*pr*a*atation*  tak*  th*  foni 


=  (ifg?<0); 

(1,5) 

• 

0.,(W— ■y,  1  (irgf  <0). 

(1,4a) 

If  arg  q  >  0,  integration  along  th*  r*al  axis  give*  an 
additional  integral  along  a  section  froa  u  a  q  to  e**xp(if)) , 

0  >  arg  q.  This  latter  integral  i*  included  to  giv* 

V 

=  f  .  +  2#*(# );  (iff  1  >  0) 

(1,3b) 

y  r  ' 

**&alagou*ly 

O-i(if) —  ,+«*.<»);  (ugy>0). 

(1,4b) 

Practical  application*  aoa*tia*s  sake  it  conv*niant  to 
ua*  int*fral  r*pr*a*ntationa  of  th*  3  in  t*raa  of  B*aa*l  func¬ 
tion*  of  th*  third  kind;  th*a*  ar*  obtained  by  replacing  th* 
B*a*«l  function*  of  th*  flrat  kind  in  (0,5)  and  (0,6)  by  coa- 
bination*  of  thoa*  of  th*  third  kinds 


Then  th*  integral  with  r**p*ot  to  th*  real  axis  split*  into  twp: 

0,(rb)  =  ^  J  M”  (rx)  -f  4"  J  ^^rx)  x  fd  ’  (1,5) 


0_,<i*)  -  -  *V'(rx) 


rx) 


*4* 


*  +  >’  (1,6) 


Mow  w*  transfer  to  integration  in  th*  plan*  of  th*  coaplex 


i  variable  z  *  x  +  iy;  Hq(z)  and  n\(z)  tjpd  to  zero  as  z  ■♦•in 
the  upper  half-plane,  as  do  l$(z)  and.  «y(z)  in  the  lower  half¬ 
plane,  so  the  integrals  containing  ff^(x)  and  H*f(x)  say  be  re¬ 
placed  by  ones  with  respect  to  the  positive  imaginary  axis  (0 
to  i«e)  and  those  containing  h%\x)  and  Imx)  by  ones  with 
reapeot  to  the  negative  laaginary  axis  (0  to  -ins). 

If  (arg  bj  *  Tt/2  (arg  q  «e  0),  the  poles  of  the  integrands 
in  (1,5)  and  (1,6)  lie  in  the  left  half-plane,  so  by  alteration 
of  the  integration  paths  we  have 


0, 

0-i(W=  —  -r|^i(y)~ji^i 


(«rf?<0), 

(irg?<0). 


(1,7a) 


(1,8a) 


If  now  arg  b  >  "^/Z  (arg  q  >0),  we  have  to  alter  the 
paths  for  the  integrals  containing  h|  and  to  avoid  a  pole 
in  the  fourth  quadrant.  The  integrals  of  the  K  functions  are 
thereby  altered,  the  integral  representations  becoming 


(1,7b) 


=  (,rg?>0)  (1>8b) 

0 

It  is  readily  seen  that  (1,5)  and  (1,4)  are  equivalent  to  (1,7) 
and  (1,8);  in  fact,  if  in 

y<«)=4f*,0-)-jr&r;  »=w 

we  replace  K  by  its  integral  representation  (/57,  p.  206) 


<y>* 


_ C  cot  yt 


at, 


and  alter  the  order  of  integration,  we  have 


j  TtTP"  |  "7+*  dy' 


But 
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Re  *  >  0; 
Res<0; 


and 


ir  Relf  >0; 
if  Re  if  <  0. 


Tha  aubatitutlena 


give  ua  that 


Reif>0, 
tftxs-  Xwhere  Rflf<0 


which  aakes  it  obvioua  that  (1,3a)  is  equivalent  to  (1,7).  The 
equivalence  of  (1,4a)  to  (1,3a)  follows  froa  (0,6)  and  the  equiv- 
alenee  of  (1,3)  and  (1,7) *• 

2.  Aarnptotlo  Expansions 

Asymptotic  expressions  for  the  0  follow  froa  the  asympto¬ 
tic  expansion  for  Struve  funetloas  (/V ,  p.  163)  t 


H,<s)  =  K,(i)  + 


/•V 

(t; 


*— la- 1 


+  0(*-*  ‘). 


•(1,7)  and  (1,8)  iaply  for  Q*  f  *  H/2  that 

♦  **•  to/** 

Oj(l**)  -  °i  (5^*)  -  **1 

In  whioh  the  bars  denote  ooaplex  conjugates, 
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in  which  Jarg  zj  <  7T.  Putting  z  *  iq  and  *  0  and  -1,  we 

hava 


m-o  f 

CL.flf)  =2  a* 1  +0(g-»-«). 


(2,1) 


(2,2) 


Tha  asymptotic  expressions  are  readily  derived  from  the  integral 
representations  if  we  observe  that  the  main  contribution  for 
large  Be(q2)  comes  from  parts  of  the  integration  path  for  which 
the  variable  of  integration  is  much  less  than  q.  We  expand  the 
(x2  -  q2)“*  of  (1,3)  and  the  1/(y2  -  q2)  of  (1,7)  as  power 
series  in  the  variable  and  integrate  the  series  term  by  term  to 
get 

o.(i,)=!2  + **W*> + o  (»-*•')  I 

•  -0 

0_,  ((,)  -  + 1 2.AT,  (, )  +  ; 

w-o 

HO— g(a^pL_/WftW+o(#  *  '); 

«>l 

r  (,> = ia- + '^7  (<n + o(,-»  •), 


in  which  takes  the  following  values:  (MacDonald  function): 

v  =  0  v.h6S*e  argf  <  0, 
v  =  I  2  "  arg?=0, 
v=l  11  arg?>0. 

More  convenient  approximate  expressions  that  converge  for  all 
values  of  q  may  be  derived  as  follows.  The  representations 


(2,1)  and  (2,2)  ware  deduced  via  integrals  of  the  type 


«r 

P(a)  =  j  p(x)/(jf  +  a)fa 

in  which  p(x)  is  a  weighting  function  that  decreases  rapidly  as 
<  increases;  it  is  exp(-x)  for  (1,3)  and  (1,4),  KqCx)  for  (1,7), 
and  xK^(x)  for  (1,8),  Approximate  values  for  such  integrals  aay 
be  obtained  by  replacing  f(xP  ♦  a)  by  a  few  of  its  early  teras 
in  a  series  around  the  aero  point;  but  this  is  not  the  best 
approximation,  for  more  accurate  expressions  are  obtained  if 
instead  of  n  tsras  of  a  simple  power  series  for  f  we  use  a  poly¬ 
nomial  of  the  same  degree  that  provides  quadratures  of  the 
highest  degree  of  accuracy  /V*  Then  the  approxiaate  value  of 
the  integral  is  determined  by  the  quadrature  sum 

e 

F(a)s*£Ak/{*t  +  o), 

in  which 

A  —  *•  1 

*  ‘ 

eg  being  the  coefficient  to  the  term  of  highest  degree  in  the 
orthonoraalised  polynomial  of  degree  n,  P„(x)  *  a^xa  ♦  bjjX11-1  + 

•  ••••;  and  the  x^  being  the  roots  of  Pn(x), 

The  larger  a,  the  more  precise  the  value  given  by  this 
sum,  and  so  the  expression  may  be  considered  as  asymptotic. 

If  f  is  f(x  ♦  a)  and  p(x)  »  xmexp(-x)  Cm  *  1,  2,  3,  ...)i 
the  polynomials  we  require  are  the  Tschebyscheff-Laguerre  poly¬ 
nomials,  so  ths  coefficients  in  the  sum  and  the  roots  are  known* 
If  p  A  1  or  if  a  is  not  an  integer,  we  require  a  new  system  of 
orthogonal  polynomials  containing  powers  of  x*1* 

If 

we  again  need  to  construct  a  suitable  system  of  orthogonal  poly¬ 
nomials,  whose  roots  must  be  found;  the  coefficients  in  the  sum 
must  fdao  be  calculated.  The  methods  to  be  used  here  are 
standard. 

The  components  of -the  electromagnetic  field  in  the  far 
aone  involve  integrals  whose  weighting  functions  are  (2/ir)K0(y), 

.  (2y/n*)K1(y) ,  and  (2yVn)Kj(y) ;  there  is  also  the  function 
f(y2)  ■  i/(y2  -  q2)  to  consider.  I 


t 
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The  root*  of  polynomials  in  y2  that  ara  orthogonal  on 
tha  half-axis  (0,w)  and  the  coefficients  of  the  quadrature  sums 
for  integrals  of  the  form 


(y)y*f{y*)dy 


are  given  in  the  table . 


• 

y, 

• 

y . 

.  • 

A, 

K*(y) 

i 

1 

,y|  0.483094 

»y*  =  36.316006 

At  0,967671 

,4,  »  0.0123286 

K  (y)  i 

3 

1  • 

1 

,ri-  1,972244 
,y£  -  38.027756 

,4,-0.9714952 

y4j  -  04)286048 

a r.ryiy* 

15 

1 

,y|«=  94)21970 
.X5  *  74.078030 

,4,  -  0,9208498 

,4,  —  0,0791512 

3.  Power  Series 

The  usual  definitions  of  cylindrical  functions  /J7 
readily  show  that  Q0(z)  and  G^U)  are  defined  by  power  series: 

0.W  -£<-!>•  (■»  *  +  C  -  |  ‘  )}  ; 

|  ^  _  —  1 

^j(—  t)"1!  (2m  ~  l)ri(l*i -*  i)W  21"  *(I«  — 1)!«! 

m-l  ' 

■- 1  i 

x(ln*  +  C-^-]g-r)}' 


C  =  —  0,11593152. 


| For  *  «  iq  we  have 


3<» 


0 * = 2 {kI-TW  - ('"  f  + 1 T  +  c‘  -  2  -i)j ; 
0|  “  “  7  “  +  ?siC-ii!rt' 

x(l««+lJ+c'-sr-Sr)}- 

k  - 1  > 


•ad  for  F(q)  and  f  (q) 


^  *  *  +  +  ?^^(,n^  +  /T  +  c“-SF)}‘* 

^  (9 )  =  1  +  -j  +  S  {  (im  -r  l)C«  -  D«  +  *3?»=^ -  |)!«t 

x(taf+l^  +  C-aL-2i)|. 

fc-i  l 


The  functions  have  many  values  (on  account  of  the  presence 
of  the  logarithm  of  the  independent  variable) ,  which  difficulty 
can  be  removed  by  assuming  that  the  z  plane  is  cut  off  along  the  ■ 
negative  real  axis. 

4.  field  Components 

The  dipole  is  magnetic  and  vertical;  it  lies  at  the 
exposed  surface  of  an  N-layer  horizontal  structure  in  which 
only  the  base  is  nonoonducting  /  0  for  l  *  1,  2, 

*  -.1*.  «N  ■  0) .  The  field  components  are  expressed  in  terms 
of  #(q)  and  F*(q)  as  follows  /l7 ; 


in  which 


* - 

hr  =  —  t<tol?'F,(q), 

. 


3,5 


;.=**?+ w-ofl 


j_  JL 

1-2  N*t 


I- 1 


the  effective  impedance  at  the  top  edge  of  the  first  layer,  in 
the  N-layer  structure  is 


r,  =  cth  [%,</,  -f  arcth  j/j  ^z,] , 


and  the  corresponding  impedances  for  the  other  layers  are  given 
by  analogous  expressions.  The  apparent  resistivities  are 
given  by 


in  which 


Pi  **  /  r  V 

ft.  _  h,  in'( ArV-^-'i5  A 

v  =  »7i;=,wUj  U/*«; 

P"‘  _  2x  /  i|  \~  1 


]f ll- 

f  Pi 


(*)(*)• 


>r» 


ef*.,=  -p-;  x,  =  2,81 1 7*7  km  *; 


h  -  18  . 

«/  h.  R  —  -rt- : 

itjr* 


•f  La  c/s. 


-1  -1 

hrH  B  =  6/x^r,  and  ar ^  is  in  ohm  .m  . 

Tables  of  and  ^  as  functions  of  \/d1  and  of  ?  and  f*  as 
functions  of  q  enable  one  to  construct  frequency-sounding 
curves  for  the  far  zone  by  means  of  elementary  calculations. 

The  figure  shows  the  apparent  resistivity  as  deduced 
from  the  electric  field  of  a  vertical  magnetic  dipole  at  the 
surface  of  a  three-layer  structure  whose  parameters  are 


ft.  1  d.  _  n 

T*  dT  =  3’  P*  =  0°- 


The  full  lines  represent  curves  derived  from  asymptotic 
expressions  for  the  far  zone;  the  broken  lines,  curves  for 


£ 
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the  S  cone;  and  the  dot-and-dash  lines,  wave  curves.  The 
points  represent  \«lues  derived  from  quadratures  of  the  integral 
solutions,  and  the  encircled  figures  give  the  ratio  of  the 
distance  to  the  depth  of  the  base. 

Numerical  integration  has  been  used  to  test  the  ampli¬ 
tude  and  phase  curves,  and  also  the  wave  curves  and  the  curves 


IT' 


J 


3*7 


for  the  S  zone;  it  appe&ra  that  tha  asymptote  for  the  far  zone 
gives  much  better  approximation  than  that  for  the  wave  zone  at 
high  frequencies  and  than  that  for  the  S  zone  for  low  fre¬ 
quencies,  the  range  being  that  used  in  frequency  sounding.  The 
frequency-sounding  curves  for  r  =  6i>  and  r  «  9D  show  that  the 
amplitude  and  phase  curves  for  the  magnetic-source  electric- 
detector  method  may  be  calculated  from  the  formulas  for  the  far 
zone  provided  that  the  distance  from  source  to  detector  is  more 
than  six  times  the  distance  down  to  the  insulating  base.  If  a 
magnetic  detector  is  used,  the  distances  must  be  much  larger 
before  the  asymptote  for  the  far  zone  gives  acceptable  accuracy. 
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SOME  FUNCTIONS  OF  THE  COMPLEX  ARGUMENT  WHICH  APPEAR  IN  THE 
LOW-FREQUENCY  FIELD  THEORY 

L.  B.  Gasanenko,  G.  P.  Sholpo,  and  E.  I.  Terekhin  \ 

pp.  78-109 

Problems  of  the  fields  produced  by  a  low-frequency  source 
often  have  integral  solutions  that  can  be  expressed  via  func¬ 
tions  G0  and  Gi  of  complex  argument  (e.g. ,  a  vertical  magnetic 
dipole  in  air  ,  an  infinitely  long  cable  A/»  remote  dipole 
sources  at  the  surface  of  a  horizontally  bedded  structure  itf). 
Bursian  remarked  on  the  value  of  tables  of  these  long  ago  £\J , 
but  they  appear  not  yet  to  have  been  published. 

Recent  rapid  developments  in  the  geophysical  uses  of  low- 
frequency  fields  have  made  such  tables  vital  in  geophysics;  but 
specialists  dealing  with  the  publication  of  tabulated  functions 
have  been  forced  to  concentrate  on  tables  of  greater  use  in 
mathematics.  For  this  reason,  workers  in  the  Geoelectricity 
Laboratory  of  the  Physics  Research  Institute  at  Leningrad  Uni¬ 
versity  in  1959  decided  to  tabulate  these  functions  and  related 
ones. 

The  functions  Gn(z)  are  A7  defined  by 

0*  (*)  as  y  (Hu  (*)  Yn  (*)l  • 


The  properties  of  Struve  and  Weber  functions  give  rise  to  the 
following  properties  of  G„U>! 
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g  (in  f  I  >  !ll  * 

On(P,  —f)  =  0*  to,  9h 


in  which  Q,  (f,f)  denotes  a  function  of  z  --  and  'Gn(/*,-J^ 

denotes  a  function  of  "Z  =^*xp(  -if: 

^"to,  f) — R* 0H (f,  yJ+ZIraO*^,  ?); 

°h  to,  f)  =  Re  Onto.  f)-nm  om  to,  t). 

The  relation  between  Hn(z5  and  Ln(z)  is  /" 2/ 


0.  (iz)  =  7  l<  L*  (z)  —  K,  (/*)) ; 

0i(k)=-J[M*)+K,  (<*)!• 


The  functions  F(z)  and  5t,(z)  that  are  used  to 
components  in  the  far  zone  of  a  low-frequency 
given  by 


express  the  field 
dipole  /%/  are 


The  following  tables*  are  for  Go(iz)  and  G-^(iz)  (Table  1),  and 
also  for  z2f'(z)  (Table  2)  as  functions  of  z  =  fexp(i^)  for.  the 
following  values  of  the  modulus  end  phase  of  the  independent 

variable : 


P  =  0(0,1) 6;  f=±0(5®)45#; 


00 ,  (lz)  =  Re  (1 \z )  -f  i  lm  Ool  (/*); 

in  ’,',lch _ '  ,-P(*)=U»P(»)|<“.  __  _ 

*'.Ve  were  assisted  in  the  wort  of  tabulation  by  i,,  h . 
O'yachenko,  M.  A.  Skorman,  0.  A.  Lherdy  •  ,avi,  and  1..  II.  holot- 
kova.  The  calculations  were  performed  at  Leningrad  hech-uii.-i.-d 
Accounting  .Vorks  with  desk  calcui. -icin  •  •.•chines. 
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Table  1 
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Future  plans  include  the  publication  of  tables  of  Struve 

functions  of  ordsra  zero  and  one  for  f  -  0  (0,01)  and  0  (0.1)  6, 

with  tf  •  0  (5*)  90*  |  and  also  G.(iz)  and  G,  (iz)  for  /r  0  (0.1) 

1,  tO  (5*)  90*  and  £*  0  (0.1)  6,  145*  (5*)  90*. 
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APPROXIMATE  FORMULA  FOR  A  DIPOLE  FIELD  ABOVE  A  HORIZONTAL 
LAMELLAR  STRUCTURE  WITH  A  POORLY  CONDUCTING  BASE 
pp.  110-113  L.  B.  Gasanenko 

Simple  asymptotic  expressions  are  available  for  the  wave 
sons  of  a  low-frequency  dipole  the  scale  distances  r/d^ 

(relative  to  the  depth  of  the  base)  used  in  frequency  sounding 
are  usually  less  than  ten,  and  for  these  they  give  a  fairly  sa¬ 
tisfactory  description  of  the  general  trend  if  the  base  is  a 
nonconductor  or  has  a  conductivity  much  higher  than  that  of  the 
overlying  layers.  However,  if«^/er««l,  the  wave-zone  curves 
are  very  different  from  those  found  by  numerical  solution,  as 
Fig.  1  shows  by  reference  to  ft,  for  w2/w^  =  y  of  1/9  and  l/39« 
The  wave  curve  deviates  from  the  exact  one  because  the  wave  ex¬ 
pressions  are  deduced  without  allowance  for  the  energy  propating 
in  the  ground;  formally  speaking,  this  amounts  to  neglecting  m 
in  expressions  containing  nt  =  (or  +  x^)^.  For  example,  the 
f(m)  of  the  integral  solution 


e? = r%  j /,  (mr)  2f(m)  m'dm; 

/("I=^~7T;  "*■  =  ' '«■[*.*,  +  Arcd. 

is  expanded  as  a  power  series  in  m  on  the  assumption  that  the 
region  of  integration  that  introduces  the  largest  contribution 
is  such  that  the  moduli  of  the  wave  numbers  for  all  layers  (in¬ 
cluding  the  base)  are  much  greater  than  m.  But  if  1 

(X  =  1,  2,  ...,  N  -  1),  direct  calculations  for  the  region  ha¬ 
ving  Xjr  **  8  and  x^r  <  k  show  that  the  part  with  Ix^l  <.  m  |xj| 
makes  the  largest  contribution. 

For  this  reason,  we  derive  an  expression  that  accounts 
approximately  for  the  energy  propagating  in  the  base;  to  cal- 
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Fig.  1.  Field  of  a  vertical  magnetic  dipole  over  a 
stratum  having  a  weakly  conducting  base:  a)  rela¬ 
tion  of  the  modulus  of  the  apparent  resistance  (as 
determined  from  the  tangential  component)  to  wave¬ 
length  for  «|/*7  =  1/9  and  r/d^  =  9;  b)  the  same 
for  ' ®  1/39  and  r/d^  =  9*  c)  the  referred  ver¬ 

tical  impedance  for  ®^/«»T  =1/9  and  r/d^  of  h,  6, 
and  9t  1)  exact  curve  (from  quadratures  of  inte¬ 
gral  solutions);  ?.)  wave  curve;  3)  approximate 
curve.  The  encircled  numbers  are  the  r/d-^. 

eulate  the  derivatives  we  assume  that 


«/ *^0;  <  1  <  /  <  A/) ; 

1  »* 

si  -f-  -j-  0; 


Then  the  apparent  resistance  is  found  as 
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in  which  all  0z_  relate  to  a  structure  with  a  nonconducting 
base.  In  particular,  for  a  layer  above  a  poorly  conducting 

base , 


j_  f  •  0*1 — cth 


Approximate  values  of  the  field  components  are  given  by 


hr 


6 

*,r  »*>’ 


and  the  reduced  vertical  impedance  is 


T Z*  ^  : 

The  wave,  approximate,  and  numerical  curves  are  compared  in  Fig, 
1;  if  y  =  1/9  (Fig.  la),  the  IyCu>I  Given  by  the  approximate 
formula  for  5*5  <  \j/d  <.  11  is  in  good  agreement  with  the  quadra 
ture  curve  for  r/d^  =  9.  Again,  if  y  =  1/39  (Fig.  lb),  the 
approximate  curve  lies  closer  to  the  exact  curve  than  it  does  to 
the  wave  curve.  IVe  may  say  that  for  5  <■  A]/dj.  <c  10  the  approxi 
mate  curve  represents  the  behavior  above  a  poorly  conducting 
base  as  well  as  (or  as  poorly  as)  the  wave  curve  represents  the 
behavior  above  an  insulating  base. 

Figure  lc  shows  the  vertical  field  impedance  for  a  verti¬ 
cal  magnetic  dipole  above  a  two-layer  structure  for  which  = 
•t/9*  The  approximate  curve  is  close  to  the  |z||  curve  for 
r/&i  =  6.  The  shape  is  almost  precisely  that  of  the  exact 
curve  for  r/d^  =  9,  but  the  approximate  curve  is  displaced,  the 
relation  between  the  two  being  roughly 
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The  Approximate  formula  is  applicable  for  general  purposes  at 
sodium  wavelengths  and  at  medium  values  of  r/dp. 

More  accurate  expressions  could  be  derived  by  expanding 
tho  formal  solutions  as  power  series  in  (m  -  mQ),  in  which  m0  is 
a  funotion  of  X,  r/ dj.,  and  the  structure  of  the  section. 
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DIRECTIONAL  ELECTROMAGNETIC  SOUNDING 
pp.  118-128  q.  V.  Molocbnov 

Induction  methods  are  now  being  extensively  used  in 
deteraining  the  depths  of  high-conductivity  bodies  covered  by 
low-conductivity  rocks.  The  physical  principles  are  the  saae, 
no  aatter  what  apparatus  may  be  used;  a  primary  electro- 
aagnetic  field  (perhaps  from  several  sources)  produces  eddy 
currents  in  the  conducting  medium,  and  these  produce  their  own 
field,  which  is  called  secondary.  The  structure  of  the 
primary  field  is  governed  by  the  shape  and  size  of  the  source 
(usually  a  fraae  bearing  several  turns  of  wire  connected  to  a 
generator),  while  that  of  the  secondary  field  is  governed  by 
the  source,  the  frequency,  the  shape  and  size  of  the  conductor, 
the  distance  froa  the  source,  and  the  electromagnetic  para¬ 
meters  of  the  aaterial.  The  secondary  field  is  thus  the  only 
one  of  geophysical  interest. 

It  is  usual  to  examine  the  magnetic  component  by  means 
of  a  frame  probe;  measurements  are  made  of  the  amplitude  and 
phase  for  the  components  along  the  coordinate  axes  or  of  the 
elements  of  the  polarization  ellipse  (mainly  the  inclination 
of  the  major  axis).  Such  methods  are  used  at  ground  level, 
underground,  and  in  aircraft. 

The  source  coil  is  set  either  horizontal  or  vertical  in 
current  methods.  The  secondary  field  becomes  appreciable  only 
at  some  distance  froa  the  source,  so  it  is  usual  to  make  the 
distance  froa  source  to  body  several  times  the  depth  of  the 
body  below  the  surface.  This  is  an  essential  feature  of  all 
electrical  methods,  and  it  requires  an  operator  at  the  detect¬ 
ing  end  as  well  as  one  at  the  source  end  with  some  means  of 
communication  between  them;  but  it  makes  it  difficult  to 
measure  the  phase  of  the  field,  although  this  can  provide 
information  on  the  structure. 

The  secondary  field  is  always  observed  together  with  the 
primary  field,  the  latter  forming  a  background  that  increases 
the  error  of  the  result.  Attempts  have  been  made  to  balance 
out  some  or  all  of  the  primary  field,  but  these  have  not  found 


n 


application  in  practical  field  work.  1 

In  principle,  it  would  be  possible  to  operate  induction 
aethods  froa  moving  vehicles,  but  this  is  made  difficult  by  the 
peed  to  have  source  and  detector  a  considerable  distance  apart* 
For  the  same  reason,  the  source  and  detector  must  be  in  differ¬ 
ent  shafts  or  boreholes  when  the  method  is  used  underground,  and 
even  then  the  region  surveyed  is  only  a  small  part  of  that 
between  source  and  detector. 

X  have  developed  a  method  of  directional  electromagnetic 
iopnding  that  is  free  from  these  disadvantages.  The  directional 
sensitivity  with  respect  to  the  conducting  object  is  available 

for  any  frequency  and  any  shape  of 
object.  The  method  is  essentially  as 
follows.  An  alternating  current  is 
passed  through  a  coil  wound  on  a  frame 
(the  source  coil);  a  second  coil  (the 
receiving  coil)  is  rigidly  set  up 
withip  (or  near)  the  source  coil,  the 
planes  of  the  two  being  perpendicular. 
Then  the  primary  field  of  the  source 
is  not  detected  by  the  receiving  coil,, 
whereas  the  secondary  field  will  pro¬ 
duce  an  emf  in  the  receiving  coil  if 
the  system  is  suitably  placed.  The 

„  _  ..  .  ..  theory  is  as  follows.  An  inclined 

Fig.  1.  Inclined  magnetic  ..  w 

®  magnetic  dipole  has  a  moment  M  = 

dipole  above  a  two-layer  „  „  , _ ,  J 


structure. 


M0exp(i6rt),  in  which  M0  is  amplitude, 
Oi  -  rf ,  f  is  frequency,  and  t  is 
time.  (The  time  factor  and  the  subscript  zero  are  henceforth 
ommitted.)  The  xOy  plane  is  that  of  the  surface  of  the  Earth, 
the  Os  axis  bsing  vertical.  The  dipole  lies  ,at  the  point 
(O,’0,hQ),  the  projection  on  the  xOy  plane  coinciding  with  the 
Ox  axis.  The  angle  between  M  and  the  Ox  axis  is  ^  (Fig.  1). 

Wq  take  the  simple  case  of  two  homogeneous  and  isotropic 
layers  having  a  horizontal  interface;  let  and  h-j  be  the 
conductance  and  thickness  of  the  top  layer,  the  corresponding 
quantities  for  the  lower  layer  being  and  h%  (=«»).  The 
aagne^if  permeability  is  taken  as  unity;  the  displacement 
currents  in  all  media  are  neglected.  The  conductance  of  the 
air,  0*0,  is  taken  as  zero. 

The  field  of  an  inclined  dipole  can  be  considered  as 
made  up*of  the  field  of  a  Vertical  dipole  Mg  =  Mainland  of  the 
field  of  a  horizontal  dipole  MP  =  Mcos  Standard  solutions 

for  the  fields  of  vertical  /I/  and  horizontal  /2/  dipoles  f<Jr 
points  on  the  line  x  =  0,  z 
•nts  of  a  vertical  dipole 


h0  give  us  for  the  field  compon- 


mm 

tf*  =  —7T  (nr)  an'dn. 


(1) 
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Hy — M,j  /,  (nr)an*dn, 

Hx=0; 

and  of  a  horlsontal  dipole 

m 

Hx=-~t^~yt(nr)  Mdn, 

H,  =  HgLo, 


(2) 


(3) 


in  which  r  is  the  distance  to  the  point  of  observation,  flit 
J0  ^(nr)  are  Bessel  functions  of  the  first  kind  and  of  the  sero 
an&  first  orders, 


'«  = 


i  +  «**!*' 


*3T 


(4) 


•*“■5+3-*  «?-<*$ ft. 

and  c  »  3  *  10 10  ca/sec.  The  first  terns  in  (1)  and  (3) 
represent  the  primary  field;  they  become  infinite  when  r  =  0. 
The  second  terms  in  (1)  and  (3),  and  the  Hy  of  (2),  represent 
secondary  fields  resulting  from  conduction.  Now  J-jCnr)  ■*  nr/2 
when  r-*0,  and  JQ(nr)*»1,  so  Hy-#»0  in  (2),  but  the  second 
terms  in  (1)  and  (3)  remain  finite.  The  total  fields  of  (1) 
and  (3)  tend  to  infinity  as  r-*»0,  so  the  conducting  medium 
has  relatively  little  effect  at  short  distances. 

An  inclined  dipole  gives  at  the  same  points 


H, = M  cm  *  |  -  ±-  +  I  jj,  (nr)  m «/«  J , 

ss 

Hy^M  sin  Jl  (nr)  xn'dn, 

Ht ss M sin <|» ~  +  Jj0 (nr) *n*dn  J. 


(5) 


(6) 

(7) 


The  field  strength  normal  to  the  moment  and  the  Oy  axis  is 


(apart  from  sign) 


in  which 


Hxg  —  Hxt\n\  —  Hgcosty  **  -y-*in2<j» jjy4(/ir)a«a<fc  -  (8) 

**  V 

—  T  J J\  (nr)  a  ndn  |  —  M  sin  2*  F(r), 

^(r) =  \-\^J<>(nr)*n%dn — ~“^J%{nr)  <x  ndn  |.  (9) 


Tha  other  component  normal  to  the  dipole  lies  along  the  Oy  axis 
it  is  given  by  (6).  These  components  of  (6)  and  (o)  are  deter 
mined  only  by  the  secondary  field;  the  Hy  of  (6)  is  caused  by 
the  vertical  component  of  the  magnetic  field,  so  it  is  largest 
for  ^  «  90°.  I  do  not  consider  Hy  here,  for  this  is  discussed 
in  existing  papers.  As  regards  the  Hxz  of  (b),  the  principal 
feature  is  that  Rxz  remains  finite  as  r-**0,  being 


^jrjr— 1/4  Af  sin 


an'dn. 


(10) 


Another  important  feature 


Fig.  2.  Polar  diagram 
of  the  apparatus. 

numerical  integration  for 
be  given  in  approximate  fo 
especially  one  consisting 
on  a  base  whose  conductivi 
then 


is  that  Hxz  varies  as  sin  2^,  so 
HXz  s  0  if  ^  is  0  or  90°  (the 
secondary  field  is  not  detected  in 
these  positions).  The  maximum  Hxz 
occurs  when  =  45°,  the  value  being 
governed  by  F(r).  Figure  2  shows 
the  polar  diagram,  the  two  perpen¬ 
dicular  lines  being  the  projections 
of  the  planes  containing  the  source 
and  receiving  coils.  The  points 
represent  measurements  made  in  air 
above  a  sheet  of  aluminium;  here 
the  coils  were  10  cm,  in  diameter 
(source)  and  8  cm  (receiver).  In 
what  follows,  subscript  xz  is 
omitted,  so  HxZ  =  H.  The  integrals 
in  F(r)  can  be  calculated  only  by 
any  medium  generally;  but  they  can 
rm  for  certain  layered  media, 
of  an  insulator 'of  thickness  h^  lying 
ty  is  <7^  =  <r.  Let  hQ  +  h^  =  h  here; 
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(11) 


a-Zlir+jT 

« +  V'V  T»* 


e-w. 


As  asymptotic  series  in  negative  powers  of  Z  is  obtained  for  F(r) 
if  { Z )  «  2hjxj  =0  1: 


f(r) 


1 

frIF 


4-  3 


'-T** 


*('  +  W 


! -/?»  I 

TT‘T 


(12) 
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in  which  R  =  r/h.  The  subsequent  terms  contain  odd  powers  of 
Z;  only  the  first  terms  need  be  considered  if  Z  is  large,  the 

error  then  not  exceeding  the  abso- 
lute  value  of  the  first  discarded 
term. 


Fig.  3«  Relation  of  to 
distance  for  an  inclined 
dipole  above  an  ideally 
oonducting  half-space. 


Consider  the  first  term 
within  the  bracket  in  (12),  which 
we  denote  by  Y;  this  corresponds 
to  a  perfect  conductor.  Now  Y  * 
for  R  =  0,  and  Y  at  first 
decreases  as  R  increases,  becoming 
zero  when  R  =  V2.  The  maximum 
negative  value  of  -0.101  occurs 
when  R  =  2-f}/2 ;  thereafter  the 
value  tends  monotonically  to  zero 
(Fig.  3)»  Measurements  (shown  by 
points)  have  been  made  above  a 
sheet  of  aluminittm  with  coils  3  cm 
in  diameter.  In  this  case  H  is 
maximal  when  r  =  0,  so  the  receiv¬ 
ing  coil  must  be  placed  at  r  *  Oj 
then  for  <J>  =  45°  we  have  that 


H 


M 


(13) 


For 


finite  conductivity 


6  ,  24  fO  .  630  \ 

?  +  7f—  zr  +  zr*"*--]- 


(14) 


|  The  number  of  terms  needed  in  (14)  for  a  given  Z  can  be 
>  > 
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determined  if  the  permissible  error  in  H  is  specified.  The  x 
that  appears  in  Z  is  best  expressed  in  terms  of  the  wavelength 
within  a  layer: 


x  =  i  VT2*  V2 3/=  Vi  2*  V2  } , 


in  which 


X  = 


/  8/' 


(15) 


Then  (14)  can  be  put  as 


(16) 


in  which  denotes  the  expression  within  braces.  Table  1 
gives  the  Re,  Im,  a  (modulus),  and  (phase)  of  <|>  for  several 
J/h  and  for  several  termination  points.  The  error  in  the 
amplitude  is  only  0.7%  for  A/h  =  0.75,  1.5%  for  A/h  =  1,  and 
5%  for  Vh  =  1.5  if  the  first  two  terms  in  (16)  are  taken,  so 
these  terms  are  sufficient  if  >/h  is  small. 

The  converse  problem  (determination  of  h  from  known  M 
and  H)  is  usually  solved  by  frequency  or  radial  sounding,  the 
resulting  curves  being  compared  with  theoretical  ones.  This 
method  is  the  usual  one  and  is  widely  used  in  electrical  surveys. 
The  interpretation  can  be  made  very  much  easier  if  we  determine 
the  effective  penetration  depth  .  Here  we  can  also  use 
(13)  to  determine  the  depth,  tmd  from  this  subtract  the 
effective  penetration  depth  £ih,  which  is  deduced  as  follows. 

Let  H  be  the  field  when  the  conductivity  is  infinite;  then 
(13)  gives 

*=tj/  2 wr-  (i?> 

(17)  will  give  a  depth  D  =  h  +  flh  when  the  conductivity  is 
finite,  and  6h  is  given  by 


(18) 

dH 

in  which  'Is  the  real  part  of  the  second  term  in  (l6).  Then 
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Table  1 


Ne.  of 
vttms 

X 

* 

Ke 

Im 

a 

? 

1 

1 

0 

1 

0 

2 

0,8807 

0,1198 

0,890 

7  45' 

3 

0,5 

0.8807 

0.1008 

0388 

630 

4 

03816 

0.1012 

0390 

635 

6 

0,8816 

0,1012 

0390 

635 

l 

1 

0 

1 

0 

2 

0,821 

0.179 

0,84 

3 

0,75 

0.821 

0,1363 

0,834 

9  25 

4 

0324 

0,1395 

0.846 

9  35 

3 

0,824 

0,1396 

0,846 

9  35 

1 

1 

0 

1 

0 

2 

0,7613 

0,2387 

0,798 

17*25' 

3 

1.0 

0,7613 

0,1627 

0.779 

12  05 

4 

0,7689 

0,1703 

0,786 

12  30 

5 

0,7687 

0,1705 

0.786 

12  30 

1 

1 

0 

1 

0 

2 

0342 

0358 

0,735 

29*10' 

3  ' 

1.5 

0,642 

0,187 

0,670 

16  15 

4 

0,6676 

0,2126 

0,700 

17  40 

5 

1 

0,6657 

0,2145 

0,700 

17  50 

1 

I  1 

0 

1 

0 

2 

0,5226 

0,4774 

0,707 

42  45' 

3 

2,0 

03226 

0,1734 

0,550 

18  20 

4 

0,5831 

0,2309 

0325 

21  35 

o 

• 

0,575 

0,2389 

0,623 

22  35 

which  ie  the  value  giveu  previously  C$J •  Table  2  gives  values 
calculated  from  (18)  by  means  of  Table  1;  for  simplicity  M  has 
been  put  as  16  and  H  as  1,  so  h  =  1.  The  results  are  clearly 
satisfactory  even  for  >»/h  *  2. 

It  has  been  shown  C*>7  that  &h  can  be  determined  by  means 
of  measurements  at  two  frequencies  f^  and  f 

-  -£*c£l. 

j/4-I  ’  (20) 


8l 


Table  2 


1 

I 

D  (n) 

AM  -fro**  01) 

D-  AH 

Er»or, 

% 

IJOOO 

0,000 

1,000 

0,0 

OJS 

1.040 

0,040 

1.00 

00 

1,76 

1,068 

0,058 

1,000 

02 

1,0 

1,082 

0,08 

1.002 

0,2 

1,125 

0,119 

1,006 

0.6 

2.0 

(1.17) 

0,150 

(1.011) 

(M) 

1b  which  Df  and  Dp  are  the  depths  given  by  (1?)  for  these 
frequencies  and  is  the  Ah  for  f-j.  Then,  since  Ah  is  now 
known,  we  can  determine  9  from  (19)  as 


°  ~  (4xAA)*/  ' 


(21) 


The  effeots  of  finite  conductivity  in  the  first  layer  on 
the  determination  of  h  can  be  established  as  for  the  dipole 
eleetromagnetio  method  • 

It  can  be  shown  that  the  method  is  applicable  to  bodies 
of  finite  site;  the  case  of  a  sphere  illustrates  this. 

An  ideally  conducting  sphere  of  radius  r0  is  embedded 
in  a  poorly  conducting  medium;  the  center  lies  a  distance  h 
below  the  surface.  A  dipole  inclined  at  0/  to  the  horizontal 
and  having  a  moment  M  is  placed  at  the  epicenter.  We  assume 
that  the  field  of  the  dipole  would  be  homogeneous  throughout 
the  volume  of  the  sphere,  the  contribution  from  the  vertical 
somponent  being 


•jp-Sln^; 


sad  that  from  the  horizontal  component  being 

u  _  M 

"or  —  ""  *7pr  cos  <|». 

This  is  rather  a  gross  simplification,  but  the  general  nature 
Of  the  result  is  correct.  No  more  exact  solution  is  necessary, 
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tor  real  bodies  are  merely  more  or  lees  isometric . 

The  secondary  field  at  the  epicenter  given  by  a 
vertical  dipole  is 


—  2M slmj»  -jj- ; 

and  by  a  horizontal  dipole 

Mr~ — -y  M  CO»  <J>  . 

The  field  strength  in  a  direction  perpendicular  to  the  dipole 
is 

//r3ln«i»  — //,co«t=  4/Ws,n2+'F‘>  (22) 

The  directional  feature  is  retained  here  also. 

The  above  discussion  shows  that  a  directional  equipment 
can  be  operated  at  any  suitable  frequency  if  the  two  coils  are 
mutually  perpendicular  and  have  a  common  center  (they  can  be 
displaced  by  a  distance  small  relative  to  the  depth  of  the  con¬ 
ducting  object,  as  Fig.  3  shows).  The  signal  from  the  receiv¬ 
ing  coil  is  amplified  and  fed  to  a  meter  or  recorder,  which  is 
calibrated  directly  in  field  strength.  The  presence  of  a 
nearby  conductor  can  be  detected  by  rotating  the  coils  about  a 
horizontal  axis  in  various  azimuths  until  the  maximum  reading 
is  obtained;  the  direction  is  determined  rather  more  precisely 
from  the  points  of  silence,  as  the  polar  diagram  (Fig.  2)  shows. 
The  depth  is  determined  from  the  maximum  reading;  the  position 
of  the  upper  surface  can  be  determined  if  the  body  is  suffi¬ 
ciently  large. 

The  apparatus  is  compact,  which  makes  the  method  con¬ 
venient  and  economical  to  use;  the  receiver  and  generator  can 
be  contained  in  a  single  case,  so  the  phase  can  be  measured  as 
well  as  the  amplitude,  and  this  provides  more  detailed  informa¬ 
tion.  The  method  can  be  used  underground,  in  aircraft,  and  in 
moving  land  vehicles.  Moreover,  it  can  be  applied  to  frequency 
sounding  of  layered  structures  of  email  depth.  In  every  case 
the  operations  are  simple  and  the  interpretation  is  more 
complete* 


83 


Bibliography 


1. 


2. 


4. 


5, 


Moloohnov,  a.  V.  Vajbikaj.'^vy  magnitnyv  diool*  _na 

poverkhnosti  gorizobtH aloia toy  struktur. _ Ocoflz- 

vartical  magnetic  dipole  on 
the  surface  of  a  horizontal  layered  structure.  Geo¬ 
physical  Prospecting  Methods].  Moscow,  Gosgeoltekhizdat , 
1955. 

Gasanenko,  L.  B.  and  0.  V.  Molochnov.  Elektromagnitnoye 
pole  gorisontal'nogo  nagnitnogo  dipolya  nad  gorizontal' 
nosloistoy  strukturoy  fthe  electromagnetic  field  of  a 
horizontal  magnetic  '3ipole  over  a  horizontally  layered 
structure].  Ueh.  Zap.  Leningrad.  Gos.  Univ.  [Archives 
of  Leningrad  University],  No.  249,  1953. 

Molochnov,  G.  V.  Dipol'nyy  elektromagnitnyy  me tod 
opredeleniya  glubiny  zaleganiya  sloya  konechnoy  provodi- 
aosti  [A  dipole  electromagnetic  method  of  determining  the 
depth  of  a  layer  of  finite  conductivity].  Vcstnik. 
Leningr.  Un-ta  [Herald  of  Leningrad  University],  No.  10, 
I960. 

Idem.  0  vliyanii  provodimosti  pervogo  sloya  pri 
ppredelenii  glubiny  zaleganiya  provodyashchego  osnovaniya 
dinol'nva  elel^troaagnitnvm  metodoa  [Effects  of  conductivity 
of  the  first  layer  in  the  determination  of  the  depth  of  a 
conducting  base  by  the  dipole  electromagnetic  method]. 

Pch.  Zan.  Leningrad.  Gos.  Univ.  [Archives  of  Leningrad 
University],  No.  2&6,  I960. 

Izd.  Inost.  Lit.,  1954. 


84 


r 


COMPARISON  OF  THE  INDUCTION  METHOD  WITH  THE  DIPOLE 
ELECTROMAGNETIC  METHOD  ON  A  THIN  CONDUCTING  PLATE 
/MODEL) 

G.  V.  Molochhov  and  G.  M.  Grebennikov 


pp.  129-134 


Induction  methods  are  now  much  used  to  detect  steeply  dipping 
strata.  The  dipole  method  /I , ?J  has  given  good  results  for 
the  depths  of  conducting  strata;  here  the  inclination  of  the 
magnetic  vector  is  measured  for  the  field  of  a  vertical  magnetic 
dipole.  The  method  has  been  applied  to  measure  the  thickness 
of  floating  ice  (0.5  to  13  m)  on  the  drifting  stations  'North 
Pole  6*  and  'North  Pole  9*,  as  well  as  along  the  northern  ska 
route  in  1958-61  fbj*  The  error  here  is  less  than  y%. 

Vertical  dipoles  have  been  used  in  experiments  with  con¬ 
ducting  bodies  (plate,  cylinder,  sphere)  in  the  Geoelectricity 
Laboratory  of  the  Physics  Research  Institute,  Leningrad  Univer¬ 
sity  /V*  It  has  been  found  that  the  vertical  component  of 
the  magnetic  field  enables  one  to  determine  the  sise  and 
inclination  of  a  plate  in  horisontal  projection  (i.e.,  the 
shape  and  position  in  space). 


Here  we  compare  the  induction  method  with  the  dipole 
method,  the  source  and  receiving  sections  having  identical 
parameters.  The  inclination  of  the  total  field  vector  to  the 
horisontal  was  measured  in  both  cases.  The  field  was  produced 
by  a  screened  frame  (a  circle  of  diameter  5«3  cm)  fed  via  a 
screened  cable;  the  detecting  coil  had  a  screened  frame  4.7  cm 
in  diameter  and  worked  into  a  tuned  amplifier  feeding  a  MVL-1 
vacuum-tube  voltmeter.  The  model  was  a  brass  sheet  126  x  33 
x  0.1  cm,  which  was  supported  on  an  insulator  and  was  fitted 
with  a  means  of  adjusting  the  inclination  y,  which  was  varied 
from  0°  (vertical)  to  90®  (horisontal)  by  15°  steps.  The  top 
edge  of  the  sheet  was  always  horizontal. 


Measurements  were  made  at  distances  of  20,  30,  W),  50, 

|  and  60  cm,  with  the  top  edge  at  depths  of  10  and  20  cm.  The  _ j 
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Table  1 


r  (cm) 

Actual 

*<«*> 

Celctiilutcd 

A  <«*) 

Ml 

aT. 

7 

5 

4 

3.1 

3.1 

fWw!a 

m«tW 

■free* 

0) 

S|  +  *2 

SO 

30 

40 

80 

40 

14 

5 

0 

11 

18 

15 

i 

8.4 

9.3 

114 

13.7 

15.1 

1 

80 

59 

50 

43 

43 

- 3 - 

19 

— 

r- 

80 

SO 

— 

19.4 

—  ! 

— 

40 

— 

20,1 

— 

— 

model  wee  moved,  the  ooila  being  kept  fixed.  In  the  dipole 
method,  the  sheet  waa  inelined  towards  the  source  or  towards 
the  detector.  Figure  1  shows  profile  curves  for  y  of  0, 

45.  and  90a. 

The  induetion  results  were  worked  up  to  give  the  angles 
as  usual  Table  1  gives  the  depth  of  the  top  edge  h  as 

dsdueed  by  the  parabola  aethod  ft]  and  froa  ft] 

h  —  (VT  +  Ctg’a,  -  l),  ( 1 

in  whieh  x  is  the  abscissa  corresponding  to  the  extreme  angle 
a..  Table  1  also  gives  (Xi  ♦  Xz)/(a-\  *  op),  which  is  used  to 
calculate  B^/Hp,  which  appears  in  formula  (2)  of  ft]: 


*1±*L. 

*i  +  •* 


(2) 


in  which  X-t,  and  Xj.ay  *r*  th#  coordinates  of  the  angle 
curve  to  the  left  and  to  the  right  of  the  origin,  H1  and  Hg 
being  the  strengths  of  the  primary  and  secondary  magnetic 
fields.  We  dee  that  the  parabola  aethod  gives  satisfactory 
results  for  r/hfti  3-4,  while  (1)  is  satisfactory  for  r/h 
2*4, 


The  position  of  the  electric  axis  is  determined  in  the 
usual  way  in  the  induetion  aethod;  in  the  dipole  method,  the 
|  position  of  the  electric  axis  of  a  vertical  plate  is  determined  1 
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Profile  Curves  over  the  Flate  using  induction 
method  (1)  and  .Mi pole  electromagnetic  method  (2 


r 


Fig.  2.  Profile  curves  for  the  dipole  method  with 
1)  source-detector  setting;  2)  detector-source  setting. 


by  the  point  at  which  the  profile  curve  meets  the  axis  of  abscis¬ 
sas.  The  profile  must  be  measured  twice,  as  shown  by  curves  1 
and  2  of  Pig.  2,  if  the  plate  is  inclined.  The  point  half-way 
between  the  points  where  these  curves  meet  the  axis  is  then  the 
position  of  the  electric  axis. 

■  The  depth  of  the  top  edge  aay  be  deduced  from  the  observa¬ 
tion  that  the  distance  l  between  the  two  points  where  the  curve 
meets  the  axis  is  almost  independent  of  y  if  the  plate  is 
inclined  towards  the  detector,  provided  that  y  <  60°,  as  Table 
2  shows.  Figure' 3  shows  l  as  a  function  of  r  (the  distance) 


Table  2 
r 


1 

10 

40 

50 

60 

f 

40 

50 

60 

*  =f  10  CM 

A- 20  c* 

8 

15 

24 

33 

42 

15 

27 

33 

15 

15 

25 

34 

43,5 

15* 

26 

35,5 

30 

15,5 

25 

34,5 

43* 

12* 

25 

34 

45 

14,5 

24 

34,5 

43* 

14 

26 

33 

60 

13,5 

24,5 

34,5 

44* 

— 

22 

32,5 

75 

12 

23,5 

34  ’ 

42* 

— 

20 

31 

90 

9,5 

25 

34,5 

44 

17,5 

30 

for  depths  of  10  and  20  cm;  the  relation  is  a  straight  line 
of  unit  slope  which  meets  the  line  l  s  -5  at  r  *  h.  This  can 
be  used  to  determine  h. 
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Fig. 
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3. 

to 


Relation 

r. 


uaad  In  datailad  work,  especially 


Table  3  gives  the  ratio  of 
tha  aaxiaun  angles  of  inclina¬ 
tion  of  the  magnetic  field  in 
the  two  nethods;  it  ehowa 
that  >  1  if  Y  >15°,  the 

ratio  increasing  with  y.  This 
neans  that  the  dipole  aethod 
is  best  if  y  >15°. 

The  position  of  the  elec¬ 
tric  axis  is  thus  determined 
reliably  in  both  aethoda. 
Induction  given  the  position 
of  the  top  edge  satisfactorily 
only  for  r/h  of  3-4,  but  the 
dipole  method  is  suitable  for 
any  r.  Both  methods  can  be 
since  the  apparatus  is  the 


Table  3 


7 

n 

r 

0 

15 

30 

45 

60 

75 

90 

20 

1 

14 

143 

2.35 

246 

346 

7.0 

30 

04* 

1.13 

141 

143 

2,34 

4.18 

6.8 

10 

40 

046 

046 

1.1 

1.4 

2,28 

3,05 

4.25 

50 

1.12 

143 

146 

1.46 

1,77 

2,11 

2.8 

60 

148 

1.42 

1.43 

1.43 

14 

1,78 

2.47 

40 

04 

14) 

147 

1.71 

2,28 

8,66 

5,33 

90 

50 

047 

14 

I.M 

14 

3,12 

4.0 

64 

60 

04 

1.11 

Ml 

1,73 

24 

3.26 

347 
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REPRESENTATION  OF  OBSERVATIONAL  DATA  OF  ALTERNATING- CURRENT 
PROSPECTING  AS  AN  APPARENT  SPECIFIC  RESISTANCE 

A.  V.  Veshev 

pp.  167-186 

There  are  many  types  of  electrical  prospecting  methods, 
which  use  alternating  and  direct  currents;  there  are  also  many 
types  of  line,  methods  of  examining  the  fields,  and  so  on.  This 
makes  it  of  value  to  consider  which  are  the  best  parameters  to 
use  to  characterize  the  field. 

The  apparent  specific  electrical  resistance  (apparent 
resistance*)  is  commonly  used  in  d.e.  methods,  in  addition  to 
directly  measured  quantities  (potentials  and  gradients).  This 
is  convenient,  for  it  is  independent  of  the  apparatus  and 
method  if  the  medium  is  homogeneous,  though  it  is  in  general 
only  a  function  of  the  field,  being  the  product  of  some  measured 
quantity  (reduced  to  unit  current)  by  a  coefficient  dependent 
on  the  geometry  of  the  system  (in  particular,  on  the  form  of  the 
grounds,  if  these  are  not  point  grounds).  However,  the  function 
is  such  that  it  corresponds  to  the  specific  resistance  if  the 
medium  is  homogeneous;  the  coefficient  is  adjusted  accordingly. 

The  apparent  resistance  is  a  function  of  the  component  and 
method  if  the  medium  is  not  homogeneous,  but  it  is  still  a  con¬ 
venient  parameter,  for  the  effect  of  inhomogeneity  may  be  ex¬ 
pressed  as  a  ratio  to  some  parameter  common  to  all  systems, 
which  parameter  corresponds  to  the  specific  resistance  of  a  homo¬ 
geneous  medium.  This  makes  the  apparent  resistance  a  con¬ 
venient  parameter  for  geologic  mapping  and  for  structure  studies. 

Alternating  currents  are  used  mainly  to  detect  anomalies  in 
conductance,  the  results  being  given  directly  as  the  amplitude 
and  phase  of  the  field,  as  the  angles  of  inclination  of  the  mag¬ 
netic  vector,  and  so  on.  Sometimes,  as  in  ratio  methods,  the 
results  are  presented  merely  as  changes  in  field  parameters. 

This  is  a  rather  narrow  use  of  a.c.  methods,  and  it  makes  it 
difficult  to  use  the  results  for  structure  mapping  and  resis¬ 
tance  measurement;  but  the  methods  are  actually  at  least  as 


•This  abbreviation  will  be  used  in  future. 
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[valuable  as  d.c.  ones  for  thesa  purposes.  Before  they  can  be 
so  used,  though,  there  is  need  for  sobs  work  on  apparatus  and 
oa  methods  of  presenting  the  results.  Here  I  shall  deal  only 
with  the  latter  aspect  of  the  matter. 


It  is  best  to  present  the  results  from  a.c.  methods  as 
reduced  parameters  of  apparent-resistance  type.  The  parameter 
to  be  used  is  to  be  chosen  on  the  basis  that  its  dimensions 
and  numerical  value  must  be  the  same  as  f (the  specific  res¬ 
istance)  for  a  homogeneous  and  isotropic  half-space;  then  an 
inhomogeneous  half-space  gives  a  parameter  analogous  to  the 
speeific  resistance.  Here  we  must  point  out  that  the  ap¬ 
parent  resistance  will  be  dependent  on  the  source  and  on  the 
dimensions  of  the  system  if  the  space  is  inhomogeneous.  It 
is  convenient  to  denote  the  d.c.  and  a.c.  values  of  the  ap¬ 
parent  resistance  by  J  and  ,  by  analogy  with  f ,  the  frequency 
being  indicated  (if  necessary)  by  a  subscript  to  f. 

There  are  great  difficulties  in  calculating  f  for  artifi¬ 
cial  electromagnetic  fields,  for  the  situation  differs  from 
that  in  d.c.  methods,  where  the  potential  and  the  gradient  are 
linearly  related  to  f  and  bear  a  simple  relation  to  the  geo¬ 
metry  of  the  system  (the  coefficient  k  has  the  dimensions  of  a 
length  and  is  defined  Hy  purely  geometrical  relationships). 

The  source  in  an  a.c.  method  is  an  electric  or  magnetic  dipole, 
or  a  very  long  cable;  here  the  vector  potential  and  the  field 
components  have  a  complicated  relation  to  the  reduced  distance 
kr,  in  which  r  is  the  distance  from  the  source  and  k  is  the 
wave-number  of  the  medium.  If  displacement  currents  are 
negligible,  and  if  k2  can  be  treated  as  imaginary,  we  have 

in  which  r  is  in  km,  f  is  in  c/s,  and  f  is  in  ohm-m.  For  field 
aources  of  finite  Else,  we  have  an  equivalent  quantity  : 

§'  =  5/-0,4nr/y^.  (2) 

Formula  (2)  shows  that  the  length  l  (km)  of  the  source  line 
affects  the  field  pattern;  the  effects  of  this  length  are-some- 
times  described  by  a  parameter  P,  which  for  k2  imaginary  is 

given  by 

P = 0,4*/  ]/^,  (3) 

in  which  l  is  in  km,J  and  fare  in  c/s  and  ohm-m,  resp. 

Mow  r  and  f  are  always  known,  so  §  and  5’  give  us  f,  provided 
that  there  is  a  unique  relation  between  the  measured  field 

L  .  J 
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[  components  and  §  or  5'**  The  apparent  reaietance  as  ao  found  , 
enables  one  to  separate  field  changes  associated  with  features 
of  the  section  from  changes  associated  with  the  apparatus.  The 
apparent  resistance  will  be  the  specific  resistance  for  a  homo¬ 
geneous  half-space  no  natter  what  the  system  and  wavelength. 
Moreover,  the  asymptotic  values  will,  for  a  two-layer  structure, 
correspond  with  the  resistances  of  the  two  layers.  The  apparent 
resistances  as  found  by  d.c.  and  a.c.  methods  for  a  given 
section  enable  one  to  conpare  resolving  powers;  the  problea  can 
also  be  discussed  theoretically  to  a  certain  extent,  but  the 
a.c.  fields  nust  then  be  represented  in  terms  of  parameters 
corresponding  to  the  apparent  resistance**. 


The  apparent  resistance  is  a  function  of  many  parameters, 
so  symbols  ere  needed  to  indicate  the  type  of  apparatus  and 
field.  These  may  be  subscripts  to  f  and  such  as  3,  f. 
(horisontal  electric  dipole),  M.B.  (vertical  magnetic  dipole), 
M.r.  (horisontal  magnetic  dipole),  b.k.  (infinite  cable), 
k.k.  (finite  cable),  ns  (vertical  magnetic  component) «  mr 
(radial  nagnetic  component ), if  (azimuthal  electric  component), 
and  so  on.  The  symbols  for  the  source  will  be  placed  in  the 
right  upper  position,  and  those  for  the  field  in  the  right 
lower.  for  example,  pjj;*  and  f  denote  respectively  the 
use  of  the  vertical  magnetic  component  from  an  infinite  cable 
and  of  the  seme  from  a  horisontal  nagnetic  dipole.  Only  the 
electric  field  is  used  in  d.c.  methods,  so  there  is  no  need  to 
put  a  subscript  for  electric  to  f  .  Here  we  can  us-.*  symbols 
such  as  (potential  system),  ft.  (gradient  system),  ft..Cp 
(nean  gradient  system),  and  so  on;  the  superscript  indicates 
the  source  type  (e.g.,  dipole  source).  These  numerous 
subscripts  and  superscripts  are  somewhat  inconvenient;  in  any 
given  case,  they  nay  be  omitted  after  the  type  of  source  and 
so  on  have  been  stated.  This  question  of  notation  requires  a 
separate  study. 


•The  field  source  is  represented  as  a  plane  incident  wave 
when  natural  fields  are  used,  as  in  magnetotelluric  sounding. 
Here  ? *  4«z/y,  in  which  z  »  E/H  and  cc  is  frequency.  The  ratio 
of  any  two  components  will  be  dependent  on  J  for  any  other  form 
of  source,  so  f  can  then  be  deduced  only  from  5  • 

••Van'yan's  &*and  fr  for  frequency  sounding  and  for  field 
equilibration  with  an  electric  dipole  and  analogous 

quantities  for  magnetic  dipoles  /67,  do  not  satisfy  these 
requirements  fully;  the  fa  and  fr  for  a  homogeneous  half-space 
to  f  as  limit  only  for  J  very  small  (S  zone)  and  very 
(wave  zone). 


It  is  possible  to  deduce  f  in  any  given  case  only  if  § 
or  io  known;  the  conditions  for  this  vary  with  the  type  of 
source  and  can  be  deduced  from  the  type  of  variation  of  the 
normal  field.  Aa  an  example  I  consider  the  use  of  magnetic 
dipolea.  Here  there  are  many  different  forms  of  source  and 
detector  /V,5,7“9,1‘1  «127,  but  the  main  differences  relate  only 
to  the  frequency  and  to  the  orientation  of  the  dipoles  with 
roapeet  to  the  surface  and  to  one  another.  The  real  and 
imaginary  parts  are  recorded  for  one  or  more  field  components, 
and  also  the  magnitude  and  inclination  of  the  total  field  (as 
in  induction  methods).  The  normal  fields  of  magnetic  dipoles 
on  the  surface  of  a  homogeneous  conducting  space  have  been 
examined  the  field  components  are  represented  as 
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is  which  M/r^  is  the  field  of  a  dipole  in  air  (sometimes  called 
the  primary  field);  the  upper-case  letters  represent  compon¬ 
ents  of  the  electromagnetic  field,  while  the  lower-case  h  and  e 
demote  changes  resulting  from  the  use  of  a.c.  (these  arise  from 
conductivity)  and  are  dimensionless  quantities  called  the  mag¬ 
netic  and  electric  numbers  (Krayev's  notation). 

It  has  been  shown  that  these  numbers  have  the  following 
relations: 

A?' r = cos  9  (4  —  e*'  *  —  hg'  *), 

a?  r= sm  ?(*;;•• -2),  (5> 

hr — cospAr*. 


In  particular,  (5)  shows  that  the  radial  magnetic  number 
of  a  vertical  magnetic  dipole  is  equal  to  the  vertical  magnetic 
number  of  a  horieontal  magnetic  dipole  if  ~  0.  The  real  and 
imaginary  parts,  and  also  the  amplitudes  and  phases,  have  been 
given  for  the  h  for  5  up  to  20  ^3,^/;  this  range  covers  all 
the  main  practical  cases  for  this  type  of  equipment.  The 
asymptotic  expressions  and  tabulated  values  /3,4/  show  that  the 
magnetic  numbers  are  usually  multivalued  functions  of  5;  lor 
example  In(hjf)  has  a  maximum  at  $  •  1.076  and  a  minimum  at 

4.808,  while  Pe(h£*)  has  several  maxima,  the  first  two  turn- 
j  ing  points  being  at  $  of  2.382  and  8.4413,  and  the  first  zero  j 


at  7.21,  the  aubaequent  turning  pointa  baing  aeparatad  by  intar-  j 
▼ala  of  about  TT.  Again,  tha  Re  for  hr  and  hs  (for  Of  90*) 

have  aaxiaa  at  5  *  4.6,  while  tha  la  have  aaxiaa  it  j*  2.6, 
ainiaa  at  3  -  8.0,  and  saroa  at  ^  *  4.5. 

It  ia  poaaibla  to  dataraina  $  froa  ona  or  othar  of  tha 
aagnatie  nuabara  in  aoat  eaaaa,  provided  that  tha  raal  and 
iaaginary  parta  ara  aaaaurad  aaparataly  (with  dua  ragard  for 
aign).  Such  aaaauraaanta  ara  vary  tadioua  to  parfora  and 
intarprat,  and  tha  aathod  in  thia  fora  ia  applicabla  only  in 
apaeial  eircuaatancaa,  a.g.  ovar  ragiona  known  froa  othar 
evidenee  to  ba  auitabla  for  datailad  examination.  It  will  ba 
ahown  that  a  aiaplar  aathod  of  determining  $  oan  ba  uaad  in 
nearly  all  eaaaa. 

Tha  relation  of  f  to  J  for  a  hoaogeneoua  aadiua  ia 

p-7,89  (-£)*/.  (6) 

hut  raal  conditiona  ara  thoaa  of  an  inhonogeneoua  hqlf-apace,  in 
which  caaa  it  ia  batter  to  uaa  tha  eff active  value  5 ,  which 
corraaponda  to  a  fictitioua  hoaogeneoua  half-apace.  Than  (6) 
with  5  give a  ua  f .  Tha  aathod  of  dataraining  f  can  ba  greatly 
aiaplified  if  tha  field  work  ia  performed  with  diatancea  and 
frequeneiea  auch  that  tha  entire  range  in  f corraaponda  to 
rangea  of  aonotonic  change  in  tha  aagnatie  nuabara.  Tha  firat 


Fig,  1.  Palationahipa  fox'  dipolaa  1)  on  tha  aurfaca; 
at  a  height  of  50  a;  5)  a t-  a  height  of  100  a. 


■ueh  rang*  ia  tha  baat,  for  this  corresponda  to  low  J  and  thus  i 
to  fairly  low  fraquaneiaa.  Figures  1-4  give  curves  for  the 
first  ranges  in  tha  aagnetic  numbers  as  functions  of  #r2f; 
these  are  convenient  for  the  determination  of  The  relation¬ 

ship  is  nearly  rectilinear  over  much  of  the  range.  Electrical 
profiling  is  performed  with  fixed  r  and  f;  the  graphs  enable 
one  to  determine  f from  the  magnetic  number. 

The  practical  use  of  these  relationships  may  be  established 
by  examining  tljj  rangea  in  $ they  cover.  The  table  below  gives 
the  ranges  of  f  that  correspond  to  the  frequencies  used  in  the 
API  and  PEMK  equipments  (75,  125.  375,  1125,  3375,  and  7500  c/s), 
the  r  being  0.05,  0.1,  0.2,  and  0.4  km. 


Fig.  2.  Relation  of  Im(hr)  to  £7r2f;  notation  as  in  Fig.  1. 


Theae  ranges  are  such  that  the  practical  f' very  seldom  fall 
outside  them.  An  increase  in  f  displaces  the  range  in  £ 
towards  higher  values.  The  least  variation  in  ? (about  a 
faotor  of  100)  for  given  f  and  r  occurs  for  Im(h2);  the  other 
components  correspond  roughly  to  factors  of  1000.  On  the 
other  hand,  the  real  parts  of  h«  and  hr  are  large,  which  would 
make  it  very  difficult  to  isolate'  the  imaginary  parts  of  these.  I 
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f  Figures  1-i*  also  show  these  functions  of  f/r2 f  for  dipoles 

at  heights  of  $0  and  100  m,  positions  that  are  used  in  aerial 
surveying  by  induction;  these  curves  have  been  computed  from 
relationships  given  by  Waite  /13,1^7.  The  absolute  values  of 
the  Magnetic  numbers  decrease  as  the  height  increases,  while 
the  sene  of  maxima  shifts  to  smaller  ?/r2f.  This  means  that 
the  usable  ranges  are  extended  towards  lower  in  aerial 


Fig.  3.  Relation  of  Re(hr)  and  IhjJ  to  see  Fig.  1 

for  symbols. 


(The  subject  of  aerial  surveys  generally  requires  a 
study. ) 


surveys 
j _ special 
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Similar  relations  between  the  field  components  and  f  may 
be  derived  for  electromagnetic  dipole  profiling,  in  which  two 
mutually  perpendicular  dipoles  are  used  as  sources  in  conjunc¬ 
tion  with  two  identical  receiving  dipolea  (the  rotating-field 
method).  In  the  rotating  magnetic  field  method  /§/,  two 
identical  but  mutually  perpendicular  source  coils  are  used  in 
conjunction  with  two  detecting  coils  set  at  right  angles  to 
their  direction  of  displacement.  The  two  source  coils  are  fed 
with  currents  differing  in  phase  by  90® ,  which  produces  a  rotat 
ing  magnetic  field  in  any  homogeneous  insulating  medium.  The 


—  p 

Fig.  4.  Relation  of  Im(hr)  and  Im(hj>)  to  f/r  f;  see 
Fig.  1  for  symbols. 


signals  received  by  the  detecting  coils  also  differ  in  phase  by 
90°.  The  orientation  is  such  that  the  field  component  sensed 
by  the  vertical  magnetic  dipole  (a  horizontal  coil)  is  Im(H_), 
while  the  one  sensed  by  the  horimontal  dipole  (vertical  coil) 
is  Hf.  The  signal  in  the  vertical  dipole  is  displaced  in 
phase  Jby  90°,  so  the  sum  of  the  two  signals  is 


J 


38 


Fig.  ?.  Relation  of  J hr|  to  ^7r^f  for  dipoles  on  the 
surface  of  the  ground. 


resultant  signal  will  be  zero  for  a  homogeneous  insulating 
space,  but  a  conducting  half-space  upsets  this.  The  real  and 
imaginary  parts  of  the  unbalance  signal  are  dependent  on  f, 
and  r,  as  well  as  on  the  height  (for  aerial  surveys). 
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**  p  1 

Figures  5  end  6  show  these  as  functions  of  £/r  f  for  heights  of  i 
25,  50,  sad  100  a  on  the  basis  of  white's  results  /i 3,14/  in 
conjunction  with  results  for  the  surface  /3,**7*  Figure  5  shows 
that  the  real  part  has  a  maximum  at  ^/r2f  *4.5  for  the  surface, 

there  being  a  rapid  rise  to  this  point  from  the  zero  at 
l^r^f  ■  2.25*  The  real  part  is  almost  constant  at  about  2.0  for 
( yr^t  ■  0.15  or  so.  The  relationship  is  of  the  same  fora  for 
all  heights,  but  it  shifts  to  smaller  ?/r2t  as  the  height 
increases;  e.g.,  the  maximum  lies  at  0.17  when  the  height  is 
25  ■«  the  sero  being  at  0.052.  The  corresponding  values  for 
50  a  are  0.09  for  the  maximum  and  a  value  outside  the  limits  of 
the  figure.  Both  points  are  outside  the  limits  for  a  height 
of  100  a. 


The  imaginary  part  is  negative  at  the  surface  of  the  ground 
and  haa  its  minimum  at  [/r^t  =0.42  (Fig.  6);  the  minimum  lies 
outside  the  figure  for  a  height  of  25  m.  The  positive  branch 
appears  on  the  right  for  a  height  of  50  a;  this  has  a  maximum 
at  f/r^-t  ■  0.55,  the  zero  be  ins-  at  0.11  and  the  minimum  outside 
the  figure.  Only  the  positive  branch  (maximum  at  0.24)  is 
seen  for  100  a. 


These  curves  show  that  the  imaginary  part  is  of  the  main 
interest  for  ground-level  studies,  because  this  shows  a  single- 
valued  relation  to  distance  for  f/r^t  >0.42;  for  aerial 
studies,  the  real  part  is  more  useful,  for  this  gives  a  single¬ 
valued  relation  for  f/r^t  >0.17  (25  m),  >  0.1  (50  m),  or 
>0.07  (100  a). 


Induction  results  can  also  be  represented  in  apparent- 

resistance  form;  all  that  is  needed  is  a  certain  modification 
of  the  field  observations,  for  the  existing  apparatus  enables 
one  to  measure  the  amplitudes  of  the  components  as  well  as  the 
orientation  of  the  magnetic  vector.  These  amplitudes  do  not 
vary  rapidly  over  fairly  highly  differentiated  sections,  so 
the  method  would  be  of  comparatively  poor  resolving  power;  this 
is  a  result  of  the  distribution  of  the  normal  field.  The 
source  coil  acts  as  a  horizontal  magnetic  dipole  in  the  induction 
method,  and  the  field  components  are  examined  in  the  equatorial 
plane  (^»  90°).  Now  (5)  shows  that  the  corresponding  hr  and 
hg  are  sero  no  matter  what  f  may  be  if  there  are  no  vertical  or 
inclined  inhomogensities;  only  ht  is  finite,  and  this  here 
coincides  with  h*.  The  f  used  in  induction  commonly  range 
from  2.5  to  7.0,  and  the  modulus  of  ht  varies  little  in  this 
range,  so  the  resolution  in  %  is  low.  A  much  more  useful  case 
is  that  in  which  the  dipole  is  oriented  for  measurements  along 
its  axis  (p=  0°),  for  here  hf  is  zero,  whereas  hr  and  hz  are 


LT1-1 


,  if  the  medium  is  homogeneous.  Figure  5  shows  the 
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Fig.  8.  Curves  for  dot  training  r*  by  the  inclined  coil 

not hod. 


absolute  values  of  hr  and  hz  as  functions  of  ^Vr^f;  here  the 
real  and  imaginary  parts  increase  together  for  2  >  >0.2. 

The  frequencies  used  in  induction  (18.75,  37.5,  and  75  kc)  give 
hr  such  that  f  in  excess  of  100  ohm-a  can  be  measured,  but  the 
deteraination  is  not  unambiguous  for  lower  for  a  maximum 
occurs  in  hr. 

Measurements  of  /hjJ  can  be  used  to  calculate  "f  directly, 
for  this  increases  regularly  froa  1.966  for  J  *  1.4  to  4.076 
for  1  *  8.0;  Fig.  7  shows  this  for  hr  as  a  function  of 
A  range  of  10-fold  variation  in  (  is  accessible,  which  for  r 
of  0.1  ka  corresponds  to  40-400  ohm-rn  at  18.75  kc,  75"750  oha-a 
at  37*5  kc,  and  150-1500  oha-a  at  75  kc.  Further,  f  can  be 
found  froa  angle  measurements;  I  oait  details  and  merely  state 
that  existing  recommendations  require  operation  at  fixed  f  and 
variable  r;  the  results  are  used  to  determine  the  r  corres¬ 
ponding  to  ]  >  1.  Now  f  is  known,  so  can  be  found,  but  the 
.operation  is  very  tedious  and  is  seldom  used.  In  principle,  j 
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1  though,  9  »Bd  ?  can  bo  determined  from  tho  inclination  for  one  i 
r  only}  tho  uoo  of  several  r  only  impairs  tho  interpretation, 
for  f  ii  a  function  of  r  if  tho  medium  is  inhomogeneous.  Tho 
inclined-coil  gethod  is  simplified  if  r  is  kept  constant;  tho 
calculations  can  be  simplified  or  eliminated  if  f  is  replaced 
as  absoisaa  by  f/r2f,  or  even  by  Z  itself  for  the  fixed  r. 

Figure  8  shows  the  angle  of  inclination  du  a  function  of  f 
for  f  e  18.75  kc;  the  point  a  corresponds  to  $  *  1.  Then  the 

f  for  18.75  kc  is  read  directly  as  the  abscissa  of  the  point  on 
the  curve  for  the  given  a  whose  ordinate  is  tp.  The  Z  so  read 
off  must  be  doubled  at  37.5  he  and  quadrupled  at  75  kc.  Of 
course,  a,  the  inclination  of  the  source  coil,  can  be  varied 
during  the  measurements. 

i 

In  detailed  surveys,  in  which  ^must  be  determined  as  a 
function  of  depth,  it  ie  best  to  use  several  frequencies.  The 
above  curves  show  that  Z  is  found  the  more  precisely  the  larger 
the  difference  between  a  and  The  differences  are  large  only 
for  the  lower  values  of  Z*  though  the  above  frequencies  give 
reasonable  accuracy  for  r  of  50-100  m  for  Z  up  to  500  or  600  ohn-eu 

Wide  ranges  in  Z  can  be  covered  by  existing  standard  induc¬ 
tion  equipment  only  if  the  frequency  range  is  extended  downwards. 
The  precise  frequencies  required  should  be  derived  from  field 
trials* 


The  main  conclusions  to  be  drawn  are  as  follows. 

1.  A  function  of  apparent-resistance  type  should  be  used 
in  a.c.  methods;  the  value  and  dimensions  of  this  are  the  same 
as  those  of  the  specific  resistance  fur  a  homogeneous  medium. 
This  would  facilitate  the  use  of  electrical  methods  in  structure 
mapping  and  would  provide  a  fuller  picture  of  many  geological 
sections.  Special  subscripts  and  superscripts  are  needed  to 
Indicate  the  type  of  equipment  and  field  used. 

2.  The  apparent  resistance  is  calculable  from  the 
observations  only  if  the  reduced  distance  can  be  determined  unam¬ 
biguously;  the  conditions  that  ensure  this  are  dependent  on  the 
fields  used  and  on  the  nature  of  the  section. 

3.  The  value  of  Z can  *>•  determined  in  principle  for  any 
distance  from  the  field  of  a  vertical  or  horizontal  magnetic 
dipole;  the  methods  of  observation  can  be  greatly  simplified 


for  short  distances.  For  example,  only  aaplltudea  aay  ba 
naadad  for  soae  ooaponanta  (hr  and  hs,  for  inatanea).  Tha 
maximua  uaabla  dlatanea  in  that  eaaa  la  dapandant  on  tha  flald 
ooaponant  aaployad.  Tha  uaual  diataneaa  between  dipolaa  (50  to 
200  a)  and  resistivities  (1  to  100  oha-a  or  ao)  raquira  frequen- 
elaa  that  ahould  not  axcaad  7-8  ke.  Thaaa  fraquenclaa  can  ba 
lncraaaad  by  faetora  of  1.5  or  2  for  aarlal  uaa. 

4.  Induction  apparatus  employing  higher  fraquenclaa  can 
alao  ba  uaad  to  aeaaure  f ;  rocks  of  raaiatanca  a  few  hundred 
oha-m  can  ba  axaainad  by  tha  inclinad-coil  method,  but  tha 
technique  auat  ba  altered.  High  valuea  of  &  enable  one  to 
operate  merely  with  lhrl  or  lhBl  ,  or  sometimes  only  with  lh,l. 

Thia  diacuaaion  haa  bean  concerned  only  with  theoretical 
poaaibilitiaa  for  tha  uaa  of  apparent  raaiatanca.  Special 
atudiaa  are  naadad  to  evolve  detailed  methods  and  techniques. 


Bibliography 


1.  Van 'yen,  L.  L.  Nekotoryye  voprosy  chaatotnykh  zondirovaniy 
gorizontal'nykh  naplaatovaniy  [Soaa  aspects  of  frequency 
sounding  for  horizontal  strata].  Prikladnaya  Oeofizika 
[Applied  Geophysics],  Ho.  23.  Gostoptekhizdat ,  19^9. 

2.  Idea.  Elementy  teorii  stanovleniya  elektroaagnitnogo 
polya  [Elements  of  the  theory  of  build-up  of  an  electro¬ 
magnetic  field].  Ibid.,  No.  25,  I960. 

3.  GasanenkO,  L.  B.  Noraal*noye  pole  vertikal ' nogo  garmoni- 
cheskogo  nizkochastotnogo  mognitnogo  dipolya  [Normal  field 
of  a  vertical  harmonic  low-frequency  magnetic  dipole]. 

Och.  Zap.  LQU  [Archives  of  Leningrad  University],  No.  24 9, 
issue  10,  1958. 

4.  Gaaananko,  L.  E.  and  G.  V.  Molochnov.  Elektromagnitnoye 
pole  gorizontal'nogo  magnltnogo  dipolya  nad  gorizontal *no- 
sloistoy  strukturoy  [Electromagnetic  field  of  a  horizontal 
dipole  over  a  horizontally  layered  structure].  Ibid. 

5.  Mikhaylov,  I.  G.  Elektromagnitnoye  pole  naklonnoy  ramki 
[The  alactroaagnetic  field  of  an  inclined  coil].  Mater. 
TaHIORI.  Oeofizika  [Archives  of  the  Central  Geophysical 
Prospecting  Research  Inatitute,  Geophysics],  coll.  2, 

Moscow- Leningrad,  1936. 

L  J 


104 


Moloehnov,  G.  V.  Chaatotnoye  elektromagnitnoye  zondirovanly*] 
aa  aedali  [Frequency  electromagnetic  sounding  on  models]. 

Och.  Zap.  LOP.  No.  278,  laeua  11,  1959. 


7#  Tarkhot,  A.  G.  Gaoflsichaakaya  razradka  mstodom  induktsll 
rOaophysioal  prospacting  by  the  induction  mat hod].  Moscow, 
Gosgeolisdat,  1954. 


8.  Tarkhov,  A.  G.,  E.  M.  Ershov,  A.  I.  Katskov,  S.  R.  Maron, 
and  8.  F.  Fadorov.  Instruktivnyye  ukazaniya  po  provadaniyu 
geof lzicheaklkh  rabot  natodoa  induktsii  [Instructions  on 
gaophysical  work  by  tha  induction  method].  Moscow, 
Goagsolisdat ,  1959. 

9.  Shaub,  Tu.  B.  0  primenenii  aatoda  vrashchayushchego 
magnitnogo  polya  dlya  elektrokartirovaniya  [Use  of  tha 
rotating  magnetic  field  method  for  electrical  mapping]. 

1st.  AN  SSSR  [News  of  tha  Academy  of  Sciences  of  tha  USSR], 
Ho.  10,  I960. 


10.  Sheynman,  S.  M.  and  L.  V.  Labadkin.  Apparatura  dlya 
mnogochastotnoy  elektrorasvedki  paremennym  tokom.  Novoye 
v  aetodike  i  takhnika  geolograzvedochnykh  rabot  [Apparatus 
for  multifraquancy  a.c.  electrical  prospecting.  New 
developments  in  methods  snd  techniques  for  geological 
prospacting  work].  Trudy  VITR  [Transactions  of  the  All- 
Union  Prospecting  Technique  Institute],  coll.  2,  Gostop- 
takhisdat,  1959. 

11.  Lemberg,  A.  E.  A  new  electromagnetic  device  for  the 
surveying  of  subsurface  conductors  and  its  application  in 
Canada.  Can.  Min.  Met.  Bull.,  No.  485,  Sept.  1952. 

12.  Heltner,  L.  and  P.  Lasfar^ues.  Le  developement  actuel  da 
la  glophysique  at  de  la  geochimie.  Application  a  la 
recherche  des  gisements  metallif ires.  Rev.  Ind.  Min.,  40. 

Ho.  3»  1958.  ~  ~ 

13.  Wait,  I.  R.  Mutual  electromagnetic  coupling  of  loops  over 
a  homogeneous  ground.  Part  I.  Geophysics,  20.  No.  3,  1955* 

14.  Whit,  I.  R.  Mutual  electromagnetic  coupling  of  loops  over 

a  homogeneous  ground.  Part  II.  Geophysics.  21 .  No.  2,  1956. 


J 


105 


r 


i 


THE  USE  Or  LOW- FREQUENCY  ALTERNATING  CURRENT  IN  ELECTRICAL 
PROFILINO  AND  SOUNDING 

Yu.  I.  Bulgakov  and  A.  V.  Vaahav 

pp.  187-192 

Industrial  intarfaranca  causae  serious  difficulty  in 
alaetrical  prospecting;  the  intarfaranca  is  especially  strong 
around  ora  bodies  under  exploitation  and  in  mining  areas  in  the 
Ural,  Rudnyy  Altai,  and  so  on.  It  is  then  usually  iapossible 
to  use  d.c.  methods.  This  interference  has  many  different 
sources  and  shows  a  vary  complex  time  course,  which  we  have 
studied  in  detail  at  the  Belousovo  and  Grakhovo  deposits  in  the 
Rudnyy  Altai.  The  low-frequency  intarfaranca  was  recorded 
with  an  EAK-3  electronic  potentiometer  /V*  noise  of  frequency 
$-8  c/s  and  above  was  recorded  with  a  small  portable  oscillo¬ 
graph  type  MPO-2  fitted  with  a  wide-band  amplifier,  which  pro¬ 
vided  coverage  up  to  300  c/s.  The  detecting  line  was  20  a 
long,  the  largest  potentials  being  30-40  mv .  Figure  1  shows 
a  typical  curve,  the  tine  marks  being  of  frequency  300  c/s 


is  the  width  of  the  peak. 


(the  abscissa  is  graduated  at  intervals  of  0.003  sec).  The 
frequency  spectrum  is  very  wide. 

Figure  2  gives  frequency  analyses,  which  were  made  by  a 
graphical  method  designed  for  seismograms  fli/.  Curve  a  is 
a  general  one,  while  the  others  represent  the  behavior  for 
limiting  frequencies  of  36,  100,  and  373  c/a.  The  intensity 
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Fig .  2.  Frequency  analysts  of  interfersncs;  h(f) 
is  tbs  aaplitude  at  frsquaney  f,  ths  unit  being 
0.01 A*,  in  which  A^  is  the  aaxiaua  amplitude. 


is  lowest  at  17.23  e/s,  with  other  ainiaa  at  37  and  73  c/s. 
Ths  intensity  is  highest  at  0-3  c/s  and  at  30,  100,  and 
150  e/s.  The  17-23  c/s  region  is  clearly  the  best  for 


|  electrical  prospecting, 


for  the  frequency  is  reasonably  low 


10? 


and  on*  may  axpect  that  tha  results  for  the  short  lines  usually  J 
adequate  for  ore-bodies  will  be  virtually  the  ease  as  those  for 
d.c.,  so  the  methods  for  d.e.  can  be  used  in  interpreting  the 
results. 

Theoretioal  relationships  for  a  homogeneous  half-space  and 
for  two-  and  three-layer  media  are  applicable  (if  the  upper 
layers  are  of  very  high  resistance)  in  order  to  compare  the 
fields  and  f jj  for  d.c.  and  a.c.  (dipole)  methods.  Stefanescu 
/6J  has  discussed  the  electromagnet  field  of  a  grounded 
horisontal  electric  dipole  on  the  surface  of  a  homogeneous 
half-space;  his  tables  show  that  the  amplitude  of  the  radial 
component  of  the  electric  field  is  as  for  d.c.  if  ^  <  1.0; 
the  same  applies  to  the  asimuthal  component  to  5K  if  5  *  0.7 • 

Here  $  ■  xr,  in  which  r  is  the  distance  from  the  center  of  the 
dipole  (km)  and  X  *  0.4?(5 f/f)#,  (  being  the  resistivity  in 
ohm-m  and  f  the  frequency  in  c/s.  The  above  limits  for  ^ 
imply  that  a  dipole  00'  operating  at  20  c/s  over  rocks  of 
specific  reaistaneee  20  and  200  ohm-m  has  maximum  working 
ranges  of  350  and  1100  m  respectively  (radial),  or  2 50  and 
800  m  (azimuthal). 

Similar  estimates  can  be  made  for  a  horizontally  layered 
medium  (also  for  20  c/s);  Van'yan's  three-layer  frequency 
sounding  curves  can  be  used.  An  unfavorable  section  is 

one  having  hi  ■  10  e,  (i  «  200  ohm-m,  fc2  *  50  m,  ■  20  ohm-m, 
fx  « OB,  and  Si  ,2/^1  *  51.  The  wavelength  for  the  first  layer 
is  Ai  ■  3.162  x  105(^/f)#  ■  10®  cm  •  10  km,  and  that  for  the 
second  is  X2  «  3.162  x  105(f/f)H  .  3.162  x  10®  cm  *  3.162  km, 
so  Ai/hi  m  1000  and  A2/h2  ■  63.2.  The  00’  for  the  determin¬ 
ation  of  hi  +  h2  may  be  as  much  as  1000  m  [00'/(hi  ♦  h2)  «  16.73. 
The  above  relations  enable  us  to  use  S  curves.  The  largest 
differences  in  the  fk  will  occur  for  the  maximal  distances,  so 

we  perform  our  estimate  for  this  case. _ The  S^A  curves  give  us — 

that  ffc.20  i»  for  the  axial  dipole  system  and  is  1*05fk 

for  the  equatorial  setting.  No  curves  are  available  for  the 
case  of  /,  finite,  but  we  give  below  some  examples  from  field 
work  of  1957  at  20  c/s  for  ^3  large  but  finite. 

We  used  an  EP-1  electronic  potentiometer  and  a  direct- 
reading  ESK-1  compensator  in  the  d.c.  measurements;  a  vacuum- 
tube  millivoltmeter  with  a  selective  amplifier  was  used  in  the 
a.c.  measurements.  The  millivoltmeter  has  already  been 
described  /V .  Results  for  a  profile  at  Kamyaha  (AO  *  OB  •  100  a 

00'  ■  100  m,  AB  *  MN  «  10  a,  the  last  being  the  pitch)  gave  almost 
identical  fk  curves,  the  largest  discrepancies  being  about  5X 
and  the  mean  difference  2.6%,  which  is  within  the  error  of  the 
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Fig.  3>  Curves  for  ^for  symmetrical  and  dipole  methods 
performed  with  d.c.  and  a.c.  in  the  Kamysha  area  in  the 
Rudnyy  Altai]  1)  oyaaetrical  electrical  aethod,  d.c. 
(profile  1,  PK-29 ) !  2)  the  same  with  a.c.  (MM  »  2  a); 

3)  (k  b7  d.c.  log  in  borehole  47;  4)  the  eaae  by  an 

a.c.  aethod. 


observations.  Figure  3  shows  ^curves  from  syaaetrical  and 
dipole  aethods  of  vertioal  electrical  sounding;  AB  *  800  a 
(aaxiaua)  for  the  first  and  00'  ■  ^50  n  maximum  for  the  second 
(axial).  Bare  the  curves  differ  by  only  1.2$.  In  addition, 
curve  2  is  for  AB  *  800  a  with  MM  *  2  a;  this  short  distance 
was  possible  on  account  of  the  high  sensitivity  of  the  vacuua- 
tubs  voltaeter  (30  uv  for  full-scale  deflection  on  the  M-24 
aster),  so  reliable  readings  could  be  made  at  10-13  uv,  which 
is  oae*twentieth  of  the  ainiaua  potential  that  can  be  measured 
with  adequate  precision  in  d.c.  methods. 

There  were  no  parts  to  be  matched  up  (on  account  of  the 
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use  of  differing  values  of  MN),  so  interpretation  of  the  fa 
facilitated.  For  exanple,  there  is  a  22%  reduction  in  toe 
nuaber  of  aeasureaents  for  the  curves  of  Fig.  3  (14  instead  of 
18)  if  only  one  MN  is  used. 

These  results  show  that  the  line  lengths  usual  in  ore- 
bearing  regions  (about  1  ka)  give  fa  at  20  c/s  virtually  the 
saae  as  those  froa  d.c.  methods,  so  relationships  established 
for  d.c.  can  be  used  in  the  interpretations. 

Measureaents  have  also  been  aade  in  other  regions  of 
strong  interference;  in  every  case  good  and  reproducible  results 
were  obtained  at  20  c/s. 

The  advantages  of  low-frequency  a.e.  are  as  follows! 

1.  The  work  is  aore  rapid  and  detailed  while  retaining 
the  saae  accuracy. 

2.  Measureaents  can  be  aade  in  the  presence  of  strong 
industrial  interference. 

3.  The  apparatus  can  be  light  and  portable,  for  it 
deaands  less  power  in  the  measuring  circuits,  so  the  power 
sources  are  lighter. 

It  would  be  advantageous  in  aany  cases  to  transfer  froa 
d.c.  aethods  to  a.c.  ones. 
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KBTHOOi  OP  CALCULATING  A  POINT- SOUfiCE  FIELD  IN  THE  PUESEWCE 

OF  A  VERTICAL  BED 
N.  A.  Sakharnikov  and  D.  M.  Volkov 

pp.  193-202 

It  is  necessary  to  solve  for  the  potential  distribution 
for  the  field  of  a  point  souroe  in  the  presence  of  two  interfaces 
(Fig.  1)  in  order  to  calculate  apparent-resistance  curves  from 
eleotrioal  profiling  above  a  stratum  of  low  or  high  resistance? 
The  potential  at  any  point  in  the  field,  U(r,Z),  can  [\]  be 
represented  ae  follows  for  media  I,  II,  and  III  respectively! 


«,<r.  (<•*)«-“ X 

«.v.  0-'£0 + *,,)j  AW-" 1 
«,<<•.  O-fe (!+*,.)(!  +  *«) 


in  which 


fc  —-fr-TA  .  t  — lift. 

**  H  +  Pi  ’  **  P»  +  H  * 


(1) 

(2) 


*A.  V.  Veshev  posed  the  problem  of  the  best  method  of 
[calculation. 


_1 
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FTere  Z  is  the 
1  resistivities 


IN  always  lass 


ourrant  from  the  souroa  electrode)  tha 
and  pj  oan  vary  widely,  while  k 


ftjk’ 


one.  The  potential  gradient  is 


* found 


by  differentiating  the  expreaeions  in  (l)  with  respect  to  r  or  I. 


Fig.  1.  Coordinate  system  for  oaloulating  the 
field  of  a  point  souroa  in  tha  prasanoa  of  two 
interfaces.  A,  position  of  souroa  electrode) 
lit  position  of  deteeting  eleotrode;  1  ,  distance 
between  aleotrodas)  d,  distance  from  A  to  nearer 
interfaoe)  r,  distanoe  from  Li  to  surface  of  ground) 
f,  resistivity  of  corresponding  medium* 


For  the  surface,  r  *  0  and  J  (rx)  •  1,  so  the  formulas  for 
U(r,l)  and  B(r,l)  are  simplified)  the  apparent  resistances  for 
potential  and  gradient  systems  respectively  are 

p—  2*1*  -j~. 

It  is  usual  to  express  the  apparent  resistanoe  by 
formulas  [\J  oontaining  infinite  series,  which  for  media  I  to 
III  ta*.e  tue  following  form  for  potential  systems 


FAi3**!  [l  +  HnT/  +  ^  Id  +*] (a*  *«)]  » 

F/H^FlO  +  *«)  +  hi'+iA  (a»  *|)]» 

=  h  0  +  ^i»)0  +  bt); 


L 


_l 


r. 


and  for  gradient  systems 


[t  -  s*<a>  *»>  • 

P»t  -  Pi  0  +  *»)  “ (M+tkZfy  S*(a>  *»)  « 

P*i ==  Pi (1  +  *»•)(!  +  ^n),,S(fl»  ^»)» 


in  whioh 


(4*) 


s‘<‘-  ^-SlfTT&r  <1-1,  <-2). 

•nd  ^  2* 

d"^n‘ ■  ^JT+TT— T 1  ^  “T  ‘ 

Series  (5)  converges  for  (a|  <  1,  b  *  0  and  for  a  *  1, 
b  >  0  for  i  ■  1  and  i  ■  2,  as  well  as  for  a  «  -1,  b  >  0  for 
i  *  2.  The  ranges  of  praotioal  intersst  are  as  follows i 

0.6<|*„|<O,96.  0,0B<-J-<8, 

0.6  <  |  *„  |  <0,96,  0<  >  <  10. 

It  is  shown  below  that  replacement  of  the  infinite 
series  in  (4)  by  particular  sums  results  in  a  large  volume  of 
oomputation  if  the  thickness  h  is  small  and  18  01088  to 

el  or  -1.  Direct  a u.nmat ion  for  p.  involves  tne^caloulation  of 
the  sum  in  (5)  from  the  formula 


(5) 

(6) 


#-» 

5|<«.  *)=*=  ‘ 

’■*  (7) 

Now  we  estimate  the  number  N  of  terms  in  (5)  that  must 
be  used  in  direct  summation.  Let  A  be  the  absolute  error 
allowed  in  the  approximate  value  Of  S. (a,b)j  then  a  theorem 
on  sign-varying  series  fcj  gives  us  that  the  N  for  a  >  0  must 
satisfy 

(I  A» 

L  ,  _1 
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*hiehf  for  b  small, is  equivalent  to 


If  •  <  0,  the  terms  in  (5)  are  positive, 
the  V  terms  is 


(8) 


so  the  remainder  after 


R< 


yN 


0  +  wv)'(i-l«l>' 


To  make  H  we  must  make  N  such  that 

which  for  b  small  is  equivalent  to 

log  1  a  |  ~  tog  I  a  | 


(9) 


L 


Pig.  2.  Graphs  of  N  *  f(a)  for  various  A,  for  use 
with  (7)}  aw  ^.-k  ?  ie  *he  *east  ®  that  will 

satisfy  (8)  or  (97,  -^is0the  absolute  error  allowed  in 
S.(a,b)}  values  of  &<  l)  0.00001}  2)  0,0001} 

37  0.001. 
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Figure  2  shows  the  relation  of  a  to  No,  the  least  N  that 
satisfies  (8)  or  (9),  for  various  6.  Thesegraphs  enable  one  to 
determine  N  for  the  sum  of  (5)  from  (7)  of  0  <  b  <  0.1$  the 
curves  for  I  as  a  funotion  of  a  run  somewhat  below  those  of  Fig.  2 
if  b  >  0.1,  but  the  differenoe  is  unimportant,  for  N  varies 
little  for  Ja)  small  and  deoreases  as  |a|  beoomes  large,  but  not 
so  muoh  so  that  (7)  beoomes  useful. 

It  is  practicable  to  use  (7)  for  |a|  small,  i.e.,  for 
-0.6  4  k12k23  *  0.6j  th*n  *  ^20. 

Some  ways  of  simplifying  the  oaloulations  are  illustrated 
below.  Consider,  for  example,  the  integral 


l,c.  <<=1. ' 1=2,1 

which  converges  for  |a|  4  1,  b  >  Oj  is  suoh  that 

if  l  >  0,  as  oan  be  seen  if  we  replaoe  x  by  x/l. 


(10) 


(ID 


te  oan  represent  f.  in  the  following  integral  form,  as 
(l)-(3)  imply  for  r  •  0,  if  we  use  (ll)  in  conjunction  with  the 
definition  of  L^ayb).  For  potential  systems 


-  L  ,  .  (‘ -*?»)*»'  .  *  J 

+  jjttH — W7%r=rrL'(a’b')\' 

=  *i)J,  (12) 

Pm  =  P»(*  ■+■*!*)  (I 


J 
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|  and  for  gradient  ones 

f,  *,.<• 

P*|—  P«Ll  (2d- t)'  ~(2d  I  2A-/)* 

Pm  —  fcO  ^i)|> 

For  example,  from  the  first  formula  in  (12), 

ee 

P,i  =  2*/-^  =p,/-j70(r.r)<r,^l  -j-*Ite-2W-0jr 


For  r  =  0 


p»i — pi  *  [7  +  afer  +  amcrr/  (g>  *>)]■ 


(12') 


Formulas  (4)  can  be  derived  from  (12)  by  replacing  L.(a,b)  by 
S^(a,b){  standard  transforms  f\J  enable  us  to  show  that  these 
functions  are  identically  equal,  for  which  purpose  we  replace 
the  integrand  in  (10)  by  a  convergent  series 

£,(«,  *)=J[^(- a)«.«-"+,*>* jrfx. 

The  sequence  of  integration  and  summation  can  be  reversed  in 
this  case  [2],  so 

m  «c 

1,(0,  *)=S(-c)"  f«-|U*“,'£<x=S1(a,  *). 

*-0  6 

In  the  same  way  we  can  show  that 


dx  = 


for  all  X  >  0}  differentiating  k  times  with  respect  to  I 
(which  is  permissible  /£/),  we  have 

L  -  _l 


(jZl L 

l  +  bn 


M-0 


P  t  ,X 

J  1  +  at-*1 
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f 


x**'u  r  v,  V  <-•>* 

—rr  ax  — k\  y ,  ---  —X- 


—  - 


whioh  applies  for  all  l  >  0,  b  >  0,  and  )a|  *  1.  In  particular, 
if  k  9  1  and  1  *  1,  wa  have  from  (5)  and  (10)  that 

i|(®i  — 5j(®»  ^)« 

The  integral  of  (10)  oan  ba  rapraaantad  exaotly  as  a 
finita  sum  of  elementary  functions  for  all  a  and  for  rational  b 
if  i  *  1,  as  may  ba  shown  by  replacing  tha  variable  of  integration 
in  (10)  by  means  of  the  substitution  exp(-x)  s  yi 


*)  = 


J 


*y 

i  +«/  ’ 


The  oorrasponding  indefinite  integral  oan  [l]  ba  represented  in 
finita  form  for  any  rational  b  (the  substitution  is  t'  s  y,  in 
whioh  b  *  p/e),  so  L^(a,b)  is  an  elementary  function  of  a  for 
each  rational  b.  In  partioular 


Li(a,  l)=~ln(I-f  -«), 
Lt  (a,  2)  =  -~r  arctg  Va%  if 


Lx  (a,  2)  =  - 


1 

2 /Til 


In 


1-/TTT 

i  +  /T*\  ' 


a>  0 

if  a  <  0. 


Let  b  ■  l/n,  in  whioh  n  is  a  positive  integer}  n  >  1.  Then 
for  all  a  we  have  exaotly  that 

L*  v** 

V  V  (14) 

Tha  following  transforms  of  the  sum  of  (5)  serve  to 
demonstrate  (14)> 


_1 
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Proa  this  (14)  follows  dlrsotly  If  we  use  the  expansion  of  the 
logarithm  function  as  a  power  series. 


(14)  is  suitable  for  any  a  and  b  >  0.05,  and  also  for 
smaller  values  of  b  if  )a)  is  sufficiently  small.  The  series  on 
the  right  in  (14)  converges  more  rapidly  than  the  series  in  (5)» 
so  it  is  simpler  to  use  (14)  for  pk  than  it  is  to  use  (7)* 

The  integral  in  (10),  or  the  equal  sum  of  (5)»  for  a 
gradient  system  can  be  operated  on  with  the  transforms 


«r 


n *  c*4* 

«*■  /!(«  :  *)’ 


for  i  *  2,  b  «  l/n  <  1,  and  a  <  0  to  give 


fc-J 

in  which  c  =  -a  and 

*-1 


(15) 


-  (16) 

Series  (16)  converges  for  all  |c|  <  1}  f(c)  satisfies  the 
functional  equation  (2 ] 


/(*)+/(! -c)  =  £  if  0<e<\, 


whioh  enables  one  to  calculate  f(c)  for  c  close  to  1  as  simnly 
as  for  c  close  to  zero,  for 


lnc  ln(l  -c)-f(\  — c). 
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|  direct  summation  of  (16)  oan  oa  used  for  o  oloso  to  aero)  (16)  | 

oonverges  aora  rapidly  than  (5),  so  p.  oan  ba  oaloulatad  froa 
(1$)  and  (16)  for  a  gradiant  system  aora  aimply  than  froa  (?)• 
Formula  (15)  applias  alao  for  c  *  -a  <  0,  i.a.,  if  k.gk_,  >  0. 

For  b  e  1  and  any  a  it  oan  ba  ahown  that  5 

Formula  (15)  oan  ba  uaad  for  any  a  and  b|  tha  worst.  _ 
value  of  a  aa  regards  tha  oooputation  la  -0.9|  but  than  a ^/25 
:  0.00016,  so  tha  aaoond  tarm  in  (15)  will  oontain  not  mora  than 
25  oomponanta  no  mat tar  what  b  may  ba. 


Tha  intagral  of  (10)  oan  ba  oaloulatad  approximately  to 
any  apaoifiad  aoouraoy  for  all  a  and  b  by  means  of  quadrature 
formulas  £}>] 1 


Lt 

Lt 


•  ft 

(a,  b)  —  j  e~*-f(x)  dx—^Ak  <p(xk) , 

(«.  b)=jxe~M-  f(x)dx—J^A*hf{x'k), 


(18) 


(19) 


in  whioh 

,w“TT3=*r- 

(20) 

The  n  of  (16)  and  (19)  is  Ohoaen  to  suit  tha  accuracy  needed 
in  p.j  n  5  usually  sufficient.  A  table  [$]  gives  tha 
x.  and  A.  of  (18)  for  various  nj  tha  x  and  x*  in  (l8)  and  ( 19) 
are  tha  roots  of  polynomials  of  degree  ns 


£?V)=<-  #•*--*£{*"**) 

for  dl  s  0  and  dL  >  1  raspactivaly.  Tha  A^  and  a£  are  given  by 

m  a!  r  («  t  a  -J-  1) 

H-WW 

also  for  si  •  0  and  d  *  1.  Table  1  gives  x£  and  a£ 
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Table  1 


Values  of  x*  and  A*  for  Use  with  Formula  (19) 


ft 

X* 

4 • 

I 

Xi  -  2,000000 

Ai  -  1.000000 

2 

x,  — 1,367  049 
x, -4,732061 

4,-0.788676 

4,-0.211324 

3 

x,-  0,036  822 
x,  -  3306  407 
x,  —  7,783  771 

4,  —  0.588681 

4,  -  0.391  216 

4j  —  0,(1)201  024 

4 

x  =  0.743291 
x,=  2371635 
x,-  5,731179 
xi- 10.963994 

4,-0.446869 

4*  =  0,477  635 

4,-0.11)741777 

4,-0.(2)131585 

5 

x,-  0.617011 
x,-  2,112966 
x\  =  44110833 
x,-  8.399071 
x  .  --  14,260099 

4,-0348014 

4,  =  0302  280 
4,-0.140916 

0.(2)871 983 
ft,  ■»  0.  (4)  889  738 

Note.  The  numbers  in  parentheses  indicate  the  number  of 
noughts  between  the  decimal  point  and  the  first 
significant  figure. 


The  integrals  of  (10)  can  be  calculated  for  a  close  to  -1 
means  of  slightly  altered  formulas;  we  put  f  ■*  1  ♦  a  to  get 


•o 

"C 

i‘Lt (a,  b)  —  jxe~Mf9(x)dx. 


(18') 
<19')  _j 


r 


Formula  (18)  can  be  used  to  find  the  integral  of  (id1),  and 
•ioilarlyt  (19)  oan  be  U3ed  for  (19* )»  here^(x)  is  replaced 
by  ^f(x),  which  is 


?o(*)  = 


I  r(S-l)s 


-»* 


and  which  decreases  monotonically  from  1  to  £  as  x  goes  from  0 
to  eo.  The  values  of  f£(x)  should  be  oalculated  to  a  high  ord*r 
of  aeouraoy,  for  the  integrals  in  ( 18 • )  and  (19*)  are  divided  by 
the  small  quantity  £.  For  example,  f  (x)  must  be  calculated  to 
a  further  significant  figure  for  f  r  071  (relative  to  (  •  0.9, 
that  is). 

Formulas  (18)  and  (19)  do  not  beoome  more  difficult  to 
apply  to  (10)  when  h  becomes  small,  although  (3)  converges  less 
rapidly  (for  a  fixed  a)  and  although  the  N  required  in  (7) 
increases.  In  fact,  even  when  k.-k-  =  -0.9  (the  worst  case), 
the  integral  of  (10)  converges  so  rapidly  that  its  value  differs 
by  less  that  0.01$  from  the  definite  integral  having  limits  of 
0  and  13* 

quadrature  formulas  can  also  be  used  for  r  /t  0;  ve 
merely  have  to  replaoe  x  by  x/l  in  (l),  which  gives 


in  which 


as 

U\(r,  j*  (x)dx, 

0 

so 

U,(r,  +  *„)  f  e--f,(x)dx. 

SC 

u*(r,  l)  -  -f  (1  +  *■*)(!  +  *,1)  je-'-KWdx, 


-m-n  T  -2(st*_/)4 

*1*  '  r-M  ‘ 


1  -f  *1 


.2lk-i 
/ 


L 


».  <■*)  =  /.  (y) :  (l  +  *„*»«" “ ' *) ,  q  =  . 


(20) 


_J 
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Eaoh  ^  (x)  (i  «  1*  2,  3)  decreases  monotonically  in  absolute 
value,  so  the  integrals  of  (20)  converge  rapidly,  for  exp(-x) 
is  present  as  well.  It  is  usually  sufficient  to  put  n  =  5  in 
( 18 )  in  order  to  calculate  the  potential. 
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INFLUENCE  OF  THE  EARTH-AIR  BOUNDARY  ON  THE  ELECTRICAL 


i 


FIELD  OF  A  POLARIZED  SPHERE 
pp.  222-225 


Q.  A.  Tarasov 


Buraian  ftj  and  Patrovskiy  ft,  3/  have  eonsidtrad  tha  affects 
of  tha  air-ground  boundary  on  tha  electric  fialds  of  polarised 
bodies,  but  they  have  not  aade  astiaataa  of  tha  aagnitude  of  tha 
affaet.  Tha  present  discussion  relates  to  tha  field  of  a  vert¬ 
ically  polarised  and  ideally 
conducting  sphere  of  radius  a 
at  a  depth  s« f rom  the  surface 
(Fig.  1).  The  polarisation  is 
such  that  the  external  field  in 
in  an  isotropic  and  unbounded 
aedium  may  be  replaced  by  tha 
field  of  two  charges  equal  in 
aagnitude  but  opposite  in  sign 
syaaetrieally  placed  with  respect 
to  the  center  of  the  sphere  a 
distance  2l  apart,  this  distance 
being  auch  less  than  a,  which 
aeans  that  the  polarisation  of  the 
sphere  is  as  nearly  uni  fora  as 
possible.  The  axis  of  polaris¬ 
ation  is  directed  vertically 
downwards.  The  field  and  the 
polarisation  are  affected  by  the  presence  of  the  boundary;  the 
potenetial  in  the  space  external  to  the  sphere  can  be  represented 
as  the  sua  of  two  eoaponents.  The  first,  U«,  includes  the  potent* 
ials  from  tha  two  real  and  two  iaage  charges;  the  value  of  !!• 
at  the  interface  is  twice  that  for  an  unbounded  aediua.  The 
second  includes  the  potentials  fiUn  from  a  system  of  fictitious 
charges  resulting  froa  successive  reflection  of  the  lntial  ones 
in  the  sphere.  The  potential  is  put  as 


Fig.  1.  Scheae  for  the 
electric  field  of  a 
polarised  sphere. 


(l) 

J 
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i  i  i  \ 

-  2<7o|  (js,  +  (H  _  [jt*  +  (*«  +  Z)*!175”  J 


or,  after  transformation,  as 


//  _ 

U*~~  a(P+1?W‘ 


(2) 


(3) 


in  which  X  •  x/a,  J  ■  s«/a,  and  X  ■  1/a..  The  quantity  in 
Xy(J*  +  f  * )  is  small,  so  it  is  neglected.  The  fields  of  the 
real  and  reflected  spheres  are  accounted  for  by  means  of  suc¬ 
cessive  reflections  of  the  charges  in  the  spherical  surfaces; 
the  correction  to  Dj  from  the  first  reflection  may  £\J  be  put  as 

iw*  l 


X 


[** 


['•“-E&T 


\  l  J  (4) 

rri75-  +  [-g37+  5^+7  | '  (*•  +  x*)'* 


7) 


or,  after  simple  but  tedious  transforms,  in  which  the  second- 
order  terms  are  discarded,  as 

*Ul==  sc(?tcy*  ‘  (p+c*-ij  • 

Then  the  potential  of  the  sphere  as  corrected  for  the  first 
reflection  is 


(5) 


£/,— - ^ 

*  -  I  P+\ 


Then 


**&  r»  » 


At/l  _ _ 

1 

»  \ 

p+c*  ; 

(6) 


(7) 


How  (7)  implies  that  the  quantity  on  the  left  is  maximal  when 


if/, 


! 


I 


U*  lt-o_4C^l— 


4C(P-I) 


and  minimal  when  S-^00! 


if/, 

f/e 


oc 


:4^s  • 


(8) 


(9) 


That  is,  the  maximum  AUi/U#  occurs  at  the  epicenter,  so 
the  fields  of  the  real  and  reflected  spheres  are  of  the  greatest 

L  interest  as  regards  this  point. 

Graphical  solution  of  (8)  shows  that  the  correction  to  Up j 
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for  the  first  rtf lection  is  lsss  than  j 
5^  if  7>  1.9  and  is  lsss  than  1#  if  ' 
5*3  (Fig.  2). 

By  analogy  with  U»f  tha  eorraction 
arising  froa  tha  second  reflection,  AOj, 
is  given  for  tha  turning  point  as 

—  4C4i(2C«-0(at*-3)  •  (10) 


v  *  — 


1 


(11) 


Than  tha  potantial  of  tha  sphara  as 
corraetad  for  tha  saeond  raflaetion  at 
tha  turning  point  is 

I 

4C<C“-I> 

+*»«»- w-*]- 

This  aaans  that 
I  At/il  2«»-l) 

Ixl  "tw- Dw-awtt*-  i)-n  •  <“> 

Graphical  solution  of  (11)  shows  t)iat 
tha  eorraction  to  B|  froa  tha  saeond 
raflaetion  is  lass  than  5)f  if  J  >  1.5 
and  is  lass  than  1)1  if  £*1.8  (saa 
Fig.  2).  Any  furthar  calculation  is  of  no  value,  in  view  of  tha 
saallnass  of  tha  quantitias. 

Than  wa  have  that  tha  alactric  field  is  given  within  5%  by 
tha  seroth  approximation  (simple  doubling)  if  tha  depth  of  tha 
eenter  of  tha  sphara  is  mors  than  twice  the  radius;  these  values 
are  usable  in  calculations  to  1&  if  tha  depth  is  more  than  3  times 
tha  radius.  A  eorraction  should  be  applied  for  the  first  reflect¬ 
ion  if  the  depth  is  lass;  this  correction  is  adequate  for  calcul¬ 
ations  in  which  up  to  5%  error  is  allowed  if  the  depth  is  more 
than  1.3  times  tha  radius,  or  1.8  times  the  radius  if  only  1% 
error  is  permissible.  Correction  should  be  made  for  the  second 
reflection  for  smaller  depths. 


Fig.  2,  Correction 
curves  (.%)  for  the 
extrema  value  of 
tha  potential  of  a 
polarised  sphere  on 
account  of  an  inter¬ 
face. 
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LOW-FREQUENCY  NCHMU  EQUIPMENT  FOR  PROSPECTING 
pp.  241-244  D.  R.  Barkhatsv 

The  following  is  a  brief  description  of  apparatus  designed 
at  the  All-Union  Institute  of  Prospecting  Geophysics  for  elee trical 
prospecting  at  low  frequencies;  at  present  it  is  used  mainly  for 
mean-gradient  measurements,  where  it  gives  very  good  reproducibility 
and  highly  detailed  results.  It  can  also  be  used  in  profiling  and 
sounding  with  dipole  systems  with  short  distances,  and  also  in  some 
other  instances.  The  apparatus  measures  potential  differences, 
which  can  be  reduced  to  apparent  specific  resistances  if  necessary. 
Thecomplste  set  of  equipment  consists  of  a  generator  combined  with 
a  converter  and  a  vacuum-tube  millivoltmeter. 


Fig.  1.  Circuits  of  the  generator:  Rt  to  R*  ,  160  ohms;  R»  , 
470/1;  14,  100  kit;  R,*,~10  kit;  R$,~2.5  kll;  R,*,~80it;  res¬ 


istors  R|  to  Rj  are  of  type  MLT-0.5;  R?  to  R,  are  adjusted 
during  calibration;  14  0  ,  4 -ft  (wire);  R,  t  ,  }0  Si  ( variable 
wire);  Rt*-R»,  30-A;  Rrt-Rtj,  300  Jl;  R*>,  1  +  0.01  fl  (wire); 
to  C*  ,  0.05  pF;  C*,  5  pF;  B,,  VS-4.8  vibrator. 
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[  The  main  part  of  the  converter  is  the  VS-4.8  vibrator  (Fig.  j 
1),  which  produces  an  alternating  current  at  100  +  10  c/s  with  a 
maximum  voltage  of  300  v  and  a  maximum  current  of  300  ma  (power 
limit  70  va).  The  current  in  the  supply  circuit  is  controlled  by 
resistors  R*  to  Rn,  whioh  are  inserted  in  series  with  the  supply 
AB  by  the  switches  TT}  and  H*  .  The  current  is  measured  by  the  M-24 
meter,  which  is  connected  via  a  bridge  rectifier  (silicon  diodes 
type  D-205),  The  same  instrument  serves  to  check  the  exciting 
voltage,  which  is  kept  between  4.3  and  6  v.  The  current  in  the  AB 
circuit  may  also  be  measured  with  the  millivoltmeter,  which  for 
this  purpose  is  connected  across  R*  (lit).  Spark-quenching  RC 
circuits  are  fitted  in  order  to  stabilize  the  operation  of  the 
synchronous  vibrator  (R»  to  %  and  C*  to  CjX  The  exciting  battery 
consists  of  four  1.5-TMTs-29.5  cells  (3s-L-30)  or  four  1.48-PMTs-9 
cells  (2s-L-9),  which  are  housed  in  the  generator  case.  The  sources 
of  current  are  batteries  such  as  29-GRMTs-13  (B-30),  69-PRMT8-6 
(B-72),  102-AMTs-U-l. 0  (BAS-80),  and  so  on,  which  are  connected 
externally.  The  converter  is  contained  in  a  metal  ease  170  x  200 
x  310  mm;  the  weight  without  batteries  is  2  kg. 

The  vacuum-tube  voltmeter  (Fig.  2)  consists  of  the  input 
divider  I,  the  amplifier  II,  the  calibrator  III,, and  the  measuring 
instrument  IV.  The  input  divider  (range  switch)  provides  ranges 
of  1,  3,  10,  30,  100,  300,  and  1000  mv  full  scale,  the  input  res¬ 
istance  being  constant  at  2.4  M  .  The  resistance  of  the  MN  circuit 
is  estimated  by  means  of  button  Kf,  which  throws  in  a  shunt  resistor 
of  240  kJl(R?X  The  amplifier  is  resistance-capacity  coupled  and 
is  built  around  subminiature  valves;  the  over-all  gain  is  about 
2000.  The  output  is  connected  to  a  bridge  rectifier,  which  feeds 
the  M-24  meter;  the  scale  is  graduated  0-10  and  5-30,  to  correspond 
with  the  ranges  of  measurement. 

The  calibrator  serves  to  check  the  gain;  it  is  simply  an  LC 
oscillator  built  around  the  P-13  transistor. 

The  instrument  is  fed  from  a  1.48-PMTs-9  cell  (filaments) 
and  a  GB-80  battery  (plates);  the  instrument  works  normally  over 
voltage  ranges  of  1-1.3  (filaments)  and  60-90  (plates).  The  supply 
unit  is  such  as  to  provide  ?0-80  hours'  continuous  use. 

The  millivoltmeter  is  contained  in  a  metal  case  135  x  150  x 
280  mm;  the  weight  with  the  batteries  fitted  is  about  2.3  kg. 

The  apparatus  is  simple  to  make  and  service;  it  is  also  light. 
It  should  be  usable  under  a  variety  of  conditions,  especially  in 
inaccessible  regions.  The  only  disadvantage  is  that  it  cannot  be 
used  when  the  level  of  industrial  interference  is  high. 
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CONCERNING  THE  VERTICAL  T0R8I0N  BALANCE  THEORY 
pp.  278-287  B.  V.  Oran 


1.  Much  attention  has  bean  given  to  the  Measurement  of 
the  vertical  gradient  of  the  gravitational  acceleration  in  recent 
years,  for  the  quantity  is  useful  in  geophysies,  geodesy,  and 
related  fields.  The  need  for  a  special  instrument  to  measure 
Ussr<  arose  long  ago;  the  problea  was  considered  by  E(5tv(5s,  while 
in  19^3  Haal ck  /§/  wrote  that  "the  next  task  of  instrumental 
gravimetry  is  to  produce  an  instrument  for  measuring  the  vertical 
gradient".  In  1920  Berroth  proposed  a  vertical  torsion 
balance  for  measuring  Uzz* }  Sadovskiy  /k/  has  considered  the 
theory  of  such  balances,  but  no  attempt  to  utilize  the  method  has 
been  described,  although  the  proposal  has. been  discussed  more 
than  once  /7/. 

Previous  treatments  deal  only  with  the  case  in  which 

the  center  of  ease  lies  on  the  axis  of  rotation;  the  torques 
produced  by  the  second  derivatives  of  the  potential  are  small. 

A  slight  displacement  of  the  center  of  mass  from  the  axis  would 
produce  a  torque  comparable  with  that  from  the  second  derivatives, 
so  there  is  some  point  in  considering  the  case  in  order  to  find 
the  permissible  limits  of  such  deviation.  Here  I  deduce  the 
equation  of  equilibrium  for  a  body  of  arbitrary  shape  having  a 
horizontal  axis  of  rotation  that  does  not  pass  through  the  center 
of  mass.  This  equation  is  discussed  from  the  point  of  view  of 
possible  absolute  and  relative  measurements  of  Uzx. 


2.  Consider  a  solid  of  arbitrary  shape  having  a  horiz¬ 
ontal  axis  and  suspended  on  two  elastic  fibers.  The  two  coord¬ 
inate  systems  to  be  used  are  (mobile,  rigidly  connected  to 
the  body)  and  xyz  (fixed);  the  origin  of  both  lies  at  0,  the 
point  where  the  axis  of  rotation  meets  a  plane  passing  through 
the  center  of  mass  and  normal  to  that  axis.  The  Oz  axis  is 
directed  vertically  downwards,  Ox  and  Oy  lying  in  a  horizontal 
plane.  The  Oij  axis  is  directed  along  the  axis  of  rotation  and 
so  is  horizontal;  Of  is  normal  to  0^,  and  0%  lies  normal  to  the 


| _ *1  call  Berroth* s  system  a  vertical  torsion  balance. 
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fori 


ees  is 


The  torque  about  0^  produced  by  the  gravitational 

Mi  —  ddm> 


in  which  V  is  the  volume  of  the  body,  £  and  J  are  the  coordinates 
of  an  element  of  mass  dm,  and  g«  and  g.  are  the  components  of  g 
in  the  coordinate  system.  This  M|  is  balanced  by  an  equal 
torque  from  the  fibers  when  the  system  is  at  equilibrium;  if  we 
neglect  the  sag  in  the  fibers,  we  have  that 


(l) 


in  which  is  the  torsion  constant  of  the  fibers  and  is  an 
angle  defining  the  position  of  the  body  when  there  is  no  torque 
in  the  fiber.  Transformation  of  coordinates  gives 


gt  -  Bx  cos  «,  -f  gy  cos  3,  +  gx  cos  Y, , 
Be  -  gx  cos «,  gy  cos  p,  -f  g z  cos  If,. 


(2) 


in  which  a^,  3i,  and  Vi  are  the  angles  between  the  axes  of  the 
two  coordinate  systems.  The  direction  cosines  are  found  from 
the  condidtion  that  the  0$  axis  makes  an  angle  0  with  the  Oz 
axis  and  that  0^  makes  an  angle  f  with  Oy.  Substitution  for 
the  direction  cosines  in  (2)  gives  that 

gt  ’"'B* cos  <?  cos  8  +  gy  sin  s  cos  8  —  gt  sin  8, 
g t  ~gx  cos  <p  sin  8  -f  gy  sin  5  sin  8  +  gz  cos  6.  ( 5) 


We  expand  gx,  gy,  and  gz  as  Taylor  series  on  the  assump¬ 
tion  that  the  second  derivatives  are  constant  within  the  volume 
of  the  body: 


Bx  -  (**)«  +  Uxx‘X  -  Uxy  )f  +  Uxx ■  z. 
By  :m  (gy)*  ■+■  UXy  •  X  +  Uyy  -.V  +  Uyx  •  2 , 
Bf  =(Bi)o  +  Uxz-X  +  U„-y  -f  Utzz. 


(4) 


Here  (gx)e*  (gy)«»  0  and  (gz)»*  go  •  Then  the  coordinate- 
transformation  formulas  give  us  that 

x  «  I  cos?  cos  •  —  vjsiny  -f  Ccosf  sin  8, 
y  —  $Un?cos8-f  ijcos?  -K  sin?  sin  8,  (5> 

*  =  —  (sin  8  -j-  Ceos  8. 

We  substitute  (5)  into  (4),  then  (4)  into  (3),  and  finally  (3) 
into  (1)  to  get  that 


*(*-♦0)=  —  B$ 


i 


(€cos  8  +  C  sin  6)  dm 
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*  r  1* 


_ J*-  - 


+  u,x  f  j[(,.,-«‘Jr!+  Kcoi2*]Co«V 

+  (t<l  tin  •  -  tfcoi  •)  }  dm  +  U„  j  (|<?  -  {*) 

+  ^  cos  2®  1  its*T  —  (?S]  sin  6  —  Yj’  col  I)  dm 

-  U„  f[(P-P)  *-4^+  Kcoi 

+  U„  f  |<C*  -  P)  I42  +  «  COI  2»]  tin  2t 

— (Ci) sin*  —  v col*) cos 2||  dm  -ri4#j{[(C*-l*)co,2| 

—  2(t  sin  28  cos?  +  (fy cos#  +  sln8) slnrjdm 
+  uyxj  [(C*  -  &*)  cos  26  -  2K  tin  28]  sin  ? 

—  (bjcos  8  +  if,  sin  8)cos  tj  dm. 
j(<? -{•)<*■ -AT.  \yfMm=J„ 


Now  by  definition 


^Uim  =  mU,  mt,, 


in  which  m  is  the  mass  of  the  body  and  and. Jo  are  the  coord¬ 

inates  of  the  center  of  mass.  Further,  Laplace's  equation  gives 

Vyy  «■  —  UXX  —  Un.<  (9) 

Then  (7)-(9)  enable  us  to  transform  (6)  to  give 

*(♦  “ ’V«)  -  -  *»*♦($,  cos8  +  C*sln8) - 
“~i  ux  [(ys,n  28 +  y,  cos  28)  cos  2?  -f  (A  sin  8  —  7,  cos  6)  sin  2?  j 

+  Uxy  [(y  *to  2*  +  /,  cos  28)  sin  2?  —  (7t  sin  8  —  7,  cos  8)  cos  2* j  (io) 

+  Uxm  [(* cos  28-27,  sin  28)  cos  ?  +  (7,  cos  6  +  7,  sin  8)  sin  ? 

~~  +  Uyu  |(^c  os  28—  27,  sin  28)  slh  ?  —  (7,  cos  8  +  7,  sin  8)  cos?  ~ 
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“  7  »in26  -f-  /,cos26j .  (1Qj 

U k  ■»  U yy  Uxx. 

Then  (10)  is  the  equation  of 
equilibrium  for  a  body  of  arbitrary 
shape  having  a  horizontal  axis  of 
rotation.  For  a  vertical  torsion 
balance  (Fig.  1),  the  symmetry  about 
the  and  *)0£  planes  makes  the 
central  moments  Ji  ,  Jj  ,  and  Jj  zero, 

^  so  (10)  becomes 

—  1»o)  =  “  mg*  (50  co*  6  +  ^  sin  6)  —  A  sin  26 

—  j  UkK  sin  26  cos  2?  4-  ^'UxyK  sin  26  sin  2?  ( 1 1 ) 

+  UXgK  cos  26  cos  ?  -}-  UygK  cos  26  sin  9. 

The  main  difference  between  earlier  equations  and 

(11)  lies  in  the  first  term  on  the  right;  this  term  vanishes 
if  0.  Other  differences  are  purely  formal,  for  prev¬ 
iously  the  choice  of  the  mobile  coordinate  system  was  not 

correct,  in  that  the  0£  axis  was  directed  vertically  downwards, 
so  this  system  was  not  rigidly  linked  to  the  body  and,  moreover, 
was  not  rectangular. 

Equation  (11)  contains  the  seven  unknowns  «|,o  ,  go  ,  Ufl,  U  , 
Uxy  «xz*  and  UyZ,  so  the  balance  senses  g  as  well  as  all  the 
second  derivatives.  Now  Uzz  contains  only  the  angle  6;  let  the 
observations  be  made  in  a  way  such  that  0  remains  constant  at  any 
given  point,  so  only  cfi  varies.  Then  the  first  two  terms  on  the 
right  in  (11)  also  remain  constant,  so  they  determine  merely  the 
zero  reading  C*  : 


co-~  *!»•  + W*oft>cos0  -f  Co  sin  0)  +  A  UttKi\nW.  (12) 

This  cannot  give  us  Uzz,  for  (12)  contains  f'o  and  go  as 
unknowns  in  addition;  thus  it  is  impossible  to  keep  6  constant 
in  absolute  determinations  of  Uzz,  and  observations  must  be  made 
for  various  values  of  &. 

Now  we  consider  when  the  first  term  in  (11)  can  be  negl¬ 
ected;  we  must  have  that 

I  **£o(§oCosO  ;0sln  6)|  <  —  U„K sin  26 1 ,  (13) 


in  which  Uzz  is  the  minimum  value  of  the  vertical  gradient  that 
has  to  be  measured.  If  the  balance  arm  is  weightless,  we  have 
I  from  (7)  that 

J —  K  =  —  ml*.  .  n  41. 


(14J-J 


135 


in  which  Z  is  the  distance  from  the  axis  of  rotation  to  the  load. 
It  can  be  shown  that  0  must  be  45®  if  the  sensitivity  to  tlzz  under 
these  conditions  is  to  be  maximal;  then  (13)  can  be  put  as 


Let  the  desired  limit  be  10  E  (10"8cgs  units); 


gives  us  that 

15*  +  Col  < /MO-11  cm. 


(15> 

then  (15) 

(16) 


Now  a  balance  arm  of  total  length  50  cm  implies  that  If#  +  To  I  i 
0.66  x  10”0 cm ;  (16)  also  shows  that  the  only  way  of  increasing 

the  permissible  deviation  from  the  axis  of  rotation  is  to  increase 
l,  but  the  effects  of  temperature  and  of  deflection  in  the  arm 
increase  rapidly  with  Z.  The  system  is  unworkable  for  X  *  50 
cm,  although  in  this  case  ]Jo  ♦  jo |  should  be  less  than  2. 5  x  10~* 
This  shows  that  the  axis  of  rotation  must  coincide  very  precisely 
with  the  center  of  mass  for  any  usable  value  of  [ ;  the  desired 
state  is  not  at  present  attainable.  That  is,  a  vertical  torsion 
balance  is  always  a  gravimeter  as  well  as  gradient  meter  no  matter 
how  carefully  it  may  be  adjusted. 


3.  Consider  now  a  balance  with  real  values  of  Jo  and  Jo • 

We  employ  the  usual  technique  to  estimate  the  probable  distance 
to  the  axis;  an  analytic  balance  is  a  vertical  torsion  balancfe 
if  we  take  the  masses  of  the  pans  and  weights  to  be  concentrated 
at  the  points  of  support.  Morozova  has  shown  that  the  dist¬ 
ance  from  the  axis  of  rotation  to  the  center  of  mass  of  the  beam 
must  not  be  less  than  100  if  the  balance  is  to  work  stably  for 

loads  of  20  to  200  grams;  the  distance  is  sometimes  reduced  to 
10  fi,  but  any  further  reduction  leads  to  great  instability.  It 
would  seem  possible  that  especial  care  under  exceptionally  good 
conditions  might  enable  one  to  make  Jo  =  5o  =  1  ^5  then  (15) 
gives  us  that  this  deviation  of  l  u,  with  2  =  50  cm,  gives  a 

torque  of  3*4  x  10'*  E.  N  w  (11)  shows  that  this  torque  is  reduced 
by  a  factor  of  -fit  if  6  is  reduced  from  45®  to  0®  ,  so  the  torque 
would  then  be  only  1.4  x  10* E.  Even  so,  the  balance  is  subject 
to  a  large  torque,  and  this  varieG  with  0;  t..e  small  quantity  can 
be  measured  in  the  presence  of  the  large  one  only  if  the  latter 
is  constant,  the  gravimeter  being  an  example  of  this.  That  is, 
the  first  term  in  (11)  must  remain  constant  at  each  point  if  the 
Uzz  are  to  be  measured  to  5-10  E,  so  0  must  not  alter;  but  then 
the  absolute  value  cannot  be  determined,  so  a  vertical  torsion 
balance  cannot  be  used  to  measure  Uzz  absolutely.  This  conclus¬ 
ion  is  a  result  of  purely  practical  considerations,  not  of  theory; 
it  is  merely  that  at  present  there  is  no  sufficiently  rigid  mat- 


,  erial  for  the  arm  and  no  adequate  mean  of  bringing  the  center  of 
j _ mass  close  enough  to  the  axis  of  rotation.  I 
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4,  Consider  now  the  relative  measurement  of  UZB  from  (nT1 
with  6  constant;  (12)  reduces  the  unknowns  in  (11)  to  five,  and  ! 
the  coefficients  to  these  are  functions  of  f  (apart  from  Co),  so 
readings  .in  five  azimuths  give  us  a  system  of  five  equations  to 
derive  Co*  Let  the  values  for  two  adjacent  points  be  Coi  and  Co*  ; 
then,  if  jf'o  ,  go  ,  J0  ,  and  Jo  are  constant  for  these  points,  we 
have  from  (12)  that  -  r 

U" 26  ’  U7> 

That  is,  relative  measurements  of  Uzz  are  possible  if  the  above 
assumptions  are  obeyed  and  if  $  varies  only  slowly  with  Uzz  (i.e., 
remains  virtually  constant).  Now  we  have  to  consider  the  permiss¬ 
ible  changes  in  ,  go  «  Je  •  and  Jo  as  between  points.  The  theory 
of  errors  gives  us  the  formula 


in  which  df* is  the  dispersion  of  a  function  whose  arguments  have 
dispersions  dx?  ,  dy*,  ...  .  Let  the  of  (11)  be  the  function  and 
6,  together  with  the  above  four  quantities,  the  arguments;  the 
changes  in  are  small,  so  we  neglect  the  last  five  terms  on  the 
right.  Then  (11)  and  (18)  together  give  us  that 

xW  =  t'dti  +  («5o  \ '  Wdgo  +  ( 19 ) 

On  the  assumptions  that  all  terms  contribute  equally,  we  have 

d%  —  ^, 

A  —  ^ 

(20) 

Consider  Vm;  the  equation  for  the  the  angular  sensitivity  of 
a  vertical  torsion  balnce  takes  the  form* 

M  **sln2#  (21) 

4Utt  x  +  mg0  cos  i~i,  sin  •) 

Then  (21)  becomes  in  conjunction  with  (14)  for  $ s  45°  and  50 =  Jo 


_  ImC 

dUu  T* 


We  can  derive  (21)  from  (11)  by  means  of  tne  rule  for  finding 
the  derivative  of  an  inexplicit  function  if,  in  view  of  the  very 
j  small  changes  in  t^*,  we  neglect  changes  in  the  second  derivative^. 
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j  Then 


• = 3  *U„  n 

m  4  ’  ttl  ‘ 1  ‘ 


W^-  insert  (23)  into  (20)  to  get 


* o 


/T  =  •1 1  2 

0  4  *o  ’ 


(23) 


(24) 


The  differentials  nay  be  found  by  putting  dd/dUzz  ■  1"/1  E{  dUzz 
•  10  E,  d^  *  10",  21  a  50  cm,  and  jo »  5« a  cm,  so  (24)  gives 

*♦.  =  5', 

dg0—  16.5  myj.1,  (25) 

du  0.33- 10 '•  CM. 

(25)  shows  that  the  position  of  the  center  of  mass  should  not 
vary  more  than  O.33  x  10“®cm  as  between  points;  this  distance 
is  of  the  order  of  the  dimensions  of  atoms.  Long  ago  Mendeleyev 
/!/  remarked  that  the  position  of  the  center  of  gravity  of  a 
precision  balance  relative  to  the  axi3  of  rotation  alters  with 
time  even  if  the  pans  are  allowed  to  swing  freely;  one  of  the 
main  causes  of  this  is  that  temperature  fluctuations  occur.  Also, 
the  ceneter  of  gravity  is  appreciably  affected  if  the  balance  is 
brought  to  rest.  The  examples  Mendeleyev  gives  show  that  the 
center  of  gravity  of  a  precision  balance  can  alter  by  10”4  cm  or 
more,  so  we  may  expect  that  the  jo  and  Jo  of  a  vertical  torsion 
balance  will  vary  by  far  more  than  the  above  limit.  Further, 

(24)  shows  that  the  sole  means  of  increasing  the  permitted  range 
in  jo  and  5# is  to  make  7  large;  but  7  must  be  increased  tenfold 
in  order  to  increase  jo  and  ?o  by  factors  of  100,  which  is  quite 
impossible,  for  the  total  length  of  the  beam  would  then  be  5  »• 
Thus  the  very  rigorous  demand  for  stability  in  the  center  of  mass 
makes  it  doubtful  whether  Uzz  can  be  measured  even  in  relative 
terms. 

I  do  not  propose  to  deal  with  the  equations  of  motion  and 
sensitivity,  for  (11)  alone  gives  a  sufficiently  clear  indication 
of  the  possibilities  of  absolute  and  relative  measurement  of  Uzz. 

Conclusions 

The  equation  of  equilibrium  for  a  vertical  torsion  balance 
has  been  examined;  the  following  deductions  have  been  drawn. 

1.  The  displacement  of  the  center  of  mass  from  the  axis 
of  rotation  cannot  be  neglected;  a  displacement  of  10”* cm  causes 
a  torque  that  is  several  orders  of  magnitude  larger  than  that  1 
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•rifling  from  the  second  derivatives,  and  this  torque  1*  not  even  ; 
constant  during  absolute  measurements  of  U2Z.  That  is,  the 
instrument  actually  acts  mainly  as  a  gravimeter,  and  it  is  imposs¬ 
ible  to  measure  gradients  of  the  order  of  10  E. 

2.  Relative  measurements  on  U22  could  be  made  if  the  tilt 
of  the  beam  were  to  remain  constant  when  the  balance  is  moved;  but 
the  permissible  displacement  of  the  center  of  mass  is  then  less 
than  1  a.  Experience  with  precision  balances  shows  that  such  stab¬ 
ility  is  at  present  unattainable,  so  it  is  doubtful  whether  the 
instrument  can  make  relative  measurements  of  U  . 

This  means  that  it  would  be  undesirable  to  devise  some  form 
of  vertical  torsion  balance  with  the  object  of  measuring  U22.  The 
discrepancy  between  the  requirements  of  theory  and  the  practical 
possibilities  is  so  great  that  it  would  be  better  to  seek  other 
and  more  promising  methods.  I  consider  that  measurements  with 
gravimeters  at  different  heights  and  Jolly's  method  are  suitable 
for  this  purpose. 

I  am  Indebted  to  Professor  5.  A.  Andreyev  and  Dr.  V.  S. 
Mironov  for  valuable  discussions. 
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PLOTTING  WAVE-FRONT  CHARTS  WITH  AN  ARBITRARY  LAW  OF 
CHANGE  IN  SEISMIC  VELOCITY  WITH  DEPTH 
pp.  292-300  Yu.  Ye.  Nekrasov 

Ray  diagram?  are  u.;ed  in  interpreting  seisaic  observations 
if  the  distribution  of  the  seismic  velocity  in  the  section  is 
complicated,  for  the  usual  mean-velocity  method  then  becomes  inap¬ 
plicable.  Methods  have  been  published  /l,3-fj7  for  calculating  ray 
diagrams  for  given  laws  of  variation  of  velocity  with  depth;  but 
the  laws  must  be  available  in  analytic  form  %  whereas  many 

practical  sections  have  velocity  variations  too  complicated  for 
such  methods  to  be  applicable. 


Fig.  1. 

Here  I  give  a  method  that  can  be  used  for  any  graphically 
specified  law  of  variation;  it  consists  in  approximate  calculat¬ 
ion  of  the  Integrals 
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r 


(2) 


_ dx _ 

v{t)  /T— »M(2) 


Here  (1)  is  ths  aquation  of  a  ray,  (2)  being  the  time  of  travel 
along  mray  fe] ;  x  and  z  are  the  coordinates  of  points  on  the  ray, 
v(z)  is  the  velocity,  and  a  ■  1/v*  is  the  ray  parameter.  Consider 
now  the  construction  of  rays;  here  we  use  the  accessory  function 


10 

—  i Hi* 


If  v(z)  is  given  in  graph  form,  we  can  construct  Fj(z).  For  ths 
appropriate  depths  we  read  off  v  from  the  graph  of  v(z),  and 
from  a  table**read  off  the  value  of  Fi  corresponding  to  the  given 
v* ;  this  value  is  entered  on  the  graph  of  Ft(z)  at  the  appropr¬ 
iate  t  (Fig.  la).  A  discontinuity  in  the  velocity,  as  at  z?  gives 
us  two  values  of  F«  (z).  The  points  are  joined  up  into  a  smooth 
curve;  Fig.  la  shows  two  curves  for  two  values  of  v*.  The  next 
step  is  to  integrate  Fi(z),  which  is  best  done  by  the  trapezium 
method  as  applied  to  graphs  drawn  on  mm  graph  paper  with  depth 


scales  or  1:1000  and  so  on.  The  integration  amounts  to  the  det¬ 
ermination  of  the  areas 

Sk  =  flk^zh 


in  which  z  ■  -  *k  (Fig.  lb);  if  now  z  is  100  a  (for  a 

scale  of  1:10  000;  10  m  for  1:1000),  the  area  is  found  simply  by 

sultiplying  the  ordinate  hfc  by  100  (or  10).  Now  is  numerically 
equal  to  the  increase  in  the  horizontal  path  of  the  ray  for  the 
change  in  depth  of  Zfc+i  -  z^;  the  sum  of  the  is  the  horiz¬ 
ontal  path  consequent  on  penetration  to  the  depth  zn,  so  the  prob¬ 
lem  of  constructing  the  ray  is  solved. 

Consider  now  the  time  equation: 


f _ _ 

J  V  (*)  ■/ 1  —  irV * 


(*) 


*K.  A.  Nekrasova  performed  the  calculations. 


(2) 
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f*  -50 


41) 


1.5 

1.6 
•  7 
1.8 

1.9 
2.0 
2.1 
24 
22 

2.4 
24 
24 
2.7 
23 

2.9 
3.0 

3.1 

3.2 
33 

3.4 
33 

3.6 

3.7 
33 

3.9 
4.0 

4.1 

4.2 
43 

4.4 
43 
43 

4.7 
43 

4.9 
53 

5.1 

5.2 
53 

5.4 
53 

5.6 

5.7 
53 

5.9 

6,0 


0,030 

0,032 

0,034 

0.036 
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0,040 
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0.046 
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0,058 
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0,066 
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0.074 
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0,060 

0,062 

0364 

0,086 

0,068 

0390 

0392 

0,094 

0.096 

0,096 

0,100 

0,102 

0,104 

0,106 

0,106 

0.110 

0,112 

0.114 

0,116 

0,118 

0,120 


f, 

0,6675 

0,6257 

0,5889 

0,5662 

03269 

03005 

0,4766 

0,4550 

0,4352 

0,4170 

0,4003 

03851 

0,3706 

0,3579 

03455 

0,3340 

0,3232 

0,3131 

03087 

0.2947 

0,2863 

0,2786 

0.2711 

0,2639 

0,2572 

0,2508 

0,2446 

0,2388 

0,2335 

0,2282 

0,2231 

0,2183 

03136 

0,2694 

03051 

0,2010 

0,1968 

0,1933 

0,1896 

0,1889 

0,1829 

0,1797 

0,1?65 

0,1786 

0,1707 

0,1678 


v”  -  30 


Ft 

0,050 
0,053 
0,057 
0,060 
0,063 
0,067 
0,070 
0,073 
0,077 
0,080 
0,083 
0.087 
0,090 
0,094 
0,097 
0,100 
0,103 
0,107 
0,110 
0,113 
0.1 .6 
0,121 
0,124 

0,128 

0,131 

0,134 

0,137 

0,141 

0,145 

0.149 

0,152 

0,156 

0,159 

0,162 

0,165 

0,170 

0,173 

0,176 

0,180 

0,183 

0,186 

0,190 

0,193 

0,196 

0,201 

0,204 


20 


0,6675 

0,6258 

0,6893 

0,5568 

0,5274 

0,5010 

0,4771 

0,4560 

0,4361 

0,4179 

0,4016 

0,3861 

0,3719 

0,3585 

0,3465 

0,3350 

0,3242 

0,3144 

0,3049 

0.2959 

0.2877 

0,2797 

0,2725 

0.2653 

0.2585 

0.2820 

0,2461 

0,2406 

0,2349 

0.2298 

0,2247 

0,2200 

0,2153 

0,2111 

0,2070 

0,2028 

0,1990 

0,1952 

0,1916 

0,1802 

0,1850 

0,1819 

0,1787 

0,1758 

0,1730 

0,1701 


Ft 

0,075 

0,080 

0,086 

0,090 

0,095 

0,100 

0,106 

0,110 

0,116 

0,121 

0,120 

0,131 

0,136 

0,141 

0,146 

0,151 

0156 

0,161 

0,168 

0,173 

0,178 

0,183 

0,188 

0,193 

0,199 

0,205 

0.209 

0,214 

0,220 

0,226 

0,231 

0,236 

0,242 

0.247 

0,252 

0.258 

0,264 

0,269 

0,275 

0.280 

0,286 

0,292 

0,297 

0,303 

0,306 

0,315 


«*-  15 


0,6684 

0,6270 

0.5907 

0,5577 

0,5288 

0,5025 

0,4792 

0.4672 

0,4378 

0,4196 

0.4032 

0,3876 

04737 

04607 

0.3483 

0,3370 

0,3265 

04167 

0,3073 

0,2966 

Q.2904 

0,2823 

04750 

04680 

9,2614 

0.2551 

0,2491 

0.2434 

0,2382 

0,2331 

0,2282 

0.2234 

0,2189 

0,2146 

0,2104 

04066 

0,2028 

0,1991 

9,1957 

0,1923 

0,1892 

0.1860 

0,1831 

0,1802 

0,1773 

0,1747 


Ft 

0,100 

0,106 

0,113 

0,121 

0,128 

0,134 

0,141 

0,149 

0,154 

0,159 

0,170 

0,176 

0,183 

0,189 

0,196 

0,204 

0410 

0,218 

0,226 

0,233 

0,240 

0,247 

0,254 

0,263 

0,260 
0,276 
0485 
049> 
0,298 
0407 
0418 
0422 
0,330 
0,338 
0446 
0455 
0462 
0470 
0,379 
0,386 
0 M6 
0,403 
0,411 
0,419 
0,429 
0,434 


0,6702 

04289 

0,6017 

0,5596 

0,5305 

0.5045 

0,4810 

0,4596 

0,4401 

0.4218 

0,4057 

0,3905 

04764 

04631 

04615 

0,3401 

04295 

0.3199 

0,3108 

04019 

0,2939 

0,2860 

0.2769 

0,2721 

0.2655 

0,2503 

04536 

0.2480 

0,2428 

04378 

0,2329 

04284 

0,2240 

04200 

04160 

0,2123 

0.2086 

04050 

04018 

0,1965 

0,1955 

0,1924 

0,1896 

0,1870 

0,1844 

0,1820 


Ft 

0,126 

0,134 

0,|43 

0,152 

0.158 

0.170 

0.178 

0,186 

0,195 

0404 

0413 

0422 

0.231 

0.240 

0.249 

0,258 

0467 

0,277 

0,286 

0496 

0,306 

0.315 

0424 

0,334 

0444 

0455 

0464 

0474 

0483 

0494 

0,404 

0416 

0.426 

0,437 

0447 

0457 

0,470 

0,480 

0,493 

0404 

0416 

0428 

0440 

0462 

0463 

0478 
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0,6720 

0,6305 

04042 

04618 

04333 

0407» 

0,4840 

0,4623 

0,4403 

0,4252 

0,4090 

04940 

04802 

04674 

04556 

04443 

0.1339 

04241 

04164 

04067 

0,2960 

04012 

04642 

04773 

04711 

04654 

04505 

04641 

04490 

0,2444 

04397 

04865 

0.2313 

04274 

04235 

04200 

0,2166 

04136 

04104 

04074 

0,2045 

04020 

0,1994 

0.19hS 

0.1946 

0.192$ 


f(#) 

i  v*  10,0 

\§MS32EIM 

v*  -  7,5 

v*  •  7,0 

DD 

Fi 

■a 

F, 

Ft 

Ft 

Ft  | 

1  Ft 

14 

14 

1,7 

0,152 

0,162 

0,173 

04739 

04333 

04074 

0,170 

0,181 

0,192 

ill 

0,190 

0406 

0417 

ooo 

§13 

—  —  40 

0.204 

0,219 

0,233 

0,6003 

04398 

0,0039 

0419 

0434 

0,250 

04821 

0,6414 

04069 
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tf(*) 

V* 

6.5 

V* 

6.0 

55 

v * 

■■  5,0 

»• 

4.8 

I 

Fx 

Ft 

Fx 

1  * 

F 

Fx 

Ft 

Ft 

F> 

1.9 

0.306 

0,5607 

0,335 

05553 

0.368 

0.5612 

0,411 

05721 

0,431 

0.5734 

2.0 

0,324 

0,6258 

0556 

05301 

EK21 

ril 

0,437 

05456 

0.467 

0,5301 

2,1 

0,342 

05033 

0,374 

05081 

0,414 

05155 

0,462 

05247 

0586 

05319 

2,2 

0,362 

0,4636 

0,395 

0,4888 

031 

0,4963 

6,400 

05061 

0516 

05115 

2,3 

0,379 

0,4651 

0,416 

0.4710 

0,460 

0.4769 

0516 

0,4397 

0546 

0,4965 

2,4 

0,396 

0,4464 

0,437 

0,4550 

0,484 

0546 

0578 

0,4612 

2.5 

0,417 

0,4333 

0,459 

0,4401 

0509 

0578 

ts 

0510 

0,4660 

2.6 

0.437 

0,4196 

0,482 

0,4270 

0537 

0500 

0,848 

0.4602 

2.7 

0,459 

0,4076 

0,604 

0,4148 

0565 

■Jr-t 

0541 

Pftv 

0570 

0,445ft 

2,6 

0,479 

0,3959 

05?3 

0,4040 

0591 

sir 

0575 

CEtt 

0.717 

0,4390 

2.9 

0,496 

05852 

0,563 

0,3940 

0521 

0,712 

0,759 

0,4310 

3,0 

0521 

05758 

0,578 

05849 

0,649 

h  rt 

0,751 

0500 

0,4270 

3,1 

0521 

05758 

0578 

0,3649 

0,683 

life  mV 

i  'JL, 

0546 

0,4223 

3.2 

0567 

0,3502 

0531 

0,3696 

0,715 

32 

rrr- 

0596 

0,4196 

3.3 

0590 

05520 

0561 

05634 

0,751 

05766 

0.4039 

0,948 

0,4196 

3.4 

•516 

0,3457 

0589 

05669 

0,786 

rwi 

£ 

0.4010 

150' 

0,4167 

3.5 

0542 

0.33 CT 

0,720 

0,3630 

0527 

O 

0.4009 

0.4175 

3.6 

0565 

05334 

0.751 

05472 

0567 

05674 

rr 

0.4000 

3.7 

0,694 

05292 

0.784 

05434 

0,909 

:  2 

0u46l«t 

1,212 

0.723 

05249 

0,820 

0,3406 

0,966 

EC!!1 

t*? 

0,4049 

1.279 

05314 

3J 

0.754 

05213 

0.656 

05374 

1,006 

05636 

& 

0,4097 

1592 

05396 

4.0 

0,781 

05173 

0.896 

05360 

1,063 

V  ^  . 

05:67 

1507 

0,4521 

4.1 

0,614 

05143 

0.935 

05341 

1,118 

nr 

0.4274 

1542 

0,4699 

4.2 

0546 

05116 

4981 

05334 

1,185 

>v  vm  l 

frf 

0.4690 

0,4928 

4.3 

0,884 

05105 

1,026 

05332 

1,260 

65556 

1019 

05280 

4,4 

0,922 

0,3062 

I.C80 

0.3347 

1,335 

'CX. 

0,4675 

2563 

05669 

4,5 

0,962 

0.3063 

1134 

0,3362 

1,428 

0,3672 

05099 

1668 

05313 

4.6 

I.OOJ 

05079 

1.196 

''1386 

1530 

tftt 

05556 

3541 

0,7605 

4.7 

1,046 

0,8079 

•563 

05496 

1,643 

2.766 

05&» 

4.790 

KM 

4.8 

1598 

05096 

1544 

05489 

1,792 

0,4277 

3,426 

0,7463 

4.9 

1.152 

05115 

1,419 

0.8544 

1,974 

0,4616 

4,964 

1.0363 

5,0 

1,907 

05136 

1507 

05017 

2,176 

0.4716 

5.1 

1.206 

05166 

1514 

05729 

2,491 

05272 

5.2 

1,336 

05210 

1.739 

95961 

2593 

I\JS  Vr  1 

5.3 

1,413 

05265 

1594 

0,4040 

BCSl 

5.4 

1,407 

05337 

2566 

0,4246 

5550 

1,0060 

5.5 

1503 

05486 

2,301 

0,47r 

5.6 

1,700 

05529 

2511 

0.4966 

5.7 

1527 

05665 

3,060 

05694 

5.8 

1564 

05631 

3571 

05697 

5.9 

3,170 

6.4066 

0.9264 

6,0 

2.400 

0,4320 

«**) 

v*  —  4,6 

V* 

-4,9 

0 * 

-4.0 

0  *  -  3.8 

F, 

F, 

F , 

Fx 

Ft 

Fx 

Ft 

Fx 

Ft 

1.5 

0545 

0,7057 

0503 

0,7092 

0,369 

6,7136 

0,404 

0,7094 

0,490 

0,7257 

15 

0,371 

0,6667 

0591 

0,6707 

0,412 

05766 

0.437 

05621 

0.461 

05682 

1.7 

0596 

0,6329 

0,418 

0,6369 

0,443 

05435 

0,470 

05602 

0.491 

0546ft 

1,8 

0,426 

0,6039 

0,448 

0,6086 

0,474 

05154 

0504 

05202 

0530 

05813 

1.9 

0,454 

0,5807 

0,479 

05638 

0,508 

0,0116 

0,540 

05961 

0571 

0,6079 
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v*  -  4,2 


v *- 


0,5656 
0,5353 
0,6173 
05018 
0,4883 
0.4768 
0,4664 
0,4570 
0,4501 
0,4443 
0,4400 
0,4371 
0,4348 
0,4354 
0,4369 
0,4403 
0,4460 
[0,4541 
0,4669 
1 0,4850 
0,5061 
05871 
05831 
0,6515 
0,7762 
1,0684 


05609  0,541 
0.5420  0578 
0,5249  0,616 
0,6097  0,654 
0,4975  0,697 
0,4859  0.740 
0,4766  0,789 
0.4695  0,839 
0,4625  0.H96 
0,4587  0,956 
0.4560  1,019 
0.4541  1,096 
0,4550  1,176 
0,4587  1,270 
0,4638  1376 
0,4715  1,497 
0.4673  1,662 
0,4989  1.862 
05222  2.131 
05540  2513 
05099  3,118 
0.6725  4.584 
0,7937! 

1,0029 


0,5688 

0,5498 

05342 

05195 

0,5081 

0,4975 

0,4900 

0,4836 

0,4801 

0,4789 

0,4762 

0,4394 

0,4822 

0,4904 

0,5010 

05139 

05382 

05714 

05836 

0,6930 

OJI70 

1.1455 


V(X) 

V* 

2.8 

~  2,7 

L’* 

2,6 

"  2,5 

V 0  = 

-  2.4 

F{ 

* 

Ft 

F, 

F. 

Fa 

Ft 

1.5 

0,635 

0,7899 

0,670 

0,8026 

0,706 

0,8157 

0,751 

0,8331 

0,805 

0,8562 

1.6 

0,697 

0,7622 

0,737 

0,7764 

0.784 

0.7843 

0,832 

0,8130 

0,895 

83389 

1.7 

0.763 

0,7396 

0,812 

0,7582 

0,865 

0,7782 

0,928 

0,8026 

1,003 

0,8333 

1.8 

0,837 

0,7241 

0.888 

0,7429 

0,962 

0,7704 

1,038 

0.8006 

1,134 

0,8403 

1.9 

0,922 

0,7163 

0,991 

0,7413 

1,067 

0.7692 

1,169 

0,8097 

1.297 

0.8628 

2.0 

1,019 

0,7123 

1,104 

0,7452 

1,212 

0.7862 

1,335 

0.8347 

1,507 

0,9042 

2.1 

1,136 

0,7205 

1,237 

0,7870 

1,366 

0,8058 

1,547 

0,8772 

1,806 

0.9843 

2.2 

1,276 

0,7369 

1,407 

0,7849 

1  r596 

03662 

1,853 

0,9569 

2.283 

1,1338 

2.3 

1,444 

0.7639 

1,628 

03313 

1,903 

0,9346 

2,347 

1,1468 

3,428 

13528 

2.4 

1,646 

0,8091 

1,955 

0.9158 

7,420 

1,0905 

3,428 

1,4881 

2.5 

1,984 

0.8889 

2.434 

1,0526 

3,480 

1,4493 

2.6 

2,527 

1,0428 

3,630 

1,4451 

2.7 

3,630 

1,3928 

1  l"  = 

=  2,3 

v *  = 

=  2.2 

|  w”-  —  2,1 

r*  - 

23 

F,  j 

Ft 

F  i 

1 

Fx 

1 

Ft 

Fx 

1  f-t 

1,5 

1  0,863 

0,8803 

0,935 

0,9132 

1,017 

0,9506 

1,134 

1,0081 

1,6 

|  0,972 

0.8718 

1,059 

0,9099 

1,176 

0,9643 

1335 

1.043% 

1.7 

1,100 

[  0,8741 

1,222 

0,9294 

1,392 

13091 

1,614 

1,1161 

1.8 

1,258 

i  0,8929 

1,422 

0,9852 

1362 

1,0784 

2,066 

13739 

1.9 

1.472 

1  0,9363 

1,723 

13482 

2.125 

1,2361 

3,060 

13863 

2.0 

2.1 

•2.2 

1,764 

2.240 

3,293 

1 

1,1377 

1,1669 

1.5674 

2,193 

3,224 

1,2048 

1,6077 

3.1 13 

1,6340 
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[The  table  gives  values  of 


1 


F2  ( v )  - 


V  'Z 1  —  «*V* 


as  a  function  of  v;  the  t(z)  curve  is  needed  in  order  to  construct 
isochronous  curves,  and  this  is  found  by  integrating  F* (z)  as  for 
F|(e)  by  aeans  of  tabulated  values  of  F».  Figure  2a  shows  this 
function  in  terns  of  z.  The  problem  is  then  merely  one  of  graph¬ 
ical  integration,  which  in  this  case  it  is  convenient  to  perforn 
with  scales  of  1:10  000  (or  1:1000)  for  depth  and  1  msec/m  per  10 
on  for  F*  .  Here  it  is  also  convenient  to  take  an  interval  in  z 
of  100  n  (or  10  a)|  then  the  area  Q  ■  qn  z  is  readily  found  by 
multiplying  qn  by  100  (or  10),  while  the  sum  of  the  is  t(z)  for 
the  path  of  parameter  v*  down  to  a  depth  of  zn.  All  graphs  should 
be  drawn  on  one  piece  of  paper  (Fig.  2b)  in  order  to  costruct  the 
isochronous  curves;  on  this  we  draw  a  horizontal  line  with  tfc  » 
constant  to  meet  the  t(z)  curves  for  each  ray  with  the  parameter 
v*  at  points  z^,  which  correspond  to  this  time  t^  of  travel  .along 
the  rays.  The  rays  already  constructed  for  the  depths  z^  are 
marked  with  points  corresponding  to  this  t^;  these  points  are 
joined  up  into  a  smooth  isochronous  curve. 

The  table  gives  Fi  (v)  and  F* (v)  for  thirteen  values  of  v* ; 
now  (1)  and  (2)  contain  v(z)/v*,  so  the  table  of  Ff  can  be  used  to 
construct  curves  of  Fi  (z)  for  rays  with  parameters  nv(z)/v*,  but 
the  v(z)  in  the  table  must  then  be  reduced  by  factors  of  n.  Again, 
the  table  of  Fj (v)  can  be  used  to  find  values  of  F{  for  velocities 
Vi  not  listed. 

We  can  put  Fj(v)  in  the  form 


F2(v) 


n  tr 


V**  —  V- 


nv  / n- 1>*!  —  n*r* 


so  for  the  vi  we  find  the  Ft(v)/n  on  the  nv*  curves  opposite  v(z)  ■ 
nv.  Thus  the  table  of  Vi  is  suitable  for  constructing  the  rays  and 
curves  for  sections  having  any  velocities.  The  calculations  are 
no  more  troublesome  than  those  in  the  method  based  on  representing 
v(z)  in  the  form  of  a  piecewise-linear  curve 


1. 


2. 


]_ 
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